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1. INTRODUCTION 


The need for better and more efficient gas turbines requires the availability of cheap 
and reliable design tools for blades used in compressors and turbines. Computational 
methods are the preferred choice for such a design tool, considering the cost and com- 
plexity of wind tunnel experiments. 

Among the computational methods available today, the logical choice seems to be 
a computer code that can solve directly the full Navier Stokes equations. However, given 
the state of the art in both algorithms and computer hardware, such Navier Stokes 
solvers are restricted only to supercomputers, and even then the computation time is 
quite long. 

In order to enable fast and efficient computations, the viscous inviscid interaction 
code was developed by Cebeci [Ref. 1]. The approach used in this code is to solve the 
outer flow field using potential methods, and solving the boundary layer flow using a 
boundary layer method subject to an interaction law, that couples the inner and the 
outer flows. This interaction law 1s needed because classical boundary laver methods fail 
in areas of flow reversal and separation, which are very common in real life flows. 

The viscous inviscid interaction code was originally developed, and successfully used, 
for flows about single airfoils. It was later adapted to cascade flows. 

[n this thesis the applicability of the code to cascades was investigated bv comparing 
its results to experimental data. It was found that although the code can reasonably 
predict experimental results in some cases, it still needs improvements before it can be 
applied generallv as a reliable design tool. 

A major restriction in improving the code is the lack of a wide data base of appro- 
priate experimental results. Some of the key elements in the code, like transition and 
turbulence modelling, are based on empirical correlations, and more detailed and accu- 
rate experiments should be performed, before a better understanding of these phenom- 
gnascan be achieved. 

In the following. the theoretical background of the code 15 presented in Chapter IT, 
a description of the code in Chapter III, the results and discussions are presented in 
Chapter IV and the conclusions and recommendations in Chapter V. A listing of the 


computer program is given in Appendix A. 


I. CASCADE FLOW PROBLEM FORMULATION 


This chapter outlines the theoretical background of the viscous inviscid interactive 
method used in the computer code to investigate cascade flows. Only the major steps in 
the mathematical developments will be described here. A detailed description of the 
theory and the numerical methods 1s given by Cebeci and Bradshaw [Ref. 2] and by 


Krainer (Ref. 3] on which this chapter 1s based. 


A. INVISCID FLOW MEIHOD 
Inviscid flow is the first building block of the flow and 1s solved using the panel 
method. The incompressible two dimensional outer flow must satisfy the Laplace 


equation: 
2 
УФ = 0, 


subject to the boundary conditions on the surface of the blade: 


Where the commonly used boundary condition of zero normal velocity on the surface 1s 
replaced by a specified blowing velocity v, to allow for the effect of the boundary ат 
on the outer flow. 

In addition, the Kutta condition must be satisfied, in order to prevent the existence 
of discontinuous pressure distribution near the trailing edge of the blade. 

Since the Laplace equation 1s linear, a solution to a complex flow field can be built 
by superposition of solutions of elementary flows. The elementary flows used in the 
panel method are the uniform parallel flow and flows about singularities (sources and 
vortices). 

The panel method is based on replacing each blade bv a distribution of sources and 
vortices on its surface. The surface 1s divided into a finite number of straight segments. 
called panels. 

On each panel, a uniform source distribution and a uniform vorticity distribution is 
assumed. The source strength at each panel is set to satisfy the boundary condition at 


the midpoint of the panel (called the control point), and so, in general the source 


strength will varv from panel to panel. The vorticity strength 1s assumed to be the same 
for all the panels and is set to satisfy the Kutta condition. 

The cascade is defined as an infinite row of similar blades, each one modelled Бу 
panels of source and vortex distributions. The flow at each point is found by summing 
the contributions of all the singularities on all the blades , and the uniform parallel flow. 

A useful concept in dealing with such flows is the concept of influence coefficients. 
An influence coefficient 1s defined as the velocitv at a point induced by a unit strength 
singularitv placed at some other field point. The influence coefficients are a function of 
 ШСІТу апа 50 сап бе computed for a given cascade and a given choice of panel ge- 
ometrv. 

Using the influence coefficients, the normal and tangential velocities at each control 
point can be written as a function of the unknown source strength of each panel and the 
unknown Vortex strength. Using the boundary conditions, bv equating the normal ve- 
locity at each control point to the prescribed blowing velocity v,, and using the Kutta 
condition (which requires equal velocities on the upper and on the lower panels at the 
trailing edge), a system of linear equations 1s constructed. 

By solving the system of equations, the strength of the sources and vortices 1s found. 
and the velocities (and the pressures) can be computed everywhere in the flow field. 

The velocity distribution on the surface of the blade, computed by the panel method, 
is used as the boundary condition for the boundary laver flow calculations. 

It should be noted that in the panel method as used in the present computer code, 


there 1s no modelling of the wake, and its effect on the flow field 1s ignored. 


B. VISCOUS FLOW METHOD 
Viscous flow 1s the second building block of the cascade flow and it is applied to the 
thin boundary laver near the blade surface. 
1. Boundary Layer Theory 
The boundary layer concept, first suggested by L. Prandtl, assumes that the flow 
field can be divided into an outer flow where the viscous effects are negligible compared 
to inertia effects, and a thin laver close to the surface where the viscous effects cannot 
be neglected. The complete presentation of the boundary laver theory, and the devel- 
opment of the boundary laver equations, is given bv Schlichting [Ref. 4]. 
Under the assumptions of two dimensional incompressible thin boundarv laver 


flow, the Navier Stokes equations and the continuity equation reduce to: 





Е = 0, 
Ох Oy 
Cu Cu аш, С Cu 
ИРИ. + у DM 
š C ax Cy Oy 
with the boundary conditions: 

y=0 u(x,0) 2 0,v(x,00) 20, 
V= U(X Ve) = u (x) , 


V 
Where b denotes Irma 


Writing the velocity components in terms of a stream function Y : 








This eliminates the continuity equation (which the stream function satisfies by defi- 
nition). 


Introducing the Falkner Skan transformation: 


YR 
n= ue 





the momentum equation and the boundary conditions transform to: 


r^y + BEE pes mt = а). 


z 
P CX 


Ia 0 ЈИ) = Ом Ее») == 0, 


1 = n, M = | 


where m 1s defined bv: 
Ж du, 


Ме dx 


The third order differential equation can be reduced to a svstem of first order 


differential equations by introduction of two new variables U and Y : 





(БҰ) M fV +m(l-— С°) = (t as =v). 


with the boundary conditions: 


0 О =, 


Я 
n = n, Их, я.) = 1. 


The next step 1s to use a finite difference approach to solve the equations. The 
box method is applied using central differencing in both the x and y directions, and 
satisfving the equations midway between nodes. 

Applving the box method results in a system of nonlinear equations in the un- 
known variables (which are f, U and V in each node along the y direction at the current 
X Station). 

In order to solve the nonlinear system the Newton iterative procedure 1s used, 
linearizing the equations first about the solution at the adjacent upstream station, and 


then about the preceding iteration. The linearization is performed by letting: 
,K |, к—1 с, к 
ре ру 
e тк] ayol K 
= уы, SI. 


Z 


ук pl к] ç v ri, K 
аа ы 


м. 
e {denotes location in the x direction 
e j denotes location in the v ( 5 ) direction 


* x indicates the iteration counter 


This linearization results in a svstem of linear equations for the unknown in- 
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This svstem of equations is solved repeatedlv until the changes in the unknowns 
are small enough. Since the system is block tridiagonal, Keller's block elimination 
method is used. 

The method described so far, 1s a direct boundary layer method. It can be used 
as long as the flow does not separate. Whenever separation or flow reversal occurs. and 
a zero skin friction coefficient 1s encountered, the equations become singular and the 
calculations will break down. 

2. Interactive Boundary Laver Method 

The interactive boundary laver method 1s designed to overcome the difficulties 
encountered at regions of flow reversal and separations. In such areas the external ve- 
locity is substantially changed by the viscous effects and can no longer be considered as 
a known boundary condition for the boundary laver flow. 

The general approach to the solution is the same as for the direct method but. 
since the outer flow is unknown, the velocity at the edge of the boundary layer is written 


as, 


Паж 
u(x,y,) = u, + >= rds ка , 


where: 
l. u(x,y,) 1s the total velocity at the edge of the boundary laver. 


2. и(х) is the velocity as computed by the inviscid method. 


y 


d 
i (и) Sy is the Hilbert integral. 


The numerical solution of the boundary laver equations follows the same steps 
as for the direct method, but with some changes. 


The transformation of the stream function and the y coordinate uses a constant 
velocity u, as a scaling factor, and a scaled velocity w is introduced: 





S 
unc vx J? 
Mo. * 
Nu M Их), 
af Uy VX 
Ud xy) 
w = Е 


Using this transformation, the boundary layer equations become a system of 
first order differential equations: 








f' =U, 
U = У, 
W' = 0, 
Tm кте. OU mee? 
NES ae ”-. 


with the boundary conditions: 


ооо, 
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The finite difference box method is used to solve the equations, in the same wav 
as it was used for the direct case, but with two additions: 


l. In areas of flow reversal the term ueu’ox is omited to assure stable mtesranom 
(the FLARE approximation). 


2. The edge velocity, Wj (where J denotes the edge station) which involves integration, 
1s approximated by : 


Wi = + сМ Л), 
where g, and c, are obtained from the numerical approximation to the Hilbert in- 
tegral (which will be presented in the next section ). 

By using central diffrencing to approximate the differential equations, a system 
of nonlinear algebraic equations 1s obtained for the unknown variables (which are 
fi, Ui, Vi and W:). To solve the system of equations, the system is linearized by the 
Newton iterative procedure, and the resulting linear system is solved (for the new un- 
known variables which are the increments óf , ÓU;*, ÓOV** and óW»* ). 

The solution of the system is repeated until the change in the increments 1s 
negligible compared to the preceding iteration, and the whole process is performed again 
at the next downstream station. 

3. Interactive Model 

The interactive model is used to couple the boundary Іауег to the external flow. 
It is needed in areas where strong interaction occurs, and both the boundary layer and 
the outer flow must be solved simultaneously. The interaction model provides the outer 
boundary condition to the boundary layer calculations by adding a correction term to 
the external velocity computed by the inviscid flow method. 

The external velocity is assumed to consist of a potential flow term ( u,,(x) ) and 


a correction termi due tO viscous eHh (s OE 
их) = их) + их). 


The viscous effect is obtained by a surface distribution of sources on the blade (a con- 
cept first suggested by Lighthill [Ref. 5]). The normal velocities at the surface of the 
blade, induced by these sources, displace the streamlines from the surface in the same 


way that the actual boundarv layer displaces them: 


dó (x) = v(x, 5 ) 
Gia и о 


Where у(х, д") 15 the normal velocity at the displaced surface. 


ee Suriace can be approximated to be a flat plate. the normal 


+ 


Assuming thai 
velocity wil be half the local source strength o (x). Assumung also that the inviscid ve- 


locitv does not change across the boundary laver, the local source strength will be: 





ó 
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ИЕ СОТО ОСОЛ O ACA OV the source distribution. 1s the cor- 


rection term to the inviscid velocity, and can be represented by the Hilbert integral: 


Xo id Хь E 
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= SS eS о) 
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The integration 1s carried out on all the sources on the surface. since the horizontal 


Velocity 1s influenced by all the sources. 
The Hilbert integral is then approximated by a finite series: 





f K 
Ca ee Stk 
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Whiere c, 1s a matrix of interaction coefficients which are functions of the geometry onlv 


( ¢ denotes the chordwise position where w,, 1s evaluated and & is the location of the 


Petree which effects 1, ). 
Since the computation of u,, involves values of 6° downstream of the current x 


location, which are not known vet, these terms are taken from the previous iteration 


using a rclaxation formula. 


J. Turbulence Model 
The turbulence model used here is the algebraic eddy viscosity formulation of 


Cebeci and Smith [Ref. 6]. According to the model used in the present computer code, 


the eddv viscosity v, 1s defined by two different expressions, for the inner region and for 


the outer region: 
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The distance from the wall to the point between the two regions, y., 1s chosen such that 


the viscosity will be continuous. 


The intermittency factor, y, is defined by: 


x 
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u dc 
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Where: 
e R 


• (5 is an empirical constant, originallv assigned the value 1200. 


is the Revnolds number based on external velocity and transition location. 


ХІ” 


Cebeci and Bradshaw [Ref. 2, p.246] described a different expression for the 


variable A in the inner region viscosity formula: 





yee 26y | 
(1 — пер”) (, шү | 
СУ / max 
Where: 
ка VU, du, 
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This version of the turbulence model was not implemented in the original computer 
code. During the work on this thesis, the effect of the modified turbulence model was 
investigated. 

A different intermittency distribution was implemented successfully by Rodi and 
Schonung [Ref. 7 ] for transition over separation bubbles. Thev used for G, the ex- 


pression: 


N 100 
Y exp(0.99Tu) ` 


Where Tu 1s the turbulence level in the free flow. This intermittency model was also in- 


vestigated during the work on this thesis. 


5. Transition 
The prediction of transition from laminar to turbulent flow 15 very difficult and 
has to rely on empirical correlations. The relation used here to predict the onset of 
transition 1s a combination of Michel’s method and the e? method, and 1s given bv 
Cebeci and Bradshaw [Ref. 2, p. 153] 


£x ІТ 


p I + 2290. је 


€ xtr 
Where: 


1. В, 15 the Reynolds number based on the momentum thickness at the onset of 
transition. 


2. R, is the Reynolds number based on x at the onset of transition. 


In the computer code, if a lanunar separation 1s detected before transition oc- 


curs, the onset of transition 1s assumed at the point of laminar separation. 


IM. DESCRIPTION OF THE COMPUTER CODE 


The computer code used here to investigate cascade flows was written by Cebeci, 
and is based on the numerical formulation that was outlined in the previous chapters. 
In this chapter the general structure and the major subroutines of the code vill be de- 


scribed. 


A. GENERAL STRUCTURE OF THE MAIN PROGRAM 

The main program reads in the cascade data (blade coordinates, spacing and stagger 
angle), the flow data (inlet angle and Revnolds number), and transition parameters. The 
transition onset on each surface of the blade can be computed bv the program, or can 
be input by the user. The intermittency parameter G should be specified by the user. 

The program then calls subroutine POTNL to compute the outer inviscid flow field 
Jor the first cycle. The output of subroutine POTNL is the external velocity distribution 
on the surface of the blades. This velocity distribution is then transferred to subroutine 
CASBLP, which calculate the boundary layer flow. 

Subroutine CASBLP returns the displacement thickness distribution and the blow- 
ing velocity distribution on the blades to the main program. This data is then transferred 
back to subroutine POTNL to the next cycle of calculations. 

The program repeats the cycles of calculations by calling the two subroutines, until 


the specified number of cycles is reached, or until a convergence criterion is satisfied. 


B. DESCRIPTION OF THE SUBROUTINES 
I. Subroutine POTNL 

This subroutine solves the inviscid outer flow by using the panel method. The 
subroutine calculates the influence coefficients and calculates the velocities subject to the 
boundary conditions. 

The velocities are evaluated on the displaced surface (the surface created bv 
adding the displacement thickness to the original surface of the blade). The input to this 
subroutine includes the cascade geometry, the blowing velocity and the displacement 
thickness (for the first cycle both the displacement thickness and the blowing velocity 


are taken to be zero). 
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2. Subroutine CASBLP 

This subroutine, called by the MAIN program, receives the blade geometry and 
the velocity distribution as input. 

It transforms the x.v blade coordinates to the chordwise tangential coordinates 
and smooths the velocity data (during the work on this thesis it was found that 
smoothing the velocitv data prevents the detection of the separation bubble near the 
leading edge, and therefore it was eliminated). The subroutine then calls subroutine 
COMPBL for further calculations. 

3. Subroutine COMPBL 

This subroutine finds the stagnation point and controls the generation of the 
boundary layer calculation grid for each surface (the grid starts at the stagnation point 
and includes 91 points in the chordwise direction for the upper surface and 71 points on 
the lower surface). 

The subroutine then calls subroutine BL2D which calculates the boundary laver 
parameters for each surface (BL2D is called twice, first for the upper surface and then 
for the lower surface). 

4. Subroutine BL2D 

This subroutine computes the displacement thickness and the blowing velocity 
and returns them back to the calling subroutine (COMBL) in arrays compatible with the 
potential flow calculations (one array that contains all the points of the blade, first the 
lower surface starting at the trailing edge and proceeding forward, and then the upper 
surface, starting at the leading edge and proceeding backwards). 


BL2D calls the following subroutines: 
1. Subroutine INPUT which calculates the following: 


a. NS, the switching point between direct and interactive boundary laver @aleue 
lations (this point is set at the first pressure peak when the blade is scanned from 
leading edge towards the trailing edge) 


b. NTR, transition location (only if the transition location is an input. Otherwise 
it is calculated by subroutine TRNS). 


с. GMTR. the distribution of the intermittency factor y,,. 


In addition this subroutine generates the boundary layer grid in the y direction and 
the initial velocity profile, by calling subroutine INTL. 


2. CALCIJ, calculates the c, coefficients used in the Hilbert integral approximation. 


„ә 


EDDY, calculates the eddy viscosity (called only after transition has been de- 
tected): 
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с 5 the coeliicients Of the boundary laver finite difference 
equations in transformed form ( for the direct method calculations). 


CA 


SOLVE3, solves the linearized boundary laver equations for the F,U and V vari- 
ables by computing the increments óF, óU and ôV 

The subroutine then checks the convergence of the Newton iterations and re- 
peats the calculations if needed. If the subroutine detects flow separation or if it reaches 
the switching point NS, subroutine MAIN2 is called for the interactive method calcu- 
lations. Otherwise, the subroutine proceeds to the next chordwise point of the grid (NN) 
and repeats the calculations. 
5. Subroutine MAIN2 

This subroutine calculates the boundary laver parameters by the interactive 


method. The subroutine performs the following steps: 


]. It first calls subroutines JOINT and COMGI to compute the interaction coeffi- 
clents. 


2. In regions of lamunar flow it calls the following subroutines: 


a. COEF. which calculates the coefficients of the boundarv laver finite differences 
equations. 


b. SOLV4, solves for the variables F, U, V and W by computing the increments 
OF, oL ,oVandoW. 


c. TRANS, to check if the condition for transition is satisfied (it also checks for 
laminar separation and initiates transition at the point of lanunar separation if 
It is detected). 


The subroutine then checks for convergence of Newton iterations and repeats the 
calculations as needed. 


3, In regions of turbulent flows the subroutine calls the following subroutines: 
a. EDDY, to compute the eddv viscosity parameter B. (B — 1 4 v,/v ). 
b. COEF and SOL V4, the same as for laminar flow. 


6. Subroutine OUTPUT 
This subroutine computes the boundary layer parameters. It is called with a 


parameter "INDEX" which determines the tvpe of calculations: 


l. For INDEX=1 the computations relates to transformed coordinates (direct 
boundarv layer method) using the relations: 


В) 
si Rope 
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Where V(1,2) and B(1,2) are the velocity gradient and the viscosity parameter at the 
surface, respectively, and 7(NP) and F(NP) are y and F evaluated at the edge of the 
boundary laver. 


2. For INDEX= 2 the subroutine calculates the boundary layer parameters lor semi- 
-transformed coordinates (interactive boundary laver method) using the relations: 


2 V(1.2) В(1.2) 
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go (n 5 Nx 


For NX > NTR (after transition has been detected), subroutine SMPSON 15 
called (subroutine SMPSON calculates the coefficient of the outer region eddy viscositv). 
The subroutine then prints out the velocity profiles at the required stations. 

7. Subroutine TRANS 

This subroutine calculates the transition location based on the Michel criterion 
or based on laminar separation (whichever occurs first). If transition has been detected 
the intermittency distribution is calculated for all the remaining points of the surface. 

8. Subroutine FILLUP 

This subroutine increases the number of points in the boundary layer grid (in 
the y direction) as needed. It also fills up the arrays of F, U, B, W and V between the 
edge of the boundary layer to the end of the arrays (with V=0, W,B and U, with the last 


values that they had in the edge of the boundary layer and F as the integral of U). 


9. Subroutine EDDY 
This subroutine calculates the eddy viscosity using the Cebeci--Smith two laver 
eddy viscosity formula. It receives the vectors L,V and s at a point and computes the 
viscosity vector B. 
10. Subroutine INTL 
This subroutine generates the boundary layer grid in the y direction. It sets the 


number of grid points and generates the initial velocity profile. 
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IV. RESULTS AND DISCUSSION 


The viscous inviscid interaction code was run with several cascades on which ех- 
perimental data 1s available. In order to enable a thorough comparison between exper- 
imental results and the computed results, a very detailed experimental data base is 
needed. The data should include measurements of the boundary laver development along 
the blade, velocitv profiles along the boundary laver, transition location and distribution, 
flow separation, and external velocity distribution. 

Unfortunatelv, verv few cascade experiments has been performed, which obtained 
the required data with sufficient accuracy, due mostly to the lack of appropriate meas- 
urement equipment. Only recentlv, with the introduction of non--interfering methods 
like the Laser Doppler Velocimeter (LDV), the required data can be measured accu- 
rately. 

Recently an experiment involving the investigation of a linear compressor cascade 
of Controlled Diffusion Blading (which will be referred here as the CD cascade) has been 
carried out by Elazar [Ref. 8]. Most of the work in the present thesis, involves compar- 
ison of the computer code results with Elazar's experimental results. 


Other cascades that were investigated are: 


1. A shockless, supercritical airfoil cascade, designed in 1974 bv Korn in cooperation 
with Pratt & Whitney Aircraft (referred here as the P & W cascade). The exper- 
imental results of the cascade were obtained from a report by Hobbs, Wagner, 
Dannenhofler and Dring [Ref. 9]. 


2. Stator blade of a single stage axial compressor (referred here as the C4 cascade). 
The blade profile is the British C4 section (10% thickness) on a circular arc camber 
line. The experiment has been performed by Walker [Ref. 10]. The detailed 
boundary layer measurements are not presented in the report and were obtained 
directlv from the author. 


The code failed to run with two other cascades: 


1. A highly loaded, double circular arc blade with a sharp leading edge and a sharp 


trailing edge, used in a compressor cascade that was investigated by Deutsch and 
лет Шекер 


2. V2 double circular arc blade, highly loaded cascade. This cascade was investigated 
by Hoheisel and Sevb [Ref. 12]. 


In both cases the code calculated the potential flow successfully but failed in trying 


to compute the first cycle of the boundary layer calculations. 
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ыы ОСА Е 

ее Ча Гог the CD cascade was obtained at М = 0.2%, В, = 700000 
and at three inlet angles: 40? (the design condition). 43.437 апа 46° . The spacing was 
Шо tie chord and the Stageer angle 14.27" . A general lavout of the cascade 15 shown 
im Figure ] on page 20. 


The following observations were concluded from the experiment: 


1. А separation bubble exists near the leading edge on the upper surface at all the inlet 
angles. The bubble became larger at increased inlet angles. 


2. Transition from laminar to turbulent flow occurred above the separation bubble 
(on the upper surface). 
3. Transition on the lower surface occurred at midchord. 


4. The boundary laver thickness on the upper surface increased with inlet angle. and 
КОО a thickness Of 13°>o chord at the highest inlet angle. The boundary laver 
thickness on the lower surface did not change significantly with inlet angle. 


5. The turbulent boundary laver on both surfaces remained fully attached at all the 

Inlet angles. 

l. Transition location and intermittency distribution 

The effects of the transition location and the intermittency factor were investi- 

gated. The code was first run with the transition location calculated by the code. and 
with several values of the internuttency factor G,. It was found that the code did not 
run with G, = 1200 (which is the value used usually for high Reynolds numbers). The 
highest value of G, with which the code run successfully was 900. 

The code failed to predict the separation bubble on the upper surface. and pre- 
dicted lanunar separation at 78% chord on the lower surface (which did not occur in the 
experiment). Transition on the upper surface occurred at 41°%> chord (detected bv 
Michel's criterion) and at 78% chord on the lower surface (at laminar separation). 

The shape factor computed by the code was compared to the experimental re- 
sults. As can be seen in Figure 2 on page 21 the shape factor as predicted by the code 
deviates substantially from the actual results, due mainly to the different transition lo- 
cation. 

On the lower surface, as can be seen in Figure 3 on page 22 the shape factor 
deviates even more from the experimental results. In this figure the effect of changing 
the intermittency factor G, can be seen. For both the extreme values of G, . 10 and 900, 


the computed shape factor curve is far from agreement with the actual results. 
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CD CASCADE 


Figure 1. Controlled Diffusion cascade 
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Figure 2. Shape factor comparison on the upper surface: Transition computed by 
e code (L = 40°) 
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Figure 3. Shape factor comparison on the lower surface: Transition computed by 


the code (f + 40”). 


Е вое ОККО рейс the existence of the scparation bubble on the 
upper surface, it was decided to try to eliminate the smoothing oí the external velocity 
as computed by the potential flow subroutine (originally, the velocities were smoothed 
in subroutine CASBLP prior to boundarv laver calculations). It was found that without 
smoothing the velocities a small separation bubble 1s predicted by the code at 4%o of 
chord. The onset of transition is set by the code at the beginning of the separation 
bubble. 

The code was run with unsmoothed velocities with two values of G,, 10 and 900. 
The shape factor behavior can be seen in Figure 4 on page 24. Changing the value of 
G, did not change the shape of the curve much, and generally the shapes of the com- 
puted and the experimental curves look alike. 

The elimination of the velocity smoothing in the code, also affects the thickness 
of the boundary laver. In Figure 5 on page 25 the displacement thickness is plotted for 
both cases (with and without the velocity smoothed). Without smoothing, the displace- 
ment thickness 1s much thicker, especially on the rear half of the blade, which is closer 
to the actual results. 

The effect of changing the intermittency distribution to the one used by Rodi 
and Schonung [Ref. 7] was investigated. It was found, as can be seen in Figure 6 on page 
26 that the effect of the new model is equivalent to using G, in the present model. 

On the lower surface it was necessary to run the code with transition as input, 
to get reasonable results, as can be seen in Figure 7 on page 27 for transition input at 


Ino of chord. 


At the off design conditions (inlet angles of 43.4? and 46^) a similar behavior of 
the transition has been observed. as can be seen for example in Figure $ on page 28 for 


ШЕШЕСІ асе амжил Figure 9 on page 29 for the lower surface, both at fp = 46°. 
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Shape factor on the upper surface without velocity smoothing. 
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Figure 5. Displacement thickness: The effect of velocity smoothing. 
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Figure 6. The effect of the intermittency model: Upper surface, f = 40° 




















O 

E 

ио 

p. 

LEGEND 

| o EXPERIMENTAL — 
E^ | TIRANSITION INPUT 
Eo АА 
о |! 
ar 
El - 
A, 
< 
am | 
Ае —— s 

ES 

Ы? 

E 

O 

ox 

0.0 Ше 0.4 0.6 0.8 1.0 


X/C 


Figure 7. Shape factor on the lower surface with transition input at 21%. 
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Figure 8. Shape factor at Bf = 46° on the upper surface 
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Figure 9. Shape factor at f = 46? on the lower surface. 
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2. External Velocity Distribution 
The external velocity distribution, computed by the code using the interaction 
law, was compared to the experimentally measured velocities. It was found that in gen- 
eral, the velocities measured experimentally, were higher than those computed by the 
code for all the inlet angles. 
There are two possible sources to the discrepancy in the velocities: 


1. The computer code calculates pure 2--D flows. In the experiment the flow was 
observed to accelerate due to the effect of the boundary layer on the side walls (a 
3--D effect). This effect was calculated in the experiment and is referred to as the 
AVDR correction [Ref. 8, p.33]. 


The flow accelerates due to the thickening of the boundary laver. Since the 
boundary layer as computed by the code is substantially thinner than the actual 
boundarv laver (as will be discussed in the next section) the external velocities 
predicted by the code are smaller. 


ә 


To compensate for the first error source, all the computed velocities were com- 
pared to the experimental velocities corrected by the AVDR correction. The comparison 
between the velocities can be seen in Figure 10 on page 31 for + 40^, in Figure 11 
on page 32 for B = 43.4° and in Figure 12 on page 33 for f = 46°. 

It can be seen from the figures that the difference between the computed and the 
experimental velocities is larger on the lower surface. The reason might be the method 
with which the correction to the inviscid velocity 1s computed. The assumption on which 
the interaction law 1s based, is that only sources (representing the viscous effects) on the 
surface being considered, affect the local velocity. In reality, the boundary layer on both 
surfaces affects the local velocity (because the boundary laver developed on the upper 
surface of a blade, causes a velocity disturbance that is felt on the lower surface of the 
adjacent blade). 

Since the boundary layer on the lower surface is much thinner, its effect on the 
velocity on the upper surface is much smaller than the effect of the upper surface 


boundary layer on the lower surface velocity. 
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Figure 10. External velocity at p = 40° 
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Figure 12. External velocity at f = 46° 
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3. Boundary Layer Thickness 

The boundary layer thickness as computed by the code was coinpared with the 
experimental results by comparing the displacement thicknesses. 

It was found that on the lower surface the computed and the actual displace- 
ment thickness agree quite well, as can be seen in Figure 13 on page 35 for — и 
Figure 14 on page 36 for f = 43.4” and in Figure 15 on page 37 for f = 46°. 

On the upper surface the displacement thicknesses computed by the program 
are significantly thinner than those measured experimentally. The difference between the 
computed and the actual thickness increases along the blade and it increases with 1n- 
creased inlet angle. It was found that by using a different expression for the inner region 
eddy viscosity (as mentioned in chapter II), the displacement thickness can be increased, 
but the difference between the actual and the computed thickness 1s still substantial, es- 
pecially at the higher inlet angles. Figure 16 on page 38, Figure 17 on page 39 and 
Figure 18 on page 40 shows the displacement thickness on the upper surface for the 
three inlet angles, with the original and the modified eddy viscosity models. 

The large error in the prediction of the boundary laver thickness, can be the re- 


sult of several reasons: 


1. The transition model used in the code, sets the onset of transition at the first point 
of laminar separation. It causes rapid transition to turbulent flow which reattaches 
immediately, resulting in a very small separation bubble compared to the bubble 
observed in the experiment. 


2. The turbulent model used in the code could be inaccurate. It was derived based on 
empirical data obtained in single airfoil experiments and not with cascades. In ad- 
dition the present model does not include the effects of the free stream turbulence 
(that was relatively high in the experiment). 


3. The boundary laver as measured in the experiment was quite thick. especially at the 
higher inlet angles (it reached 15% of the chord at f — 46^). Such a thick bound- 
ary laver may violate the basic assuinptions on which the boundary laver 
equations, and the interaction law, were based (especially when the spacing be- 
tween the blades is small, 60% chord in this case). 


It was suggested that one of the possible reasons to the inaccurate prediction 
of the boundary laver 1s the blunt trailing edge of the blade, that might cause difficulties 
in the computations. A modified blade, with a sharp trailing edge has been run, and the 
displacement thickness distribution can be seen in Figure 19 on page 41. As can be seen 
in the figure the sharp trailing edge affects only the boundary layer adjacent to the 


trailing edge, and therefore cannot provide an explanation to the difference between the 


actual and the computed displacement thickness. 
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Displacement thickness on the lower surface (P = 40? ) 








LEGEND 














o EXPERIMENTAL 
c .... COMPUIED . . 


ве ШЕЕ 
0.010 0.015 


0.005 


0.000 





Figure 14. Displacement thickness on the lower surface (f) — 43.4? ) 
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Figure 15. Displacement thickness on the lower surface (f = 46° ) 
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Figure 16. Displacement thickness on the upper surface (f) — 40? ) 


25 


0.010 0.015 0.020 0.023 OORS O 
! 


DIST DACE MONT HIG KNEES 


0.005 


0.000 


LEGEND 
оты ГА Т, Е 


ORIGINAL 





И 


n / 
E 
Z 
/ 


Y 
Я 22 
” 
n 4 


Al 
N 
` 
` 
` 
` 
` 


д ос Оз м Шо 06 0.7 0.8 


Figure 17. 


Displacement thickness on the upper surface (fj — 43.49) 
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Figure 18. Displacement thickness on the upper surface (B = 46? ) 
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The effect of sharp trailing edge. 
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4. Comparison to a Navier Stokes Code 

A hnuted comparison of the experimental and the computed results with a 
Navier Stokes (X.S.) code calculations has been performed TheM S coge ha E аан 
veloped and run bx S. J. Shamroth of Scientific Research Associates inc. 1n cooperation 
with Pratt and Whitnev Aircraft. 

Since the N.S. code does not compute the displacement thickness. the velocity 
profiles near the surface of the blade were compared. The comparisons were made at 
90° chord on the suction surface for all three inlct angles. 

At the design pomt. 6 = 40%, shown in Figure 20 on page 43. both the inter- 
active code and the N.S. code failed to predict accuratelv the actual velocity рош 
this case the interactive code seems to vield somewhat better results than the NC Se 

At the higher inlet angles, 6 = 43.4° and f = 46^, shown in Figure 21 on page 
44 and in Figure 22 on page 45 respectively, the N.S. calculations show significantly 
better agreement with the experimental results than the interactive code. 

From these comparisons, it can be seen that the interactive code deviation from 
the actual results increases with increased inlet angle (increased loading of the cascade), 


whereas the N.S. code deviation Seems teadecrease with mereased inlet ancies 
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Figure 20. The results of the N. S. code at B = 40? 
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Figure 21. The results of the N. S. code at f = 43.4? 


44 


1.0 

























C 
£O 
ще 
E 
то 
~ 
| LEGEND 
. o EXPERIMENTAL 
У INTERACTIVE CODE 
E | 
© +f | | 
0.00 0.02 0.04 0.06 0.08 0.10 162 
D/C 
BETA=46 


Figure 22. The results of the N. S. code at f| = 46? 


B. P& W CASCADE 

The experimental data for the P & W cascade was obtained at inlet flow angle of 
52°, at M=0.11, and Revnolds number of 478000. The cascade had a stagger angle of 
15.75° and 0.7 spacing. A general layout of the cascade is shown in Figure 23 on page 
47. 

A comparison of the computed and the measured pressure coefficients on the blade 
is Shown in Figure 24 on page 48. There is a good agreement between the computed 
and the measured C,. 

The displacement thickness was measured in the experiment only at 96.8% of chord. 
This measurement is compared to the computed results in Figure 25 on page 49. As can 
be seen, the computed and the measured data agree almost perfectly on the lower sur- 
face, and quite well on the upper surface. The difference observed on the upper surface 
is caused bv the earlv prediction, by the code, of trailing edge separation, a short dis- 
tance upstream of the actual location. This can also be observed when comparing the 
velocity profiles at that point, in Figure 26 on page 50. The computed velocity curve 
shows a small zone of reversed flow near the surface of the blade. This reversed flow 
could be the result of a too early prediction of trailing edge separation by the code, or 


it could have existed in the actual flow but not detected because of its size. 
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P&W CASCADIs 


Figure 23. Pratt & Whiteny cascade 
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Figure 24. Comparison of pressure coefficient. 
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96.8% chord. 
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Figure 26. Velocity profile at 96.8% chord on the upper surface. 
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Cie Gas CAD TL 


nw cade Was a Siaecer ancle of 2975°7a camber angle of 31.1 


° 


and spacing 
of 0.992, 1t has been tested at Revnolds numbers of about 200000, and inlet angles of 
meee) (which comesponds to memdence angles of -10.9° to 2.7°). The general 
lavout of the cascade 1s shown in Figure 27 on page 52. A computer code that generates 
the coordinates of the blade and a summary of some experimental results arc given in 
Appendix B. 

The code was run with the intermittency constant G, = 10. Higher values of G, 
(above 100) caused numerical problems in the code. The onset of transition was first 
taken at the point where 1t was observed in the experiment. At the lower inlet angle it 
seems that a better agreement with the experimental results can be obtained bv delaying 
the onset of transition but trying to implement it resulted in numerical breakdown of the 
computation. At the higher inlet angles. better agreement with the experimental results 
was achieved by initiating the transition earlier (at 2676 chord for f — 45.6” and at 21% 
for f =47.7° as compared to 44% апа 56% chord as observed in the experiment). 

1. Displacement Thickness 

Comparisons of the experimental data to the computed displacement thickness 
are shown in Figure 28 on page 53 for inlet angle of 34.1°. in Figure 29 on page 54 for 
inlet angle of 36.3°, in Figure 30 on page 55 for inlet angle of 45.6° and in Figure 31 
on page 56 for inlet angle ol 47.7°. 

‘As can be seen in the figures, there 1s a good agreement between the actual and 
the computed results at the two lower angles (6 = 34.1° and f = 36.3°, in which the 
incidence angles were negative). At the two higher angles. 6 = 43.6° and 6 = 47.7° the 
computed results agree with the actual results up to about 70% chord, and then the 
displacement thickness predicted by the code becomes much thicker than the actual one. 

The code predicted a large flow separation area starting at about 70°% chord at 
the lower inlet angles, and at about 46% chord at the higher inlet angles. This flow 
separation was not obsc-ved in the experiment. The discrepancies between the computed 
and the actual results behind 60% то 70% chord can be explained bv the inaccurate 
calculations bv the code due to the large separated areas. When the code encounters 
separation, several approximations are made (like the FLARE approximation) based on 
the assumption that the separated area 1s small. When the separated area 1s large, these 


approximations may result in inaccurate prediction of the flow field. 


CT CASUM 


Figure 27. C4 Cascade 
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Figure 28. C4 cascade at B = 34.1%: Displacement thickness comparison with 


computed results. 
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Figure 29. C4 cascade at В = 36.3°: Displacement thickness comparison with 


computed results. 
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Figure 30. C4 cascade at B = 45.6°: Displacement thickness comparison with 


computed results (G = 10). 
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Figure 31. C4 cascade at B = 47.7°: Displacement thickness comparison with 


computed results (G = 10). 
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2. External Velocity and. Velocity Profiles Comparisons 

A comparison of the external velocity on the upper surface of the biade 16 shown 
nuuc" on расе SS for inlet angle of 45.6° and in Figure 33 on page 59 for inlet 
о! tcan be seen tliat there is a good agreement between the experimental 
and the computed results up to about $0“ chord. Near the trailing edge the computed 
resuits deviate from the experimental results due to the inaccuracy in the calculations 
of the displacement thickness. 

орой об Сао рге йе за те boundary laver at 50°% chord is 
shown in Figure 34 on page 60 for inlet angle of 34.1° and in Figure 35 on page 61 for 
ШИЕ Of 30.3. The agreement between the calculated velocity profiles and the 


measured velocity profiles is very good. 
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Figure 32. C4 cascade at B = 45.0°: 
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Figure 33. C-+easeade at f = 47.7%. External velocity distribution. 
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Figure 34. СУ сазсаде а д = 34.19: Velocity profile at 30% chord. 
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Figure 35. C4 cascade at B = 36.3° Velocity profile at 50% chord. 
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У. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The interactive viscous inviscid computer code, has been investigated by comparing 


its predictions of boundary laver parameters to experimental data. 


It has been found that the code vields reasonable results for lightly loaded cascades, 


but the prediction of the boundary laver thickness on the suction surface of highly 


loaded cascades deviates significantly from experimentally measured data. In two cases 


involving highly cambered cascade blades. with sharp leading edge, the code failed to 


run. lt was also found that the prediction of external velocity distribution on highly 


loaded a лыс уле лас л 


The maim reasons to the discrepancies in the prediction of the boundary a 


thickness seen) to be: 


ЈЕ 


5 


Јпассигасу in predicting flow parameters in regions of large flow separation (due 
to inadequate transition model and approximations made m calculating the flow in 
separated areas). 


Inaccurate turbulence modelling. 


Possible violation of the basic assumptions of the boundary laver theorv in areas 
of very thick boundary lagen 


The wake is not calculated by the code. The result is inaccurate flow prediction 
near the trailing edge. 


The maccuracy in the prediction of the external velocity distribution in highlv loaded 


cascades is due to the interaction law, which does not account for the presence of adja- 


cent blades. 


B. RECOMMENDATIONS 


Lherécommiended steps in order to Iñini ose the codeine: 


Improving the interaction law by assuming a distribution of sources on the actual 
surface (instead of the assumption of a flat plate). letting the correction term to the 
external velocity vary across the boundary laver and distributing sources on the 
adjacent blade as well for better modelling of the boundary laver effect on the ex- 
ternal velocity. 


tv 


t 42 


Changes to the derivation of the boundary layer equations should be investigated 
to allow a better treatment of thick ooundary lavers (like onutting the assumption 
(И —slnacioss the boundary layer). 


Different turbulence models should be investigated. 


The wake should be included in the calculations. 
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APPENDIX A. COMPUTER CODE LISTING 


le elo afa eta eta elo als nto nts s'e sz alarma ts ‘. 
7% ғ» 7% £9 2. 


кисти VISCOUS- INVISCID INTERACTION PROGRA p CASCADE FLOWS 


беја пе оће ventas 
2% sie er 


THIS VISCOUS-INVISCID INTERACTION METHOD, CAPABLE OF COMPUTING BOTH 
WAS DEVELOPED BY CEBECI AND 
COLLABORATEURS AT LONG BEACH STATE AND DOUGLAS AIRCRAFT COMPANY. 
THE CODE APPLIES TO СОВЕТОВ 
PAST LINEAR, ARBITRARILY STAGGERED CASCADES. THE METHODS BASIC 


5 


I 


ds 
2. 


5 


4. 


V E R ТО 
JANUARY 87 


3. А 


INGLE AIRFOIL AND CASCADE FLOWS, 


NGREDIENTS INC .UDE 
À FIRST ORD;3 PANEL METHOD TO SOLVE LAPLACE'S EQUATION, 


SUBJECT TO DIRECT OR INTERACTIVE BOUNDARY CONDITIONS, 
RESULTS, AND 


A ZERO EQUATION, ALGEBRAIC TURBULENCE MODEL TO ESTIMATE 
TURBULENT SHEAR STRESSES. 


2-DIMENSIONAL, STEADY FLOWS 


A STRONCE “A ILRACTION MODEL TO COUPLE VISCOUS AND INVISCID FLOW 


кик 1М7Т0001 
Пи EE 
еки INTOOO3| 


INTOOO4 

INTO005 
ІКТ0006 
INTO007 
INTOO08 
ІКТ0009 
INTOO10 
INTOO11 
INTO012 
INTOO13 
ІМТ0014 


A FINITE DIFFERENCE SCHEME TO SOLVE THE BOUNDARY LAYER EQUATIONSINTOO 


INTOO16 
ІМТ0017 
INTOO18 
INTOO19 
INTOO2Q0 
INTOO21 


IN SUMMARY, THE CODE WILL PROVIDE, FOR ATTACHED AS WELL AS MODERATE-INTOO22 


LY SEPARATED FLOWS PAST SINGLE AIRFOILS OR CASCADES, THE FOLLOWING 
J. 
2 


E 


c 


Whe 


INVISCID AND VISCOUS PRESSURE DISTRIBUTIONS, 
DISTRIBUTIONS OF 

A. LOCAL SKIN FRICTION COEFFICTENT, 

B. DISPLACEMENT AND MOMENTUM THICKNESS, AND 
VELOCITY PROFILES ACROSS THE BOUNDARY LAYER. 


IODIFICATIONS SINCE VERSION 3. 0: 


PRECISE ASSIGNMENT OF BEGIN OF TRANSITION. 


ADDITIONAL PRINT OPTION: 
NUMBER 12) FOR THE PLOTTING ROUTINE. 


COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 

COMMON / BLOW/VN( 100) 

COMMON/BLIN/ TITLE(20) ,XC(100),YC( 100), ISG(100) ,DELS( 100), 
n XCTR,XTR,ISTRP,ICYCLE,ICYTL,XCTRS(2) , TRFIND(2) 

COMMON/CASCDE/INLET,SP,SINGLE,ALPHAA , ALPHAI , STAG 

COMMON/TRN/ PGAMTR,OMEGA ,RTHETB ,RTRANB 

COMMON/PLOT/NVP(2) ,NXVP( 20,2) , ICC 

DIMENSION X0(100),Y0(100),X(100) , Y(100), VCOM( 100) ,DLS( 100) , 
- XS(100),Y8S(100) , XSTGR( 100) , YSTGR( 100) ,DBPP( 100) 

DIMENSION CASEID(20) ,XCTRI(2) , ITRIC(2) ,NBL(2) 

LOGICAL SINGLE,TRFIND 

TRFIND(1)= . FALSE. 

TRFIND(2)= . FALSE. 
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INT0023 
INT0024 
INT0025 
INT0026 
INT0027 
INT0028 
INT0029 
INTOO3C 
ІХТ0031 


CORRECTION OF AN ERROR IN THE CALCULATION OF MOMENTUM THICKNESS. 1№Т0032 
ТР=-2 WILL PROVIDE AN INPUT FILE (UNIT INTOO3 


ІМТ0034 
INTOO35 
INTOO36 
INTOO3 

INTOO38 
ІМТ0039% 
INTOO40 
INTOO4 

INT0042 
INTO043 
INTO044 
ІХТ0045 
INTOO46 
INTO047 
INTOO48 
INTOO49 


20 


5 


7 


80 
82 


IGASE = 0 
READ(5,5,END=999) TITLE 
FORMAT( 2044) 
ICASE = ICASE + 1 
REWIND 3 
READ (5,10) 
FORMAT( 1X) 
E08 E70). ETRI(CT),ITRIC2) ,IRST, IGCYTL,IP 
FORMAT( 1615) 
READ (5,10) 
READ (5,25) INLET,ISTAG,ALPHAI,STAG, SP, PGAMTR , OMEGA 
FORMAT(215,5F10. 0) 
READ (5,27)RN,XCTRI(1) ,XCTRI(2) , ALPHAA 
FORMAT(4E10. 0) 
IF (IP.EQ.-2) THEN 
READ (5,20) NVP(1),NVP(2) 
IF (NVP(1).NE.0) READ (5,20) (NXVP(I,1),I=1,NVP(1)) 
IF (NVP(2).NE.0) READ (5,20) (NXVP(1,2),1=1,NVP(2)) 
END IF 
IF (ICASE .EQ. 1) READ (5,20) N,NI 
IREAD = 1 
TBLOW = 1 
SINGLE = . FALSE. 
IF (SP .LE.- 0.0) SINGLE = . TRUE. 
Me N - 1 
№1= N + 1 
ШІ (ІСАСЕ .GT. 1) THEN 
Ni = N1SAVE 
BEN - 1 
GOTO 53 
END IF 
IF (IREAD .EQ. 1) GO TO 40 
READ (5,10) 
READ (5,305 (XO(I), YOCI), I=1,N+1) 
FORMAT( 2Fi0. 0) 
G9 TO 50 


READ(5, 10) 
READ(5,45) (XO(I) , I=1,N+1) 
READ(5, 10) 
READ(5,45) (YO(I) , I=1,N+1) 
FORMAT( 6F10. 0) 


CONTINUE 
ip (IP. EQ. -2) THEN 
WRITE( 12,20) N+1,NVP(1),NVP(2),90,70, INLET 
IF (NVP(1).NE.0) WRITE(12,20) (NXVP(I,1),I=1,NVP(1)) 
IF (NVP(2).NE. 0) WRITE(12,20) (NXVP(I,2),I=1,NVP(2)) 
IF (INLET. NE. 1) WRITE(12,80) RN,ALPHAA 
IF (INLET. EQ. 1) WRITE(12,80) RN,ALPHAI 
WRITE(12,82) (XO(I),Iz1,N*1) 
WRITE(12,82) (YO(1),1=1,N+1) 
FORMAT(2E15. 5) 
FORMAT(8F10. 6) 
END IF 


INTOO500 
J í 0510 
1100520 
Шә ЗО 
INTOO540 
IS T00550 
INTODS6O 
INTOO570 
ІХТ00580 
и 0590 
ІМТ00600 
INTOO610 
ІКТ00620 
INTO0630 
INTO0640 
ІХТ00650 
ІМТ00660 
ІМТ00670 
INTOO680 
INTOO690 
ІХТ00700 
INTOO7 10 
А О ZO 
ІМТ00730 
ІХТ00740 
INTOO750 
INTO0760 
TN TOO 7220 
INT00780 
INT00790 
INT00800 
INT00810 
INTO0820 
INTOO830 
INTOOS840 
INTO0850 
ІХТ00860 
INTO0870 
INT00880 
INT00890 
INT00900 
ІХТ00910 
INT00920 
INTOO9 30 
1NTOO940 
ІХТОО950 
INTO0960 
INTOO970 
INTOO980 
1NTOO990 
INTO1000 
INTO1010 
INTO1020 
INTO1030 
INTO1040 
INTO1050 


Dr 


B 


C 


NRITE = (Ni 

IMIN = (NDI 

IF((N1/272) -EQO М) ПИ = R 

CALL TRGRID (N1, OO NIE 0.5, ИИ, ВАВТ) 
NISAVE = Nl 


CONTINUE 

ALPHAA = 0.017&533 * ALPHAA 
ALPHAI = 0.0174533 * ALPHAI 
STAG = 0.0174533 * STAG 


IF (INLET. EQ ОЕ 
ALPHA = ALPHAA 

ELSE 
ALPHA = 
END IF 


ALPHAI 


IF (ISTAG . М. 0D) THEN 
CALL STAGR(N1,STAG,XO,YO,XSTGR,YSTGR) 


ELSE 

DO 55 I = 1 , N1 
XSTGR( I) = XO(I) 
YSTGR(I) - YO(I) 
CONTINUE 

END IF 


C READ DATA FROM VISCOUS CAL. 


C 


60 
C 


70 


299 


+ 


+ 


ICYCLE 
ICYCLE 


0 
TC Ye i l) 


CALL POTNL(N1,IRST,ALPHA,CHORD,XO, YO, XSTGR, YSTGR , X, Y, DLS , VCOM, 
DBPP) 

ТЕ (ТОО acm 

REWIND 3 

WRITE (3)N1,(XOCI),YO(I),DLS(1I),VN(I) ,DBPP(I),I=1,N1) 

Ooms 

END IF 


ICYTL) THEN 


ІР ОБЛАС МЕ ОИЕ 


ро 70 І = 1, М1-1 

Х(І) = и “© (XO(1)+X0(I+1)) 
Ү(І) = 0.5 * (ҮО(І)+ҮО(І+1)) 
я аи 

END IF 


CALL CASBLP(N1,X0,YO,X,Y,XS, YS, DLS, VCOM, DBPP, RN 
,NBL,ITRI,XCTRI,TITLE) 

GO TO 60 

CONTINUE 

STOP 

END 


SUBROUTINE POTNL(N1,IRST,ALPHA,CHORD,XO, YO, XSTGR,YSTGR,X,Y,DLS$S, 


VOOM, DEPP) 


COMMON/BLOW/VN( 100) 
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| 
ІМТ01060 
INTO1070 
INTO1080 
INTO1090 
INTO1109 
INTO1110 


INTO112 
ІМІ011 


INTO114 
INTO1150 


INTO116 


INTO1170 
INTO1180 


ІМТ01194 
ІМТ0120 
INTO12@ 
ІМТ0122 


INTO1230 
INTO1240 
INTO1258 


INTO126 


INTO129Y 


INTO1300 
ІКТ01310 
INTO1320 
INTO1330 
INTO1340 
INTO1350 
INTO1366 
INTO1370 
INTO138 
INTO1390 
INTO1400 
INTO1410 
INTO1420 


ІМТ0143 


INTO144 
INTO1450 
INTO1460 
INTO1470 


INTO148! 


INTO1490Q 
INTO1500 


INTO15H 


INTO1526 
INTO1530 
INTO1540 
INTO1550 


INTO156( 


INTO1570 
INTO1580 


INTO15 9§ 


INTO1600 


INTO1610 


INTO1270 
INTO1280 





EU PLEO, NAQNST NP NPT ,NTIR,IT,INVRS,NS,IP 
ШЕ БЕТМ/ТІТБЕ(20),ХС( 100) ,ҮС( 1007, 186(100),ВЕІ,5( 100), 


n XCTR,XTR,ISTRP,ICYCLE,ICYTL,XCTRS(2) , TRFIND(2) 


COMMON/CASCDE; INLET,SP,SINGLE,ALPHAA, ALPHAI , STAG 


C SIMPLE SOURCE POTENTIAL CODE 


10 


л 


ÉS 


16 


20 


+ + + + + + + + 


DIMENSION AOFF(100,100), BOFF(100,100), XP(100), YP(100),X(100), 
$(100),C( 100), D(100),VTAN(3,100) , VNOR( 3,100) ,R( 3,100), 
VCOM( 100) ,SIGCOM( 100) , CP( 100) ,X0( 100) , YO( 100) 

,VNC(100) ,D1(100),D2(100),D3(100),8S0(100),SC( 100) 
,XOFF( 100) , YOFF( 100) , T( 3,100) , VTCOM( 100) , VNCOM( 100) 
,XS$(100) , YS( 100) , SOFF( 100) , COFF( 100) , XPOFF( 100) , YPOFF( 100) 
,Y(100) ,816( 3,100) ,DLS( 100) ,DLSC( 100) ,4(100, 100) ,B( 100,100) 
, XSTGR( 100) , YSTGR( 100) , VUT( 3) , VLT( 3) , VUN( 3) , VLN( 3) , DBPP( 100) 
,CPI(100),X08(100),YOS(100) , DBPPC( 100) 

REAL NUM1 , NUM2 

LOGICAL OFF ,SINGLE 


ШЕЕ = . FALSE. 
HI - 3. 141592 
CM = 0.0 

= N] - 1 


Ше (ТСУСЪЕ .EQ. 1) THEN 
ШЕ (ІКЕТ .EQ. 0) THEN 
BO 10 I=1,N1 


 ШІ5(1І) = 0,0 

VN (I) = 0.0 
DBPP(I)= 0.0 
CONTINUE 

ELSE 

10 51= 1, м 
OSI) 2 XO(I) 
MSSCr) = YO(I) 
CONTINUE 


READ (3) NT,(XSCI), YSCI),DLSCC I), VNCC D) ,DBPPC( I) , I21,NT) 
XMIN = XS(1) 

Poms T — 2 , NH 

ШЕ (Х5(1) .GT.WNMIN) GOTO. 15 
XMIN = XS(I) 

IMIN = I 

CONTINUE 

О 17 1 = 1 , NB 

IF (I .LT. IMIN) GOTO 16 
XS(I) = XS(I) - XMIN 

GOTO 17 

ХЕ(І) = XMIN - XS(I) 
CONTINUE 


XMIN = XOS(1) 

DO 20 I = 2 , N1 

IF (XOS(I) .GT. XMIN) GOTO 20 
ХМІМ - Х05(1) 


IMIN = I 

CONTINUE 

PO 22 Г- 1, NI 

E COL EM IMIN) GOTO 21 


XOS(I) = XOS(I) - XMIN 
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INTO1620 
INTO1630 
INTO1640 
INTO1650 
ІХТО1660 
INTO1670 
INTO1680 
INTO1690 
INTO1700 
INTO1710 
INTO1720 
ІКТ01730 
INTO1740 
INTO1750 
INTO1760 
INTO1770 
INTO1780 
INTO1790 
INTO1800 
INTO1810 
INTO1820 
INTO1830 
INTO1840 
INTO1850 
INTO1860 
INTO1870 
INTO1880 
INTO1890 
INTO1900 
INTO1910 
INTO1920 
INTO1930 
INTO1940 
INTO1950 
INTO1960 
INTO1970 
INTO1980 
INTO1990 
INTO2000 
INTO2010 
INTO2020 
INTO2030 
INTO2040 
INTO2050 
INTO2060 
INTO2070 
INTO2080 
INTO2090 
INTO2100 
INTO2110 
INTO2120 
INTO2130 
INTO2140 
INTO2150 
INTO2160 


20 


60022 
KXOSCI) в 
CONTINUE 


CALL DIFF3(NT, NS ӘСЕЛ е 

CALL INTRP3(NT,XS,DLSC,D1,D2,D3,N1,XOS,DLS) 
CALL AMEAN (1,N1,XOS,DLS,1) 

CALL DIFF3 (NT,XS,VNC,D1,D2,D3,0) 

CALL INTRP3(NT,XS,VNC,D1,D2,D3,N1,XOS, VN) 
CALL АМЕАМ (1,М1,Х05,УМ,1) 

CALL DIFF3 (NT,XS,DBPPC,D1,D2,D3,0) 

CALL INTRP3(NT,XS,DBPPC,D1,D2,D3,N1,XOS, DBPP) 
CALL AMEAN (1,N1,XOS,DBPP, 1) 

END IF 

END IF 

DO 30 Iz1,N1 

ХР(Т) = XSTGR(I) 

YP(I) = YSTGR(T) 

CONTINUE 


C CALCULATE GEOMETRIC QUANTITIES 


100 


575) 


С 


DO 100 J=1,N 
МООС) = OST VN J) + UNI 

A (гү, 

MENS EIA AAA 

D(J)= SQRT((XP(J+1)-XP(3))%%2 + (YP(J+1)-YP( J) )**2) 
E(J)= (XBCJ+1)-xeQ т) 

= APTO) 

CONTINUE 


IF ( INLET .NE. O . AND. . NOT) SINGER ШЕ 
SUM = D(1) 

DO SS SS 

SUM = SUM + D(J) 

CONTINUE 

ЕЕ PI = SUM aor 

Boor 

Q = 0.0 

ЕКО ТЕ 


C CALCULATE NORMAL AND TANGENTIAL MATRICES 


102 


105 


CONTINUE 
IF (SINGLE) THEN 

IF ( .NOT. OFF) THEN 

DO 120 I=1,N 

DO 110 J=1,N 

ТЕ (в. 1 058 0/41/25 

XX» (X(1)-X(J))*0(J) + (Ү(І)-ҮСЈ))*8(7) 
YY=-(X(I)-X(J))*S(J) + (Y(I)-Y(J))#C(J) 

UU= LOG(((XX+. 5%D(J))**2+YY**2)/((XX-. 5*D(J))**2FYY**2)) 
үү= 2. *ATAN2(YY*D(J), XX 2yr O 2) 

SS= S(I OCD) - C(I)*S(J) 

СС= C(1)*C(J) + S(1)*S(J) 

A(I,J)= -UU*SS + VV*CC 

B(I,J)= UU*CC + VV*SS 

COTO JIO 

А(Т,Ј) = 6.2831853 
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170246 
INTO24% 
ІМТ0248І 
INTO249 
INTO253 
INT0251 
INT0252 
INT0258 
INT0254 
INTO253 
INTO254 
INTO25] 
INTO2 5m 
INTO256 
INTO26€ 
INTO26] 
INTO262 
INTO263 
INTO264 
INTO265 
INTO266 
INTO267 
INTO268 
INTO269 
INTOZUM 
INTO27] 
INTO274 


13 


0 


50 


45 


40 


ЕЕ 


+ 


+ 


ЩЕ 


Bio) - 0.9 
CONTINUE 
CONTINUE 


DO 130 J=1,N 

ХХ= (XOFF(1)-X(J))*C(J) * (YOFFCI)-Y(J))*S(J) 
YY--(XOFF(I)-X(J))*S(J) + (YOFF(I)-Y(J))*C(J) 

UU= LOG(((XX+. 5*D(J))* — үзе 2) /((XX=. 5*D( J) )*2+YY#*2) ) 
VV2 2. *ATAN2(YY*D(J), XXv*24YY8*2-(. S*D(J) )992) 

SS- SOFF(I)*C(J) - ¿s D 

ES-OGOERCT)SO(J) + SOFECI)*SCJ) 

ADFF(I,J)s -UU*SS + VV*CC 

BOFF(I,J)= UU*CC + VV*SS 

CONTINUE 
CONTINUE 
END IF 

ELSE 

WU . NOT. 
DO 50 I=1,N 
DO 40 J=1,N 
MET EQ. I) GO TO"45 

В XDD CD (YCI)-Y(J))*S(J) 
E. meT) SCI) + x С) 


OFF) THEN 


Еј = р 
| = *CY(: )-YP( )) / E 
NO = и *(X(1)-XP(J+1)) / SP 


OO (OD БР 


R1SQ 2 (COSH(DZXi))**2 - (COS(DY1))**2 

Е250 = (COSH(DX2))**2 - (COS(DY2))**2 

UU = LOG(R1SQ/R2SQ) 

NUM1 = DX1 * COSH(DX1) * SIN(DY1) - 
DY1 * SINH(DX1) * COS(DY1) 


n 


DNUM12 DX1 * SINH(DX1) * COS(DY1) + 


Dri GoOSHODXISESS SINCDY) 
пеша - рх2 5 COSHCDX2)95S SIN(DY2) - 
ӘЛ ОУ х2 57 1606068 


DNUM2= DX2 * SINH(DX2) * COS(DY2) + 


Di Ze О CDA 2 eer oN DY 2) 
EXV = 2.0 * ATAN2(NUM2,DNUM2) - 2.0 * ATAN2(NUM1,DNUM1) 
Wee ТАМОС ОС Л) КА 2t Y Ysw»-U. S*D( J))v*2) 


VV = VV + EXV 

SS= S(])#C(J) - C(1)*S(J) 
ВЕТО) СОЈ) + 5(1)*5(7) 
A(I,J)= -UU*SS + VV*CC 
И те UU*CC 4 VV*SS 


67228901853 
OO 


DO 70 I=1,N 
DO 60 J=1,N 
XX= (XOFF(1)-X(J))*C(J) + (YOFF(1)-Y(J))*S(J) 
YY=-(XOFF(1)-X(J))*S(J) + (YOFF(1)-¥(J))*C(J) 
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INTO2730 
INTO2740 
1002750 
INTO2760 
02770 
02760 
INT02790 
INT02800 
INTO2810 
Ji 102850 
INT02830 
INTO2840 
INTO2850 
INTO2860 
INTO2670 
INTO 2880 
INTO2890 
1102300 
INTO2910 
189102920 
15:802 9:90) 
INTO2940 
INT02950 
1NTO2960 
1 102970 
INTO2980 
{102990 
INTO3000 
INTO3010 
INT03020 
INTO3030 
INTO3040 
INTO3050 
ІМТ03060 
INTO3070 
INTO3080 
005 090 
INTO3100 
INTOS IIO 
J 105120 
INTO3130 
INTO3140 
103150 
INTO3160 
INTO3170 
INTO3180 
Ji 1035150 
INT03200 
ІМТ03210 
INTO3220 
INTO3230 
INTO3240 
INTO3250 
1130005260 
INTO3270 
Перо во 


60 
70 


+ 


+ 


+ 


T 


DX1 = PI *(XOFF(I)-XP(J)) / SP 

DY1 = PI *(YOFF(I)-YP(J)) / SP 

DX2 = PI *(XOFF(I) РОО E 

DY2 = PI *(YOFF(I)-YP(J41)) / S 

R1SQ = (COSH(DX1))**2 - (pu И 

R2SQ 2 (COSH(DX2))**2 - (COS(DY2))**2 

UU = LOG(R1SQ/R28Q) 

NUM1 = DX1 * COSH(DX1) * SIN(DY1) - 
DY1 * SINH(DX1) * COS(DY1) 


DNUM1- DX1 * SINH(DX1) * COS(DY1) + 
DY1 * COSH(DX1) * SIN(DY1) 
NUM2 = Dx2..** СОЗСО ООО 
ру2 * SINHCDRZ) = 0050DY2) 
DNUM2= DX2 * SINH(DX2) * COS(DY2) + 
DY2 * COSH(DX2) * SINCDY2) 
EXV = 2.0 * ATAN2(NUM2,DNUM2) - 2.0 * ATAN2(NUM1,DNUM1) 
УМЕ 2, *ATAN2(YY*D(J), XX"*2-4YY**2-(.5*D(J))**2) 
VV = VV + EXV 
55- SOFF(I)*C(J) = COPE GIS J) 
CC- COFF(I)*C(J) * SOFF(I)*S(J) 
AOFF(I,J)2 -UU*SS * VV*CC 
BOFF(I,J)- UU*CG + VV*SS 
CONTINUE 
CONTINUE 
END IF 
END IF 


C NORMAL AND TANGENTIAL COMPONENTS OF FUNDAMENTAL SOLUTIONS 


C 


145 


IESU 


160 


DO 160 I=1,N 

SUMR= 0. 

SUMT= 0. 

ТЕ ( „МОТ. OFF) THEN 
R(1,1)= SCI)+VNC(1)/COS( ALPHA) 
То SR (I) 
Dose 
TZ) = 5009 

DO 145 J=1,N 

SUMR = SUMR + B(I,J) 
SUMT = SUMT + A(1,J) 


CONTINUE 

ELSE 

R(1,1) = SOFF(1) 
T(1,1) = COFF(I) 
R(2,1) =-COFF(I) 
Т(2,1) = SOFF(I) 


DO 150 J=1,N 

SUMR= SUMR + BOFF(I,J) 
SUMT= SUMT + AOFF(1,J) 
CONTINUE 

END IF 

R(3,1)= SUMR 

ІЗ ШЕ ШІ 

CONTINUE 


ТЕ С ОЕР) Со TORTE 
DECOMPOSITION OF MATRIX A 
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INTO339f 
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INTO3420 
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INTO3450 
INTO3460 
INTO34 7 
INTO348 
INTO3490 
ІМТ03500 
INTO3510 
INTO35@ 
INTO353 
ІМТ0354 
ІМТ0355 
ІМТ0356 
INTO35A 
INTO3581 
INTO35Y 
INTO360 
ІХТ0361 
INTO3626 
INTO363 
INTO364 
ІМТ0365 
ІМТ0366 
ІМТ0367 
INTO368 
INTO369 
ІМТ0370 
INTO371 
INTO378 
INT0373 
1№Г0374 
INT0378 
INTO376 
INTO377 
INTO378 
ІМТ0379 
ІМТ0380 
ІМТ0381 
ІМТ0382 
ІМТ0383 
ІХТ0384 






210 
220 
230 
ШІ OPE 


Бэр 
260 
270 
C BAC 


280 
EOD 
300 


A 
D 


СС СУ С?З 


305 


306 


8 


507 


C CAL 
C 
279 


DO 230 I=1,N-1 

DO 220 K=I1+1,N 

ПОШ АСЕ ІУ)/АС(Т 1) 

DO 210 J=1+1,N 

ACK, J)= A(K,J)- АСК, "АСТ, У) 
CONTINUE 

CONTINUE 

CONTINUE 

RATE ON FUNDAMENTAL RIGHT SIDES WITH LOWER TRIANGULAR 
ШІ 270 K=1,3 

DO 260 J=1,N-1 

DO 250 I=J+1,N 

ЕШЕІ)- ВКТ) - АСТ, ВСК, Г) 
CONTINUE 

CONTINUE 

CONTINUE 

K SOLUTION 

Ш 500 к=1 3 

DO 290 I=N,1,-1 


SUM= 0. 

DO 280 J=N,I+1,-1 

БЕЛЕ сым 1 ACJ SIGUN) 

CONTINUE 

eet k,1)= (R(K, D) -SOMIWACI,D) 

CONTINUE 

CONTINUE 

OFF = „TRUE. 

DD DIS-PLACE VERTICALLY TO THE BODY TO GENERATE 
ISPLACEMENT SURFACE 

DO 305 I=2 , N 

AER) = (S(1)“D(1-1)+S(1-1)*D(1))/CD(1-1)+D(1)) 
COFF(I) = (C(I)*D(I-1)+C(I-1)*D(1))/(D(I-1)+D(I)) 
CONTINUE 

SOFF(1) = 2.0*S(1) - SOFF(2) 

SOFF(Ni)= 2. 0*S(N) - SOFF(N) 

COFF(1) = 2.0*C(1) - COFF(2) 

COFF(N1)= 2.0*C(N) - COFF(N) 

be 306 1 = 1, Ni 

XPOFF(I) = XP(I) - SOFF(I) * DLS(I) 

YPOFF(I) = YP(1) + COFF(I) * DLSCI) 

CONTINUE 

lD 3071 =1 , N 


NEEP(I) - E Р * (XPOFF(I) + XPOFF(I+1)) 
YOFF(I) = * (YPOFF(1) + YPOFF(1+1)) 
DOFF = АНИ г e i 
(XPOFF(I-1)-YPOFF(I))**2 ) 
(XPOFF( I*1)-XPOFF(I))/DOFF 
(YPOFY(1+1) -YPOFF(1))/DOFF 


COFF(I) 
SOFF(I) 
CONTINUE 
G TO 102 

CULATION OF SURFACE VELOCITIES FOR THE FUNDAMENTAL SOLUTIONS 


DO 330 K=1,3 
DO 320 I=1,N 
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15703850 
ІХТ03860 
ОВОО 
1703880 
Т5 690 
1М102900 
INDOSOID 
и 2928 
IPSO OU 
INTO3940 
Jl O 
ІМТ03960 
1М103970 
1103980 
1610 990 
INTO4000 
INTO4010 
INTO4020 
INTO4030 
INTOA4040 
INTOA050 
INTO4060 
ІМТ04070 
INTO4080 
INTO4090 
INTO4100 
INTO4110 
INTO4120 
INT04130 
ІМТ04140 
INTO4150 
INTO4160 
ІМТ04170 
INTO4180 
INT04190 
INTO4200 
INTO4210 
INTO4220 
INTO4230 
INTO4240 
INT04250 
15104260 
INTO4270 
INTO4280 
INTO4290 
INT04300 
15704310 
INTO4320 
INT04330 
INTO4340 
INTO4350 
1NTO4360 
INTO4370 
INTO4380 
INTO4390 
INTO4400 


SUMT- T(K,I) 
SUMNZ-R(E,I) 
DO 310 J=1,N 
SUMT= SUMT + BOFF(1,J)*SIG(K,J) 
SUMN= SUMN + AOFF(1,J)*SIG(K,J) 
310 CONTINUE 
VTAN(K,1)= SUMT 
VNOR(K,1)= SUMN 
320 CONTINUE 


DOFF1 = SQRT((XPOFF(2)=XPOFF(1))2+( YPOFF(2)-YPOFF(1))%*2) 
DOFF2 = SQRT((XPOFF(3)-XPOEF(2))*“"*2+(YPOEECS) =1P OB УЛ 


DOFFN SQRT( (XPOFF(N+1) -XPOFF(N) )***2+( YPOFF(N+1) -YPOFF(N)) 

- 
2) 

DOFFN1= SQRTCCXPOFFQN) -XPOFF(N-1) )**2+( YPOFF(N)-YPOFF(N-1)) 
+ xn BD) 

VUT(K) = VTANCK, М) + DOFFN * (VTAN(K,N)-VTAN(K,N-1))/ 
+ ( DOFFN+DOFFEN1) 

VUN(K) = VNOR(K,N) + DOFFN * (VNOR(K,N)-VNOR(K,N-1))/ 
+ ( DOFFN+DOFEN1) 

VLT(K) = VTAN(K,1) + DOFFi * (VTAN(K,1)-VTAN(K,2))/ 
+ (DOFF1+DOFF2) 

VLN(K) = VNOR(K,1) + DOFF1 * (VNOR(K,1)-VNOR(K,2))/ 
+ (DOFF1+DOFF2) 


330 CONTINUE 
C OUTPUT FUNDAMENTAL SOLUTIONS 
IF (ICYCLE .EQ. 1 .OR. IGNGRE .GE. IGXEL-1 .OR- ІР ВЕ 0) 

+ WRITE(6,335) TITLE 
335 FORMAT(1H1,///20A4//) 
C DO 360 K=1,3 
C WRITE (6,340) K 
C340  FORMAT(////,1H ,' FUNDAMENTAL SOLUTION NUMBER "„12////) 
E WRITE(6,345) 
345 FORMAT(3X,'I',8X, X' , 11X, 'Y' , 10X, ' VT' , 10X, ' VN' ,8X, ' SIG' // m 


C WRITE(6,375) 1, XP(1) , YP(1), VLT(K),VLN(K) 

C DO 350 I=1,N 

Е WRITE(6,375) 1, XCD), Ү(І), VTANCK, D , VNORCE, D) ,SIG(K, I) 
C350 CONTINUE 

C WRITE(6,375) N1 , XP(N1),YP(N1),VUT(K) , VUN(K) 


C360 CONTINUE 
C COMBINED FLOW AT ANGLE OF ATTACK 
e 
IF ( INLET .NE. 0) THEN 
YYY = ((VUT(3)+VLT( 3) )*TAN( ALPHAI )+( VUT( 1)+VLT( 1) )*Q)/ 
+ ( (VUT(3)+VLT( 3) )-(VUT(2)+VLT( 2) )*Q) 
XXX = -((VUT(1)+VLT( 1) )+( VUT(2)+VLT(2) )*TAN( ALPHAT)) / 
+ (( VUT(3)+VLT( 3) ) -( VUT(2)4+VLT(2))*Q) 
ALPHA = ACOS(i. 0/SQRT( 1. O+YYY**2) ) 
COSAL = COS( ALPHA) 
SINAL = SIN( ALPHA) 
W = XXX/SQRT(1. O+YYY**2) 
ELSE 
COSAL = COS( ALPHA) 
SINAL = SIN( ALPHA) 
W=-((VLT( 1)+VUT( 1) )*COSAL+( VLT( 2)+VUT(2))*SINAL) / 
+ CVLTC3) OUR 695) 


INTO4410 





INT04960 


END IF INTO4970 

C FORCE COEFFICIENT CALCULATION INTO4980 
SEMIS 0: INTO4990 

Sx o INTO5000 

EEMY- 0: INTO5010 

DO 390 I=1,N ІКТ05020 

SUM1= SUM1+ D(I) INTO5030 

SUMX2 SUMX- VCOM(I)**2*S(I)*D(I) INTO5040 

SUMY- SUMY- VCOM(I)"**2*C(I)*D(I) INTO5050 

390 CONTINUE INTO5060 
C FIND MAN. CHORD LENGTH INT05070 
XOMIN + ХО(1) INTO5080 

DO 395 I=2 , Nl INTO5090 

IF ( XO(I) .GT. XOMIN) GOTO 395 INT05100 

XOMIN = XO(I) INTO5110 

395 CONTINUE INTO5120 
CHORD = XO(N1) - XOMIN INT05130 

CL1- SUM1*25. 13274*W/CHORD INT05140 

CL2= (SUMY*COSAL-SUMX*SINAL) /CHORD INTO5150° 

CD = (SUMX*CCSAL+SUMNY*SINAL) /CHORD INTO5160 

E INTO5170 
С CALCULATING PARAMETERS FOR INLET VELOCITY AS MODULUS OF NOMORIZED VELINTO5180 
C INT05190 
IF (.NOT. SINGLE) THEN INTO5200 

NUM1 = SIN( ALPHA)+CL1*CHORD/(4. 0*SP) INTO5210 
ALPHID = ATAN2(NUM1,COS( ALPHA) ) S220 

NUM1 2 SIN(ALPHA)-CL1*CHORD/(4. O*SP) INT05230 
ALPHED = ATAN2(NUM1,COS(ALPHA)) INTO5240 

NUM1 = CL1*CHORD/(2.0*SP)*COS( ALPHA) INTO5250 
DNUM1= 1. 0-(CL1*CHORD/(4. 0*SP))**2 INTO5260 
DALPHA = ATAN2(NUM1,DNUM1) INTO5270 

UOUI 2 (TAN(ALPHID)-TAN(ALPHED))*(2. O*SP/CHORD*COS(ALPHID))/CL1 INTO5280 

BET — CLI * UOUI*-2 1№Т05290 

ПЕСО = № 0/U0UL INT05300 
VEXIT = COS(ALPHA)/COS(ALPHED) INTO5310 

ELSE INTO5320 
ALPHID = ALPHA INTO5330 
ALPHED = ALFHA INTO5340 
DALPHA = 0.0 које 

ШЕШІ - 10 ІКТ05360 

MOU = 1 0 1405370 

ШЕГІ = С INTO5380 
WEXIT = 1.0 INTO5390 

END IF INTO5400 

ЕЕ. 180. О/РТ INTO5410 
ALPHID = ALPHID * FAC INTO5420 
ALPHED = ALPHED * FAC INTO5430 
DALPHA = DALPHA * FAC INTO5440 
ALPHAD = ALPHA * FAC INTO5450 

ШО (ШЕСЕ ЖО. 1 ОВ ТОМСЕН СЕИ ICYTL-1 .OR. IP .GE. 0) THEN  INTOS460 
WRITE(6,370) ALPHAD ,ALPHID,ALPHED,DALPHA,UIOU,VEXIT INTO5470 

370 FORMAT (////,1H , ‘COMBINED FLOW AT AVERAGE ANGLE OF ATTACK = ', INTO5480 
+ F8.3, 4X,'DEGREES', /,1H ,17X,' INLET ANGLE OF ', INTO5490 

+ "ATTACK = 5 \F6.394%, DEGREES ,/,1H , INTO5500 

+ 17X, EXIT ANGLE =  ',F8.3,4X, DEGREES ',/,1H ,17X, INTO5510 

+ "TURNNING «ANGLE = *,F8.3,4X,'DEGREBS',/,1H ,17X, ІКТ05520 


D 


905 


380 


374 


375 


385 


400 


420 


430 


ec c) 6j 


60 СС) 62 С) сл Со С 


ФФ 


+ "INLET VEL = *,F10%6,3X, EXE Т 
WRITE(6,365) 

FORMAT(3X, '1',8X, ' XO' , 10X, У0' , 10X, x! LAY p: 
+ ОХ, УТ ое .9X, "OPI! ///) 

END IF 

ПО Вот 

VTCOM(I)= VTAN(1,1)**COSAL+VTAN( 2, 1)*SINAL+W*VTAN( 3,1) 
VNCOM(1)= VNOR(1,1)*COSAL+VNOR(2,1)*SINAL+W*VNOR( 3,1) 
VCOM(1)= SQRT(VTCOM(I)**2 + VNCOM(1)***2) 

ТЕ (УТСОМТ) „ПТ. 0.0) COM) NES UID) 

СРСІ) - 1.0 = МЕТ p. 

СРІ(І)т 1.0 - (УСОМ(ІУ%4/001)%%2 

SIGCOM(I) = SIG(1,1)*COSAL+SIG(2,1)*SINAL+W*SIG(3,1) 
ХР(Т) = 0.5 во) 

YPCI) = 0.5 “MMM Ac 

CONTINUE 
IF (ICYCLE .EQ. 1 .OR. ICYCLE .GE. ICYTL-1 .QR. IP GE. OD IPN 
WRITE (1,374) ( XOC1),YOC1),XPCl eect ОР Ту ОР S aE 
WRITE (6,375) ( І, ХОФ ХОСТ); ХЕБ) ЕОНИ ГО 
+ VNCOMCI) , COM EPOD АЕ) пне 
FORMAT(6F10. 4) 

FORMAT(1X, I3, 9F12.4) 


WRITE (2) I,XOCI), YOCI) ,XSTGRCI) , YSTGR( I) ,DLSCI) ,XC I) , YCI) , VCOMC I) INTOS75f 


WRITE(6, 385) N41 XO(N-DD, YO(N+1) 
FORMAT( 1X, ов 212. 4) 
МЕТТЕСС 400) CHORD TOTE CLOI 

FORMAT(///3X, ' CHORD = ' ,Flüsh АХ. "ОШСАШЕ) = F10. 5 N, 
т "CLCINLET) s ',F10.5) 

END “РЕ 

ҒОҺМАТ( /3Х5,1НІ,6Х,2Н50,10Х,2НҢ5С,9Х,3Н УЫ,9Х/ЗНУКС,9Х ӘНЕ 
+ 4HDLSC) 

FORMAT(15,6E12.4) 

RETURN 

END 


DATA SET KCBCAMEAN AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCAMEAN AT LEVEL 003 AS OF 04/05/84 
SUBROUTINE ANEANCNS,ND,X,Y, IT) 


SMOOTH DATA USING 3-PTS WEIGHTING FORMULA 

NS : STARTING PINT OF DHE PACA LO BESO AED 

ND : END PINT OF THE WATA TO BE МОШЕ 

X, Y : INDEPENDENT + DEPENDENT VARAIBLES OF THE DATA 
TO BE SMOOTHED 

Lh : CYCLES OF DATA SMOOTHING 


DIMENSION X(101),Y(101) 


NM = ND -NS 
ІЕСММ МЕ. 22 ОС. TT Tl l SE OMS 
NDM1 = ND - 1 
МР? ЕР 
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INTO5530 
INTOS54 

INTOS550 
INTO5560 
INTO5570 
INTO5580 
INTO5590 
INTO5600 








INTO5 76€ 


ІКТ05780 
INTOS790 
INTO5800 
INTO5810 
INTO5820 
INTO5830 
INTO5840 
INTOS850 
INTOS86 

INTOS87 

ІКТ05880 
INT05890 
INTOS90 

INTOS91 

INTO5920 
INTO5930 
INTOS940 
INT0595 
INT05960 
INT05970 
INT05980 
INTO5990 
INTO6000 
INTO6010 
INTO6020 
INTO6030 
INTO6040 
INTO6050 
INTO6060 
INTO6070 
INTO6080 


10 


aa aa a 


aaa 


10 


с С O Г 


C3 


DO 20 K-1,IT 


DL1 = X(NSP1) - XCNS) 
Yl = Y(NS) 

DO 10 I=NSP1,NDM1 

DENEN j 1 ]) -X(1) 
у2 = ¥(I) 

YM | - (DL2 - Yl + DL1 * Y(I+1))/(DL1 + DL2) 
Ү(І) = 0.5 * (Y2 + YM) 
DLI = DL2 

її = ү? 

CONTINUE 

CONTINUE 

RETURN 

END 


DATA SET KCBCBLGRID AT LEVEL 001 AS OF 08/24/84 

DATA SET KCBCBLGRID AT LEVEL 001 AS OF 08/24/84 

DATA SET KCBCBLGRID AT LEVEL 004 AS OF 04/05/84 
SUBROUTINE BLGRID(N,X,T,D1) 


GENERATE В. L. X-WISE GRID USING MODIFIED COSINE DISTRIBUTION 


s S ION X(101).T(101),D1(101) 
IDE RADUS7. 29577957, BPI/3. 14159265/ 


NN т 

EN = FLOAT((NN-1)/2) 
THO = 10. /CRAD 

CTOl = 1. + COS(THO) 
DTH = (ВРІ - THO) / EN 
БИ = FLOAT(N - 2) 

DO 10 I=N,NN 

FI = 1.0 + FI 

II =] - N+ 1 

XII = THO + FI * DTH 
Maly = (1.0 cr1608(XII))WeTOl 
CONTINUE 

У1 = X(1) 

XN = X(N) 

CH = XN -X1 

ЕМІ = FLOAT(N-1) 

N10 = №/10 

DO 20 I=1,N 

Т(1) = FLOAT(1-1)/FN1 
ХХІ) = (X(1I)-X1)/CH 
CONTINUE 


CALL SMFIT(N10,N,T,X,D1,N10) 
ШЫЛ) ІЛ. 0 35 * X(3)) X(2) 9 0.35 *.X63) 
DONNE SMETTCISN,T,X,D1,2) 

CALL AMEAN(1,N,T,X,N10) 


RETURN 
END 


DATA SET KCBCBL2D AT LEVEL 001 AS OF 08/24/84 


INT06090 
INTO6100 
| таб 110 
INTO6120 
INTO6130 
INTO6140 
INT06150 
INT06160 
INT06170 
INT06180 
INTO6190 
INTO6200 
1106210 
1106220 
INT06230 
INTO6240 
ІМТ06250 
INTO6260 
INTO6270 
INTO6280 
INTO6290 
INTO6300 
INTO6310 
INTO6 320 
1359909350 
INT06340 
INTO6350 
INTO6360 
INTOG370 
INT06380 
1406390 
INTO6400 
INTO6410 
INTO6420 
INTO6430 
INTO6440 
INTO6450 
INTO6460 
ІМТ06270 
INTO6480 
INTO6490 
INTO6500 
1106510 
0520 
INTOGCS30 
INT06540 
INT06550 
INT06560 
ІЧТОӨ570 
INTO65380 
INTO6590 
INTO6600 
INTO6610 
INT06620 
INT06630 
INTO6640 


Сз СУ С? 


О СУС C) C9 


ИЕ С) 


РС) О) 


20 


30 


DATA SET KCBCBL2D 
DaTA SET KCBCBL2D 


SUBROUTI 


COMMON / 
COMMON / 
COMMON / 
COMMON / 
COMMON/E 


¥ 


COMMON / 


+ 


COMMON / 
COMMON / 
COMMON / 
COMHON / 


J 


COMMON / 
COMMON/P 
DIMENSIO 


DO 5 I = 
ALFAS(I) 
FFS(I) 
RTS(I) 
CONTINUE 
CALL INP 


WRITE(6, 


AT LEVEL 001 AS OF 08/24/84 
AT LEVEL 012 AS OF 04/06/84 
SE EL2D ( ТТК, ТОРТ ОМК 


PROGRAM CALCULATES VISCOUS/INVISCID INTERACTION USING HILBERT 
INTEGRAL. 


BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 

BLC1/ F(101,2),U(101,2) ,V( 1013205 W 2) ВОТОТ 2) 

BLC2/ DELF(101),DELU(101),DELV(101) ,DELW( 101) 

BLC7/ C(100,100),D(100) ,DB( 100) ,DBP( 100) ,UEO( 100) ,GI 

DDY1/ RL,RX,SQRX,RXNTR,GMTR , GMTRS( 100) , ALFAS( 100) , 
FFS(100),RTS(100),IEDY,NXSPT 

SMRY/ VW(100),ITP(100),ISL(100),DLS(100),CF( 100), 
THT( 100) , NPSTR( 100) 

GTY / X(101),UE(100),P1( 100) ,P2( 100) , CEL, CELH 

BONV/ ITMAX,EPSL,EPST,CONV 

SAVE/ FS(101),US(101),VS$(101),WS( 101) ,BS( 101) 

BLIN/ TITLE(20),XC(100),YC(100),ISG( 100) , DELS( 100) , 
XCTR,XTR,ISTRP, ICYCLE, ICYTL_XCTRS(2), TREINDOD) 

1SURF/ ISF 

ПОЛА сар МУРА O 

№ SURFID(4) 


GENERATE B. 1. 6405 + SET INLITAL CONDEMIONS 


XT 


Wo Wee 


SO O = 


0. 
i; 
Je 


UNITE ISWPT URET) 


CALCULATE HILBERT COESBS. 097) 


CALL (GEG TIONS TO) 
LOOP OF CALCULATIONS 
NSS = NS 
NXSPT = NXT + 1 
IF ( ICYCLE .EQ. 1 ) NS = NXT + 1 
NX =NX+ 1 
CEL -0.5 % (Х(МХУ + X(NX-1)) /(X(NX) -X(NX-1)) 
CELH = 0.5 * CEL 
IT = 0 
RX = UE(NX)**X(NK)*RL 
SQRX = SQRT(RX) 
ІТ =IT+1 
IF(IT .LE. ITMAX) GO TO 40 
NXM1 = NX-1 
CALL HEADER( TITLE,SURFID,ISTRP ) 
WRITE(6, 170 ) (M,X(M),CF(M) ,DLS(M),UE(M),P2  (M),THT(M), 


D(M) ,ALFAS(M), 
160 ) NX 


ITP(M) ,NPSTR(M) , M21, NXM1) 


76 


|| 
| 


INTO6650 
INTO6660 
INTO6670 
INTO6680 
INTO6690 
INTO6700 
INTO6710 
INTO6 720 
INTO6730 
INTO674( 
INTO6758 
ІМТ06760 
INTO6778 
INT06780 
INT06790 
INT0680( 
ІМТ06810 
INTO682¢@ 
ІМТ06830 
ІКТ06840 
ІКТ06851 
INTO686( 
INTO687( 
INTO688( 
INTO689( 
ІХТ0690( 
INTO691( 
ІМТ0692( 
ІМТ0693( 
INTO694( 
ІМТ0695( 
INTO696( 
INTO697( 
INTO698( 
INTO699( 
INTO7 006 
INTO701¢ 
INTO702] 
INTO703( 
INTO704( 
INTO7051 
INTO706( 
INTO/O@ 
INTO7080 
ІМТ0709( 
INTO7 108 
INTO711g 
INTO7 12€ 
INTO713g 
INTO7 14€ 
INTO715€ 
INTO7160 
INTO717( 
INTO7 18€ 
ІМТ0719% 
INTO7 20M 


ЗБОР 
40 CONTINUE 
ВЕС. NTR) GALL EDDY 
CALL COEFTR 
CALL SOLV3 
i (1520 0.0) СОТО 60 


EXTRAPOLATE CALCULATED D FOR TURBULENT SEPARATION OR LAMINAR 
SEPARATION FOR LAMINAR FLOW CALCULATION ONLY 


Cet Су 


CALL EXTRAP(NX,NXT,X,D) 

NKM1 = NX - 1 

CALL HEADER( TITLE,SURFID,ISTRP ) 

WRITE(6, 170 ) (M,X(M),CF(M),DLS(M) ,UE(M) ,P2(M) , THT(M), 


т D(M) , ALFAS(M) , ITP(M) ,NPSTR(M) , M=1,NXM1) 
WRITE(6, 180) 

WRITE(6,190) (M,X(M),D(M) ,M=NX,NXT) 

GOTO 130 


60 IRNS GT. NTKONGSERO- 70 
SEN SCDELV(1)) . GT. EPSL) GO TQ 30 
ЕШ ТО 80 
70 CONTINUE 
G S UELV(1)/VCr,2Z)) „СТ. EPST) GO TO 30 


80 CONTINUE 

C 

C CHECK FOR GROWTH 

C 
ПМР - GE. NP.) GO 10.90 
IAS VEND, 27) ТОО 0005 . AND. ABSC1. 0-UC(NP-2,2)) 

> "LIT 0000 SMS OO 

CALL FILLUP(1) 
=. | 
Cree 30 

90 CONTINUE 

C 


CALL FILLUP( 2) 
CALL OUTPUT(1) 
ШОЕ. ЕО. 0 ОЕ. МЕЕ МТК) 6ЮТО 100 
Ше ОЕТ. 3 .ОК. ITR. NE. 3) 69500700 
CALL TRNS(ICODE) 
ЦЕО ОВОЈ ОИ РО О 20) 

100 MERX -NE. N55) GOTO 120 


C 
C | STORE PROFILES AT THE STATION NS FOR INVERSE В. Г. 
C | CALCULATION 
C 
DeEIOJ- 1, NPT 
FS@J) = F(J,2) 
US(J) = U(J,2) 
Wes) = v(J,2) 
WS(J) = W(J,2) 
P J) S BC 3,2) 


110 CONTINUE 

120 IIS: LT. NSS ) GOTO 10 

C REGIO COLE NE. 1160 IO Ко 
ПАЗИ "СЕ. КУ) 6970 130 


p 


J 407210 
107220 
1 107230 
INTO72%0 
J 167250 
INTO7260 
INTO7270 
INTO7280 
INTO? 290 
INT07300 
PINTOS HO 
INT07520 
INTO7330 
INTO7340 
INTO7350 
1897 560 
14107 70 
INT07380 
INTO7390 
INTO7400 
INTO7410 
INTO7420 
INTO7430 
INTO7440 
INTO7450 
INTO7460 
INTO7470 
INTO7480 
INTO7490 
INTO/500 
INT07510 
INTO7520 
INTO7/530 
INTO7540 
[5107550 
INT07560 
INTO7570 
INTO7580 
INT07590 
INT07600 
INT07610 
INT07620 
INT07630 
INTO7640 
INTO7650 
INTO7660 
INTO7670 
INTO7680 
INTO7690 
INTO7700 
INTO7710 
INTOTT20 
INTO7730 
INTO7740 
INTO7750 
INTO7760 


130 
140 


150 
153 


C2 C2 C2 


ES 


Б С) С С 


170 


180 


NOU 


moo. 


COCO 


IF (МХ ТАШ ОЛУ СОТТОЛО INTO777 
CALL HEADER( TITLE,SURFID,ISTRP ) INTO7 78 
WRITE(6, 170 ) (M,X(M),CF(M),DLS(M),UE(M),P2 (M),THT(M), INTO77£ 
T D(M), ALFAS(M), ITP(M),NPSTR(M) ,M21,NXT) INTO78€ 
DORLOTI в INTO781 
DB(I) 9» D(I) ІМТ078? 
NS = NSS INTO782 
NX = NS ІМТ0784 
NP = NPSTR(NX) INTO785 
DO 150 Ј = 1 , NPT INTO786 
F(J,2) = FS(J) ІМТ0787 
(7,2) - 05(7) INTO788 
Үзу Е 9) a 
W(J,2) = WSCJ) INTO79( 
В(7,2) - В5(7) INTO791 
CONTINUE INTO792 
INVRS = NS + 1 INTO793 
INTO794 
CALCULATION SHIFTS TO USING PHYSICAL COORDINATES INTO795 
CALL MAIN2(ITR,ISVPT,SURFID) INTO79f 
1NTO797 
PASS DELTA-STAR BACK TO MAIN PROG. INTO798 
INTO799 
DO 158 I = 1,NXT INTO80( 
DELS(I) = DLS(I) ІХТ080: 
CONTINUE ІХТ0802 
RETURN ІКТ0803 
ІМТ0804 
---------------------------------------------------------------- ІХТ0805І 
INTO806 
FORMAT(1HO,' *** ITERATIONS EXCEEDED ITMAX AT NX = TS INTO807 
+ 1H ,' ** CALCULATIONS STOP. «и INTO80€ 
FORMAT(1HO,' *** SUMMARY OF STANDARD B. L. SOLUTIONS. #%'/ INTO80€ 
+ 1HO ,4X,2HNX,7X, 1HX,12X, 2HCF,11X,3HDLS , 12k , 2 НИ INTO81C 
+ 12X,2HP2,11X,3HTHT, 13X,1HD,10X,4HALFA, 6X, 2HIT, 2X, 2HNey ee 
+ (1H ,3X (13,610. 5, 2k. fee ee INTO812 
FORMAT(1HO,34H FLOW SEPARATES. D IS EXTRAPOLATED/ INTO813 
+ 1H0,3X,3H NX,7X,1HX, 1395 1HB7) INTO814 
FORMATCTH ,3X T3_F10.5,2X Ela 59) INTO815 
END INTO816 
DATA SET KCBCCALCIJ AT LEVEL 001 AS OF 08/24/84 INT0817 
DATA SET KCBCCALCIJ AT LEVEL 001 AS OF 08/24/84 INTO818 
DATA CET KCBCCALCIJ AW LEVEL 005 AS OF 04/05734 INTO819 
SUBROUTINE CALCIJ ( IL, L0) INTO820 
INTO821 
CALCULATE HILBERT INTEGRAL COEFFS INTO822 
INT0823 
COMMON /BLC7/ C(100,100),D(100),DB( 100) ,DBP( 100) ,UEO( 100) ,GI ІМТ0824 
COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR,GMTRS( 100) INTO825 
+ ,ALFAS(100) ,FFS(100) ,RTS( 100) , IEDY,NXSPT INTO826 
COMMON /GTY / X(101),UE(100) ,P1( 100) ,P2(100) ,CEL,CELH ІХТ0827 
DIMENSION E(2) INT0828 
D 
С----------------------------------- INTO830 
INTO831 
PI = 3.14159265 INT0832 


78 


e» cc) 


C» C) CONG) 


+ O O O 


л С 


60 


C9 50903 


СОС Са 


Fal 

111 
РОДОВИ 
E) 

L 
DOGO J 
RII 

K 

Ше 

Dx2 

DX3 
ШЕН J 
mec J 


ШЕКЕ. 1 


E 2) 
GO TO 50 


веко. I 


К1 
F (2) 
597 192510 


ENDO. I 


R1 
в 2) 


CONTINUE 
ШЕТ 1) 
Е(1) 
CONTINUE 
EN 2) 

ҮЛ 

К 

ШІ (71,1) 
CONTINUE 


RETURN 
END 


DATA SET KCBCCOEF 
DATA SET KCBCCOEF 
DATA SET KCBCCOEF 


ze 


PI*SQRT(RL) 
Ер 2-6 

A 

0. 
L 

-i 


C4 C4 S O 
US ES 


TODO 

aly = CK) 
X(K) - X(K-1) 
X(L) - X(K-1) 


ШЕ) со оо 
КОО, (Т CO AO 
ОВ 1+1 


ОВО) АБО АВО р гырха) 


C XCK+1)-X(K) ) / (X(K+1)-X(K-1) ) 
ШЕП 2 41000 АВБ( ОВЕ ( ХСЫ)-Х(К41)) ) ) ғ 2.0) / DX2 


| 


( X(K-1)-X(K-2) ) / ( X(K)-X(K-2) ) 
( R1 * ALOG( ABS( (X(L)-X(K-2)) / DX1 ) ) - 2.0 ) / DX2 


( E(1) - E(2) ) / PI 
(2) 


0. 
[| 

Ки 
E(1) / PI 


AT LEVEL 001 AS OF 08/24/84 
AT LEVEL 001 AS OF 08/24/84 
AT LEVEL 007 AS OF 04/05/84 


SUBROUTINE COEF(GAMMA1 , GAMMA2 ) 


COMMON /BLCO/ 
COMMON /BLC1/ 
COMMON /BLC6/ 
ще 
COMMON /BLC7/ 


CALOULATE COEFFS OF B. L. FINITE-DIEFERENCE EQS. 
SEMI-TRANSF VAR^TBLES( AFTER SWITCHING). 


IN 


NX ,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 

2101,2), (Оо Оо ( 101, 2) ,B(101,2) 

61(101),52(101),53(101),54(101),55(101),56(101), 
S7(101),S8(101),R1( 101) ,R2( 101) , R3( 101) , Rá( 101) 

C( 100,100) ,D( 100) ,DB( 100) ,DBP( 100) ,UEO( 100),GI 


79 


ІМТ08330 
INTO8340 
INTO8350 
INTO8360 
INTO8370 
INTO8380 
INTO8390 
INTO8400 
INTO8410 
INTO8420 
INTO8430 
INTO8440 
INTO8450 
INTOS460 
INTO8470 
INTO8480 
INTO8490 
INTO8500 
INTO8510 
INTO8520 
ІКТ08530 
INTO8540 
INTO8550 
INTO8560 
INTO8570 
INTO8580 
INTO8590 
ІКТ08600 
ІКТ08610 
INT08620 
INTO8630 
INTO8640 
INTO8650 
INTO8660 
INTO8670 
INTO8680 
INTO8690 
INTO8700 
INTO8710 
INTO8720 
INTO8730 
INTO8740 
INTO8750 
INTO8760 
INTO8770 
INTO8780 
INTO8790 
INTO8800 
INTO8810 
INTO8820 
INTO8830 
INTO8840 
INTO8850 
INTO8860 
INTO8870 
INTO8880 


COMMON /GRD / ETAC 101), DETAGI GH ye P) INTO88 


COMMON /GTY / X(101),UE(100),P1(100),P2(100) ,CEL,CELH INTO89( 
C INT089j 
C ————————————— eese o 
C 

РІН = 0.5 * P1(NX) 

DO 100 J= 2,NP 

FLARE = 1.0 

FB = 0.5*(F(J,2) + F(J-1,2)) 

UB B шыш + U(J-1,2)) 

FVB - 0.5*(F(J,2)*V(J,2)4F(J-1,2)*V(J-1,2)) 

IF(UB .LT. 0.0) FLARE = 0.0 

VB = Оо иж 

USB = 0.5%(U(J,2)%%2 + U(T-1,2)**2) 

WSB = 0.5*(W(3,2)**2 + W(J-1,2)**2) 

DERBV =(B(J,2)*V(J,2) - B(J-1,2)*V(J-1,2))/DETA(J-1) 

FB4 = 0.55080, 1) FFU T 5 

VB4 - 0.5%(У(7,1) % У(7-1,1)) 

USB4 = 0.5%(U(J,1)%*2 + U(J-1,1)**2) 

9584 = 0.5#(W(J,1)%%*2 - W(J-1,1)**2) 

FVB4 = 0.5*(F(J,1)*V(J,1)+F(J-1,1)*V(J-1,1)) 

DERBV4 2 (B(J,1)*V(J,1) - B(J-1,1)*V(J-1,1))/DETA(J-1) 

51(Ј) = CELH*(FB - ЕВА) + P1H*F(J,2) + B(J,2)/DETA(J-1) 

S2(J) = CELH*(FB - FB4) + P1H*F(J-1,2) - B(J-1,2)/DETA(J-1) 

S3(J) = CELH*(VB + VB4) + P1H*V(J,2) 

S4(J) = CELH*(VB + VB4) 4 P1H*V(J-1,2) 

SS(J) = -CEL*FLARE*U(J, 2) 

S6(J) = -CEL*FLARE*U(J-1,2) 

5709) СЕВ, 

S8(J) = CEL*W(J-1,2) 
C 

CRB = -DERBV4 + CEL*WSB4 - CEL*FLARE*USB4 - P1(NX)*FVB4 

R2(J) = CRB - (DERBV - CEL*FLARE“USB + CEL*(VB+VB4)*(FB-FB4) + 


+ CEL*WSB 4 P1(NX)*FVB) 
Ri(J) = F(J-1,2) 5 1 2 DE e 





R3(J-1)= U(J-1,2) - UCJ,2) + DEA В 
R4(J-1)= W(J-1,2) - WOJ,2) 
100 CONTINUE 
© 
C BOUNDARY CONDITIONS 
C 
R1(1) = 0.0 
&2(1) = 0.0 
R4(NP) = 0.0 
IF(NX .GE. INVRS) GO TO 120 
САММА1 = 0.0 
GAMMA2 = 1.0 
R3(NP) = 0.0 
RETURN 
120 CONTINUE 
Cl = En NX) * SQRT(X(NX)) 
GAMMA1 = 
GAMMA2 = © s - CII*ETA(NP))/CII INTO941 
R3(NP) = (GI + CII*(ETA(NP)*W(NP,2) - F(NP,2)) -W(NP,2))/CII INTO942 
C INTO943 
RETURN INTO944 


$0 


cc) c) 


СЕЛЕ С 


(op i СС 


100 


END 


SUBRO 


CALCULA 
IN TRAN 


COMMO 
COMMO 
COMMO 


COMMO 
COMMO 


DERBV 


51(Ј) 
SID) 


SEI) 


R1(J) 
R3(J- 
CONTI 
R1(1) 
R2(1) 
R3(NP 
RETUR 
END 


DATA SET KCBCCOEFTR AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCCOEFTR AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCCOEFTR AT LEVEL 004 AS OF 02/21/84 
UTINE COEFTR 


ШЕМЕЗЕКЕ OF B. L. EINITE-DIEEFERENCE EQS. 
SFORMED VARIABLES( BEFORE SWITCHING). 


N /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
Ry ered, 2101 .2).UC101,2),V(101,2),W(101,2) ,B( 107% 2) 
N /BLC6/ $1(101),$2(101) ,$3(101) ,$4(101) ,$5(101) ,$6(101), 


ІМТ09450 
1М709460 
INTO9470 
INTO9480 
INTO9490 
INTO9500 
INT09510 
INT09520 
ҚЫШУ 530 
INTO9540 
INTO9550 
INTO9560 


67(101),58(101),81(101),82(101),Е3(101),Е4(101)15709570 


N /GRD / ETA(101),DETA(101),A(101) 
N /GTY / X(101),UE(100) ,P1( 100) ,P2( 100) , CEL,CELH 


0 Ј= 2,МР 
O EI E ECU) > 
о ито 
о и) 
ОЕЕО бо) 
= (B(J,2)*V(J,2) - B(J-1,2)*V(J-1,2))/DETA(J-1) 
0.5*(F(J,2)*V(J,2) + F(J-1,2)*V(J-1,2)) 
208 F J I (J Fuad. 1)=vC-1.,1)) 
ШТІ С ЕГІП, 
ІТ 
ОБОО рете 2) 
(B(J,1)*V(J,1) - B(J-1,1)*V(J-1,1))/DETA(J-1) 


о о 


4 


CELH*(FB-FB4) + 0. 5%*P1(NX)*F(J,2) + B(J,2)/DETA(J-1) 

CELH*(FB-FB4) + 0.5*P1(NX)*F(J-1,2) -B(J-1,2)/ 
DETA(J-1) 

CELH*(VB + УВА) + 0. 5*PI(NX)*V(J,2) 

CELH*(VB + УВА) + 0. 5*P1(NX)*V(J-1,2) 

-(GEL+P2(NX))*U(J,2) 

-(CEL+P2(NX))*U(J-1,2) 


Bh PICNX=1)=EVB4 - P2(NA-1)*USB4 + P2(NX-1) 

-GIBA= ABE LSUSB% - P2(NX) 

= СЕВ - (DERBY + P1(NX)*FVB- (CEL+P2(NX))*USB + СЕТ“ 
(VB + VB4)*(FB - FB4)) 


= Pi J-1,2) - PO) DEM UE 
ОО) ОЕ т) В 
NUE 

= 0.0 

= 0.0 
) = 0.0 
N 


DATA SET KCBCCOMPBL AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCCOMPBL AT LEVEL 001 AS OF 08/24/84 


81 


1109580 
15109590 
INTO9600 
INTO9610 
Во 
INTO9630 
INTO9640 
ВО 
INTO9660 
INTO9670 
INTO9680 
INTO9690 
IN10 700 
INT09710 
INT09720 
1NTO9730 
INTO9740 
INTO9750 
INTO9760 
IN IO 7 0 
INTO9780 
INTO О О 
INTO9800 
INTO9810 
INTO9820 
ІМТ09830 
INTO9840 
119559 
INTO9860 
INT09870 
INTO9880 
INDOSSSO 
INTO9900 
INTO9910 
0109820 
ДАО О) 
INTO9940 
INTO9950 
2019960 
INI1062970 
INT09980 
ТШ ӘУ) 


WIO 
Ша 
120 
qo 
pU 


140 


DATA SET KCBCCOMPBL AT LEVEL 010 AS OF 08/24/84 


SUBROUTINE COMPBL(CASEID,XP,YP,V1T,S,DLSP ) Bits, DBRP, NBs lTRI , XC RR 


+ 


COMMON 
COMMON 


COMMON /EDDY1/ 


ЕЕ 
COMMON 
COMMON 

4 
COMMON 
COMMON 


RN ,NT, ISWPT) 


/BLCO/ 
/BLC7/ 


NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
C(100,100),D(100),DB( 100) ,DBP( 100) ,UEO( 100) ,GI 
RL,RX,SQRX,RXNTR,GMTR , GHTRS( 100) 

, ALFAS( 100) ,FFS( 100) , RTS( 100) , IEDY,NXSPT 
X( 101) ,UE( 100) ,P1( 100) , P2( 100) ,CEL,CELH 
TITLE(20),XC( 100) , YCC 100) ; 188 Q490 ) DEM: COSE 

XCTR,XTR, ISTRP,ICYCLE,ICYTL,XCTRS(2), TRFIND(2) 

/BLOW/ VN(100) 

/ISURF/ ISF 


Tem 
/ BLIN/ 


СОММОМ/РЪОТАМУРС 29 \ХУР ОС 


DIMENSION XP(100) ,DLSP( 100), YP(100) , VT( 100) ,S(100), 

+ DBPP( 100) ,DLS( 100) , CASEID( 20) 

DIMENSION XIN(100,2),YIN(100,2),S1(100,2) , VIN(100,2), 
т DIN(100,2),DELSTR( 100,2) ,DD( 100,2) ,DDD( 100,2) 
DIMENSION XB( 101) ,D1( 101) ,D2(101) ,D3( 101) , IEND(2) 

DIMENSION NBL(2),ITRI(2), XCTRI(2) 


LOGICAL TRFIND 
CHARACTER * 4 SURF(4),STITLE(2) , SURFID(4) 


DATA 
DATA 


FORMAT 
FORMAT 


FORMAT 
+ 
FORMAT 
FORMAT 
FORMAT 
+ 
FORMAT 
+ 
-- 
FORSAT 
T 
FCRMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


SURF у LE NEU RFAC OO EE. 
STITLE / 'UPPE','LOWE' / 


( 1H1,5X,' COMPUTING BOUNDARY LAYER USING HILBERT', 
" INTEGRAL ' ,/) 

( 1H0,6X,'I',9X, XP ;13X. YP .13X Sd 0E 
'DBP' 

( 180,6Х, 1',4Х, II ЗХ ТІКЕ ОТЫ. ОИК” 
13Х, 5І 125, VIN], 12 O ЕН 

( 1H ,5Х,13,5Е15.6 ) 

СІН ,ЗХЗОХОТЕ ФО 6 ) 

( 1H0,5X,'STAGNATION POINT IS FOUND BETWEEN POINT NO. ', 
I3,' AND POINT NO. ',I3 / ) 

( 1HO,5X,'INTERPOLATED STAGNATION POINT VALUES' / 
1H0O,5X,'S = ',E13.6,2X, XP = ',E13.602X. VENE 
E13.6,2X, DBP = “ils. 6, О 

( 1H0,5X,' TOTAL NUMBER OF UPPER SURFACE POINTS - ',I5, 

2X,' AND AT LOWER SURFACE = ',I5 / ) 
1H0,5X,' UPPER SURFACE DATA' ) 
1H0,5X,'LOWER SURFACE DATA' 
1H0,5X,' UPPER SURFACE CALCULATIONS' / ) 
1H0,5X,'LOWER SURFACE CALCULATIONS' / ) 
1H0,5X,'RESULTS OF POINT REDISTRIBUTION' / ) 
1H0,5X,' TABLE OF DELTA-STARS' / ( 1H ,5X,8E15.6 ) ) 
1HO,5X, TABLE OF BLOWING-VEL.' / (1H ,5X,8E15.6) ) 


ии ии“ 


1HO,5X,'NO CHANGE OF SIGN ON VT. CANNOT FIND STAG. PT.' )INT1052 


INT1000: 
INT100@ 
INT1002: 
INT100@ 
ІМТ10044 
ІМТІ00% 
ІМТІ006 
INT100A 
INT10081 
INT10091 
1ЧТ1016 
ІМТІ0СІШ 
INT10121 
INT1013! 
INT10141 
INT10159 
1ЧТ1016 
ІМТ1017 
ІХТ1018 
INT1019 
INT1020 
INT102M 
ІМТ1022 
ІМТ1023 
INT1024 
INT1025 
INT1026¢ 
INT1022 
INT1028( 
ІМТ1029 
ІМТ1030 
1ЧТ103 1 
1ЧТ1034 
INT1033 
INT1034 
INT1035 
INT1036 
ІМТ1037 
INT1038 
INT1039 
INT1040 
INT1041 
INT1042 
INT1043 
INT104 
INT1045 
INT1046 
INT1047 
INT1048 
INT1049 
INT1050 
INT1051 
























INT1053 
INT1054 
INT1053 


ACE 


230 


e 
e 


шы ыы С С O. ç 
O 
o 


600 


D) Cy C) C) G) 


СУС С С? 


DUE SERIPINPOT DATA FROM UNIT NO. 1. THE ORDER IS 
PA LOWER SURPACE T.E. TO THE UPPER SURFACE T.E. 


DO 230 I = 1,20 
mh I) = CASEIDN( 1) 
CONTINUE 

RL = RN 


DO 300 I = 1,NT 
DBP(I) = DBPP¿I) 
CONTINUE 


ERINI THE TNPZ DATA. 


ЕШІ17 - 10 

WRITE ( 6,90 ) 

WRITE ( 6,110 ) 

DO 500 I = 1,NT 

WRITE ( 6,120 ) I,XP(I),YP(1I),S(1),VT(1I),DBP(TI) 
CONTINUE 


FIND STAGNATION POINT 
DO 600 I = 2,NT 


ESROD = VT(I) * VT(I-1) 
ше УРКОШ СТ. 0. ) СО TO 600 


IS = I 
IS] = 13-1 
ЕВ ТО 700 
CONTINUE 

WRITE ( 6,210 ) 
STOP 1 


INTERPOLATE S AT VT = 0. 


6,130 ) 1S1,IS 


DS E ` (si) 

DV - ҮТ(15) - УТ(151) 

55 SEE) ES nS Any ©) 
РВВ = DBP(IS) - DBP(IS1) 

DX = XP(IS) - XP(IS1) 

DY = YP(IS) - YP(IS1) 

ШІ) T s - SGTS) 

DBS = DBP(IS) + DS1 * ( DBB/DS ) 
XPS — = XP(IS) + DS1 * ( DX /DS ) 
YPS = YP(IS) + DS1 * ( DY /DS ) 


ДИ 6,140") SS XPS, YPS, DBS 


ЕТПЕ TOTAL UPPER SURFACE POINTS. + STAG: PT. 
Ti, iS THE TOTAL LOWER SURFACE POINTS + STAG. PT. 


IU SE 
IL IS 
ПЕШТЕ ( 6,150 ) IU,IL 
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INT10560 
INT10570 
INT10580 
1110590 
INT10600 
INT10610 
INT10620 
INT10630 
INT10640 
ИД Ооо О 
INT10660 
INT10670 
INT10680 
INT10690 
INT10700 
NAO TO 
JI 10/20 
0730 
INT10740 
МИ 0 750 
INT10760 
J i 10770 
INT10780 
Ji [10790 
INT10800 
TNT OSTO 
INT10820 
MAT 10630 
INT10840 
INT10850 
INT10860 
INT108/0 
INT10880 
INT10890 
INT10900 
INT10910 
INT10920 
10550 
INT10940 
INT10950 
INT10960 
INT10970 
INT10980 
INT10990 
INT11000 
Ao To 
DISSE 
INT11030 
ІМТ11040 
INT11050 
INT11060 
JST 11070 
INT11080 
INT11090 
INT11100 
EN pl rO 


DES 


оо аз Мева осв 


900 


1000 


1100 
С 
1200 


ес) 


GROUP THE DATA FOR EACH SURFACE. FIRST, UPPER. 


DO 1200 Ly — MZ 
СО ТО (SUOMI 


L = 1 IS UPPER SURFACE 

L = 2 IS LOWER SURFACE 

M1 = IS 

M2 = NT 

IBND(L) = IU 

GO TO 1000 

M1 = 1 

M2 = IL-1 

IEND(L) = IL 

I=] 

XIN(I,L) = XPS 

YIN(I,L) = YPS 

(11) ж 

DIN(I,L) = JBS 

VIN(I,L) = O. 

Je (C JP SCE. 1 РЕШЕ 

IF ( L.EQ. 1 ) WRITE ( 6,160 ) 

IF ( L BO. 2 ) WRITE CIO” 

WRITE ( 6,112 ) 

WRITE ( 6,122 ) I,1,1,XINC1 B95 YIN(1,L),SI(1,L) ,VONC 1), 
DIN(1,L) 

END IF 

DO 1100 II = M1,M2 

I =I + 1 

IK = II 

Be ( LEO. 2) IK = E 

XIN(I,L) = XP(IK) 

YIN(I,L) = YP(IK) 

БТ СТЕ) Е ОЕ 

IF ( L .EQ. 2 ) SI(I,L) = SS=S(IK) 

VIN(I,L) = ABS(VT(IK)) 

DIN(I,L) = DBP(IK) 

IF(IP .GE. 1)WRITE ( 6,122 ) 1,11,IK,XINCI,L),YIN(I,L),SICI,L), 
VINCI, L),DINCI,L) 

CONTINUE 

CONTINUE 


RE-DISTRIBUTE POINTS ON EACH SURFACE TO A DENSER NUMBER. 


WRITE ( 6,90 ) 
DO 2000 ISF = 1,2 


NN = IENDCISF) 

ШЕК = ITRICISF) 

NXT = NBLCISF) 

XCTR = XCTRICISF) 
SURF (1) = STITLECISF) 
ICC = 1 


DO 1220 J=1,4 


84 


INT1112€ 
1NT1113(0 
INT1114( 
INT1115( 
INT1116€ 
INT1117( 
INT1118( 
INT1119€ 
ІМТ1120( 
ІМТ11214 
ІМТ1122 
ІМТ1123 
ІМТ11244 
ІМТ1125 
ІМТ1126 
ІМТ1127( 
ІМТ1128( 
INT1129( 
D 
ІХТ1131( 
га 
1NT1133( 
ІМТ1134 
ІМТ1135 
INT1136 
ІКТ1137( 
ІМТ1138( 
mel 
INT1140 
ІМТ11414 
INT11426 
INT1143€ 
INT1144( 
INT1145¢ 
INT1146( 
INT1147 
INT1148 
INT1149 
ІМТ1150 
ІМТ1151/ 
ІМТ1152( 
INT1153¢ 
INT1154( 
INT1155( 
ІКТ11564( 
INT1157€ 
ІМТ1158( 
ІМТ1159( 
ІМТ1160 
ІМТ1161 
ІМТ1162 
ІКТ1163 
ІМТ1164 
INT1165! 
INT1166 
INT1167 


SURFID(J) = SURF(J) 


1220 CONTINUE 
ДЕ (ШІЗЕ . BO. 
Е 19871267 


co es c)" 


1 ) WRITE ( 6,180 ) 
2 ) WRITE ( 6,182 ) 


SCALE = SI(NN,ISF) 


АРТ NXT,XB,D1,D2 D 


DO 1300 I = 1,NXT 
moo X (I) - XB(I) * SCALE 


C BEXNSRPODLATE S,VT,X,Y,D INTO THE NEW DISTRIBUTION. 


CALL SMFIT ( 
CALL SMFIT ( 
CALL DIFF3 ( 
CALL INTRP3( 
(ALE DIFF3 ( 
CALL INTRP3( 
ШЕШІ DIFFS ( 
CALL INTRP3( 
CALL DIFF3 ( 
MALL INTRP3( 
СТР. СЕ. 


ШІ о т тек) 21,2 7) 
IONNOSTCT, ISP). DINCI,TSE),D1,2 ) 
NN,SICI,ISF),VINC1,ISF),D1,D2,D3,0 ) 
NN,SICI,ISF),VIN(1, ISF) ,D1,D2,D3,NXT, X,UE ) 
NN,SI(1,ISF),DIN(1,ISF),D1,D2,D3,0 ) 
NN,SIC1,ISF),DIN(1, ISF) ,D1, D2, D3, NXT, X,DB ) 
КЕГСТІІ ТБ) ХІК(І ІСЕЗ 0102503,0 ) 
ANSeS PENE. TSE DD DS NET Х,ХС) 
NN,SICI,ISF),YIN(C1,ISF),D1,D2,D3,0 ) 
NN,SICI,ISF),YIN(1,ISF),D1,D2,D3,NXT,X,YC ) 
1 ) THEN | 


WRITE ( 6,190 ) 

WRITE ( 6,110 ) 

ШІ 1520 І- 1.NXT 

ШЕШЕ (G 120 ) I.XC(I),YC(I),X(I),UE(I),DB(I) 


1320 CONTINUE 


2000 CONTINUE 


rane 


END IF 
C 
C INPUT TO THE B.L. PROGRAM X,UE,DB,XC,YC ARE NOW DEFINED. 
С 
DO 1350 I = 1,NXT 
DBP(I) = DB(I) 
D(I) = DB(I) 
1350 CONTINUE 
С 
CALL BL2D( ITR,ISWPT,SURFID) 
ë 
CALL DIFF3 ( NXT,X,DELS,D1,D2,D3,0 ) 
CALL INTRP3( NXT,X,DELS,D1,D2,D3,NN,SI(1,ISF),DELSTR(1,ISF) ) 
GALL DIFF3(NXT,X,D,D1,D2,D3,0) 
CALL INTRP3(NXT,X,D,D1,D2,D3,NN,SI(1,ISF),DD(1,18F)) 
CALL DIFF3(NN,SI(1,1SF),DD(1,ISF),DDD(1,1SF),D2,D3,0) 
TRFIND(ISF) = . FALSE. 
IF(ITR .EQ. 3 .AND. МТВ .LE. NXT) THEN 
XCTRS(ISF) = XCTR 
TRFIND(ISF) = . TRUE. 
END IF 
б 


PUT THE TWO SURFACES DELTA-STARS BACK TO ONE STRIP 


DELSTR(1,2) = 0. 5*(DELSTR(2,1)+DELSTR(2,2)) 


INT11680 
INT11690 
700 
INT11710 
ТАГЫ 20 
INT117/30 
INT11740 
INT MSO 
INT11760 
INT11770 
INT11780 
ТНТ 220 
INT11800 
INT11810 
INT11820 
INT11830 
INT11840 
INT11850 
ТИЕ 50 
INT11870 
INT11880 
INT11890 
INT11900 
INI11910 
Пуа ите 20 
INT11930 
INT11940 
[ST 115 09 
INT11960 
INT11970 
INT11980 
Tq 99D 
INT12000 
ШУБЕ ОО О 
12020 
INT12030 
INT12040 
INT12050 
INT12060 
INT12070 
INT12080 
INT12090 
TS T 2 100 
182 110 
INT12120 
INT12130 
INT12140 
IN112 130 
INI 12160 
Де 
INT12180 
ТМ 190 
ТЫЛ 2 020 
ТОТ 10 
INg 220 
ЛБ 22500 


2100 


2200 


CONCE 


рер Фр 


180 


260 


8PLSTR(1,1) = php m s 


DD(1,2) = 0.5 * (DD(2,1) + DD(2,2)) 

DD(1,1) = DD(1,2) 

DDD(1,2) = 0.5 * (DDD(2,1) + DDD(2,2)) /SORTCRL) 
DDD( 1.1) ~ DILOR) 

IL = IEND(2) 

I = Th 

ПО ТОСО атии 

I I-1 


DLS(I) š DELSTR(II,2) 
VN(I) = DDD(II,2)/SQRT(RL) 


DBPP(1) = DD(11,2) 
CONTINUE 

I = 1551 

IU = IEND(1) 

DO 2200 11=2, 1U 

I =. T jl 

DLS(I) = DELSTR IID 


VNCI) = DDD(II,1)/SQRT(RL) 

DBPP(1) - DD(11,1) 

CONTINUE 

WRITE ( 6,200 ) ( DLS(I), I=1,NT ) 
WRITE ( 6,220) (VN(I) , I=1,NT) 


RETURN ` 

END 
DATA SET KCBCCOMPGI AT LEVEL 001 AS OF 08/24/84 
DATA SET KCECCOMPGI AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCCOMPGI AT LEVEL 003 AS OF 08/24/84 

SUBROUTINE COMPGI 


CALCULATE GI 


COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 


COMMON /BLC7/ C(100,100),D(100),DB( 100) ,DBP( 100) ,UEO( 100) ,GI 
SUM2 - 0.0 

N1 =N% =+] 

N2 = NXT 

DO 180 K = N1,N2 

SUM2 = SUM2 + C(K,NX)* (D(K) - DBP(K)) 

CONTINUE 

N1 = 2 

N2 = NX-1 

SUM1 = 0, 

DO 260 K = N1,N2 

SUM1 = SUM1 + C(K,NX)* (D(K) - DBP(K)) 

CONTINUE 

GI = UEO(NX) + SUM1 + SUM2 - C(NX,NX) * DBP(NX) 
RETURN 
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INT1224( 
INT1223 
INT12261 
INT1227 
INT1228 
INT1229 
INT1230 
INT1231( 
INT1232¢ 
INT12330 
INT12346 
INT12350 
ІМТ12360 














INT12460 
INT12470 
INT1248( 


INT12530) 
INT1254@ 
INT1255@ 
INT1256§ 
INT1257@ 
INT )I2 569 
INT1259@ 
ІМТ12600 
ІМТ1261С 
INT1262¢ 
ІМТ126360 
ІМТ12640 
ІКТ12650 
INT12660 
INT1267@ 
INT12680 
INT12690 
INT1270@ 
INT12718 
INTI1272g 
INT12730 
INT12740 
INT1275@ 
INT12760 
INT1277@ 
INT12780 


ada C2 


С2 07 €) 


CCC CC C) 


10 


Ra 


E C2 €0 C2 


END 
DATA SET KCBCDIFF3 AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCDIFF3 AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCDIFF3 AT LEVEL 002 AS OF 04/05/84 
SUBROUTINE DIFF3 (N,X,F,FP,FPP,FPPP,IEND) 


ШШНЕСШТЕМЕС TRESDERIVATIVE OF THE INPUT FUNCTION AT THE INPUT PTS. 


ПИШЕ КЫЛОО Х(101).ЕСІ01),ЕР(101), ЕРЕСІОІ?),ЕРРР(101) 


FIRST DERIVATIVES USING WEIGHTED ANGLES 


N1=N-1 
DX=X(2)-X(1) 
DF=F(2)-F(1) 

Awe2= ATAN2(DF,DX) 
DL2=DX 


PO 10 I=2,N1 

ANG1=ANG2 

DL1=DL2 

]i=I+1 

DN=X(11)=X(I) 

DF=F(11)-F(1) 

ANG2= ATAN2(DF,DX) 

DL2=DX 
ANG=(DL2*ANG1+DL1*ANG2)/(DL1+DL2) 
FP(I)= TANCANG) 


ШТІ. ХЕ. 2) СО Т0770 


ANGI = ANG 
ANG11 = ANGI 
DLI = DL1 
CONTINUE 
ANGF = ANG 
ANG2F = ANG2 
DLF = DL2 


IEND1 = IEND + 1 
ЕТО СЗТ) TENDI 


IEND = 0 , EXTRAPOLATE FOR END VALUES 
шиш) - 2 s (b(2)-F(1 ))/DLI - FP(2) 
БЕ) = 2.*(FiN)=PCNT))/0LF - FP(N1) 
GOSTOTO ` 


IEND = 1, DERIVATIVES ARE CONTINUOUS ACROSS ENDS 


ANG = (DLI*ANG2F + DLF*ANG1I) / (DLI + DLF) 
FP(1) = TANCANG) 
FP(N) = FP(1) 
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Id 7.20 
INT12800 
1 1123810 
INT12820 
INT12830 
INT12840 
INT12850 
Ти 12365 
ПОДОЛ 
1112650 
II O 
T1228 06 
I 2o 20 
| 112920 
NIS 
INT12940 
1112950 
ІМТ12960 
TS TD 
О 
о 
INT13000 
PNASESQTO 
ENS ZO 
INT13030 
INT13040 
050 
INT13060 
INT13070 
INT13080 
во 
PALEO TOO 
INTIS КО 
1115120 
Iun SO 
INT13140 
INT13150 
INT13160 
TIS 
INT13180 
О-о 
1 11 200 
TNES 2O 
E220 
INT 13230 
INT13240 
INTI 250 
INT13260 
INTI O 
ING 230 
1-20 
INT13300 
АДА TO 
1113320 
О 550 


GO TOT 0 
С IEND = 2, IF ИТ ТЕ 


Ш CONTINUE 


FCD = 2 КС В, 
IF (FP (1) .LT. 0.0) FP (1) = 0.0 
ЕР(М) = ОСЕМ) ЕОМ ОТЕ ЗЕР 
C 
C SECOND + THIRD DERIVATIVES USING CUBIC FITS 
C 
20 DO 30 I=2,N1 
11 =1-1 
12 =1+1 
0х1 = X (11) - X (1) 
DX2 - Х (12) - Х (1) 
DF1 = 2.0 * ((F (11) - Е (1)) / DX1 - FEP (T) PR h 
DF2 = 2.0 * ((F (12) - F (1)) / DX2 - FP (1)) / DX2 
FPPP(1)= 3.0 * (DF1 - БЕ?) / (DX1 - DX2) 
FPP (1)= DF1 - DX1 * FPPP (1) / 3.0 
30 CONTINUE 
FPPP(1)= FPPP (2) 
FPPP(N)= FPPP (N1) 
FPP (1)= FPP (2 ) + (X (1) = ж (2 s K O RK R 
FPP (N)= FPP (N1) + (X (N) - X (N1)) * FPPP (N1) 
C 
RETURN 
END 
Е DATA SET KCBCEDDY АТ LEVEL 001 AS ОҒ 08/24/84 
С DATA SET KCBCEDDY AT LEVEL 001 AS OF 08/24/84 
С DATA SET KCBCEDDY АТ LEVEL 003 AS OF 04/05/84 
SUBROUTINE EDDY 
@ 
C CALCULATE EDDY VISCOSITY USING C .S. TWO-LAYERS EDDY 
C VISCOSITY FORMULA 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
COMMON /BLC1/ F(101,2),U(101,2),V(101,2),V(101,2) ,B( 101,2) 
COMMON/BLC7/ C(100,100),D(100),DB( 100) ,DBP( 100) ,UEO( 100) ,GI 
COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR,GMTRS( 100) , 
+ ALFAS(100),FFS(100),RTS(100),IEDY,NXSPT 
COMMON /GRD / ETA(101),DETA(101),A(101) 
COMMON /GTY / X(101),UE(100) ,P1(100),P2(100) ,CEL,CELH 
DIMENSION FINT(101) 
С 
С =“ = = = = = =~ we = = = = о т тм >> 
E 


J0=1 
UDEL=0. 995**U(NP, 2) 
DO 10 J=1.NP 
IF(U(J,2).LT.UDEL) JJ=J 
10 IF(UCJ,2). ERROT O) ЕЕ 


. LT. 0.0 RESPIR ЛЕНЕ 


EDEL-ETA(JJ)T(ETACJJSIO SETACJJ ) ) СОС а айы 


n *(UDEL-U(JJ,2)) 
DO 15 J=1,NP 
ETADEL=ETA(J)/EDEL 
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15 


20 


30 


рее 


E 


70 


C 80 


90 


IF(ETADEL. GT. 1. 0) ETADEL-1.0 
FINT(J)=1. 0/(1. 0+5. S*ETADEL**6) 
CALL AMEAN(1,JJ,ETA,FINT, 2) 

RU = RL 

NND) cO TO 20 
ALFAS(NX) = ALFAS(NX-1) 

FFS(NX) = FFS(NX-1) 

RTS(NX) = RTS(NX-1) 


GMTR = GMTRS(NX) 
IF (NX . LE. NS) RU = RL * UE(NX) 


RL2 = SQRT(RU * X(NX)) 

RL4 — = SQRT(RL2) 

RL216 = 0.16 * RL2 

И Б Оо) е) 
DO 30 J=2,NP 

УМ = 0.5 * (V(J,2)-V(J,1)) 
IF(VM .GT. VMAX) VMAX = VM 
CONTINUE 


ШИЛЕР „ЕО. 0) СО TO 35 
ТЕ. 1 ОК GTR ЕТ. 0 85 ОК КК .. DE. NTR+3 ) 
> 


MODIFY ALFA USING SIMPSON'S ARGUMENTS 


CALL SMPSON 
ALFA = ALFAS(NX) 


EDVO = ALFA * RL2 * GMTR * (U(NP,2)*ETA(NP) - F(NP,2)) 
DO 40 J=2,NP 

JJ =J 

YBA _ = RLA*SQRT(VMAX)/26. O*ETA(J) 

EL = 1.0 

IF(YBA .LT. 10.0) EL = 1.0 - EXP(-YBA) 

ЕРУІ = RL216 * GMTR * (EL*ETA(J))**2 * ABS(V(J,2)) 


КЕРТ от. EDVO) GO TO 70 

B(J,2) = 1.0 + EDVI*FINT(J) 

RAS BAR 
B(J,2) = 1.0 + EDVI 


CONTINUE 
JM = 2 
BM BD. 2) 


DO 50 J=2,NP 

IF(BM. GT.B(J,2)) GOTO 50 

JM = ] 

BM = о 

CONTINUE 

GOTO 90 

DO 80 J=JJ,NP 

B(J,2) = 1.0 + EDVO*FINT(J) 
B(J,2) = 1.0 +EDVO 


CONTINUE 
В 2). = 1.0 
JJ = 1 


DO 100 J=2,NP 


89 
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EXT 13900 
Pris 9 10 
Т5 920 
250 
INT13940 
1 1125550 
Ш 10960 
т. 11-970 
INT13980 
TERES 99D 
INT14000 
INT14010 
INT14020 
INT14030 
INT14040 
INT14050 
INT14060 
INT14070 
INT14080 
INT14090 
INT14100 
INT14110 
INT14120 
INT14130 
INT14140 
INT14150 
INT14160 
INT14170 
INT14180 
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INT14200 
INT14210 
INT14220 
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INT14240 
INT14250 
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INT14330 
INT14340 
INT14350 
INT14360 
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A = 
100 CONTINUE 
ТЕСЈЈ. ЕД. 1) GO BO 110 
б 
С IN THE SEPARATED REGION, EDDY VISCOSITY IS SET EQUAL TO 
C THAT IN THE PREVOUS STATION TO AVOID NUMERICAL TROUBLE 
1213 = JJ +3 
JJP3  - MIN(JJP3, NP) 
CALL AMEAN(1,JJP3,ETA,B( 1,2) ,2) 
110 CALL AMEAN(1,NP,ETA,B(1,2),1) 


RETURN 

END 
DATA SET KCBCEXTRAP AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCEXTRAP AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCEXTRAP AT LEVEL 008 AS OF 02/13/84 

SUBROUTINE EXTRAP(NX,NXTE,X,Y) 


O O O 


(22702 


EXTRAPOLATE DATA USING PARABOLIC FORMULA 
DIMENSION X(101),Y(101) 


X(NX-1) 

X(NXTE) 

Xl -X3 

NOUS 

(ҮЗҮ Л Х2%%0) 
Ү1 - ВВ % ХІЗ>2 
-NX,NXTE 

X(I) -X3 


DO 
X1 
Y(I) 

10 CONTINUE 
RETURN 
END 


Ha 
O 


Ба 
Кә 
I l -—al H H H H l H H I 


DATA SET KCBCFILLUP AT LEVEL 001 AS OF 08/24/84 

DATA SET KCBCFILLUP AT LEVEL 001 AS OF 08/24/84 

DATA SET KCBCFILLUP AT LEVEL 007 AS OF 04/05/84 
SUBROUTINE FILLUP( INDEX) 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
COMMON /BLC1/ F(101,2),U(101,2),V(101,2),W(101,2),B(101,2) 
COMMON /GRD / ETA(101),DETA(101),A(101) 


C2 C2 C2 


aaa 


TP CNRS ОВ ОВА ОВА 
АВИ О Ре ОКО 0 


a a 


DEE TNE PROEMEES FORD. ТООККО ТЫ 
NP1 МЕ 1 
NP NP #2 
NP MINCNP, NPT) 
NM NP 
GOTO 20 


90 
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ІМТ1463 
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INT1469 
INT1470 
INT1478 
INT1472 
INT1473 
INT147&4 
INT1475 
ІКТ1476 
INT1477 
INT1478 
INT1479 
INT1480 
ІМТ1481 
ІМТ1482 
INT1483 
INT1484 
INT1485 
INT1486 
INT1487 
INT1488 
INT1489 
INT1490 
INT1491 
INT1492 
INT1493 
INT1494 
INT1495 
ІМТ1496 
ІМТ1497 
INT1498! 
INT1499! 


сате) 


a daa 


ЕТІ. UP PROFILES BEFORE MOVING TO THE NEXT STATION 


+ 


NP1 = NP + 1 
NM = NPT 

DO 30 J=NP1,NM 

ВОО ЕСТЕ] 2) + DETA(J-1)*U(J-1,2) 


U(J,2) z U(J-1,2) 
V(J.2) = 0.0 
B(J,2) = B(J-1,2) 
ит. 2) 
CONTINUE 

RETURN 

END 


DATA SET KCBCHEADER AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCHEADER AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCHEADER AT LEVEL 001 AS OF 04/05/84 
SUBROUTINE HEADER ( TITLE,SURFID,ISTRP ) 
COMMON /ISURF/ ISF 


DIMENSION TITLE(20), SURFID(4) 
FORMAT ( 1H1,20X,20A4 ) | 
FORMAT ( 1H0,15X,'BOUNDARY LAYER CALCULATION FOR ', 
РРЕВ SUPPA CE оса IC YCLE= 15%) 16Х 71(ІН-У / ) 
FORMAT ( 1H0,15X,'BOUNDARY LAYER CALCULATION FOR ', 
ШЕШЕК ӘШЕБАСЕ | ТОХ. ТОӨМЕШЕ- „пори 16X,71(1H-) / ) 


ATIE ( 6,10 ) TITLE 
(ОВ ГОМЕ) WRETE ( 6,20 ) ISTRP 
Meteors. BO. 2) WRITE ( 6,30 ) ISTRP 


RETURN 
END 
DATA SET KCBCINPUT AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCINPUT AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCINPUT AT LEVEL 009 AS OF 08/24/84 
SUBROUTINE INPUTCITR,ISWPT,SURFID) 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT, INVRS,NS,IP 
COMMON /BLC1/ F(101,2),U(101,2),V(101,2) ,W( 101,2) ,B( 101,2) 
COMMON /BLC2/ DELF(101),DELU(101) ,DELV(101) ,DELW(101) | 
COMMON /BLC7/ C(100,100),D(100),DB(100) ,DBP( 100) ,UEO( 100) ,GI 
COMMON /BONV/ ITMAX,EPSL,EPST,CONV 
COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR,GMTRS( 100), 
ALFAS(100),FFS(100),RTS(100),IEDY,NXSPT 
X(101),UE(100) ,P1( 100) ,P2( 100) , CEL, CELH 
TITLE(20),XC(100),YC( 100) , I$G( 100) , DELS( 100), 


COMMON /GTY / 
COMMON /BLIN/ 


T XCTR,XTR,ISTRP,ICYCLE,ICYTL,XCTRS(2) , TRFIND(2) 


COMMON/TRN/ PGAMTR,OMEGA, RTHETB , RTRANB 
COMMON /ISURF/ ISF 

DIMENSION D1(100),D2(100),D3(100) 
DIMENSION SURFID(4),XCS(100) 

LOGICAL TRFIND 
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INT15140 
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I 110] O 
LAARS FOO: 
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Т 5200 
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Jol l 220 
INT15230 
INT15240 
101520 
INT15260 
NO 2 О 
DS 280 
ИЕ 290 
INDIS 300 
ЈАКО 
Ti: 15320 
О 
INT15340 
INT15350 
INT15360 
ВО 
ИБ зоо 
В О 
INT15400 
INT15410 
INT15420 
INT15430 
INT15440 
INT15450 
INT15460 
INT15470 
INT15480 
INT15490 
INTIS 300 
ЛИ БИО 
1 15520 
| 115550 
INT15540 
SAS 550 


(ЕМЕС) 


G SEARCH FOR 


65 


70 


80 
C 
C 


C 
C 


25 


100 


ITHAX 
EPSL 
EPST 
NPT 
ETAE 
VGP 
DETA1 
NSS 


NXT / 4 


PRESSURE PEAK AS SWITCH POINT 


ОМАХ = ЏЕ(1) 


DO 50 I = 2 , NXT 

IF (UE(I) .LE. UMAX) GO TO 55 
UMAX = UE(I) 

CONTINUE 

GO TO 60 

NS = I - á 

IF (NS .GT. NSS) NS = NSS 

IF (NS .LT. 3) NS = 93 


CALCULATE THE PRESSURE PARAMETERS Pl + P2 FOR В. L. CALCULATION 
USING TRANSFORMED COORDINATES 


CALL DIFFS (ea, X, UE. ӘКІ ве 0 


DO 65 I = 
P2(1) 
P1(1) 
CONTINUE 
P1(1) 
XCMIN 
MIN 


| | 


DO 70 I=1, 


2 Mar 
X(T) = DICIT ONES 
A CI) 


0.5 * (1.0 + P2(1)) 
XC(1) 

1 

NXT 


ТЕСКЕМЕ ХЕ БЕНО УО 


XCMIN 
MIN 
CONTINUE 
DO 80 I = 


| | 


LE (ieee 


XCS(I) 
ISG(I) 
ELSE 
ХС5(1) 
ISG(I) 
END IF 
CONTINUE 


INVRS = 


ІТКРІ = 


GOTO (150, 


XC(I) 
iD 


pa 


, NXT 
MIN) THEN 

ХС(І) 

] 


| Il 


SMS 
-1 


| | 


NS 


SEARCH FOR TRANSITION LOCATION 


EN] 
95, 120, 150), ITRP1 


TRANSITON LOCATION IS INPUT ( = XCTR) 


ро ОО ТЕ 
TEC XCTR. LY. X€SC 1) ) боле 


CONTINUE 
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INT1571g 
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INT1573( 


INT1574( 
ІМТ157% 
ТКТ 1576 
INT15 78 
INT1578] 
INT1579] 
INT1580t( 
INT1581g 
INT1582( 
INT1583( 
INT1584( 
INT1585g 
INT1586( 
ІМТ1587( 
INT1588( 
INT1589( 
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INT1595@ 
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17 Т1598ЦЕ 
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INT16000 | 
INT16010 
INT16028 
INT16030 | 
INT1604 

INT16950% 
ТКТ16060 | 
ІМТ16070| 
INT16080 | 
INT16090 | 
INT16100] 
ІМТ16110 


ее нын ны ы = ee ee Oe d 


NTR = NAT + 1 
GOTO 200 

105 NTR = 1-1 
SR. ET. 3) THEN 


NTR = 3 
XCTR = XC(NTR) 
END 1F 
GOMO 200 
C 
Ш) TRANSITION LOCATION IS SET AT THE PRESSURE PEAK 
C 
ШО UM IE T) 
Ши = 1 
Doms) I = 1,NXT 
ОН. СТ, UE(I)) GOTO 75 
IM = 1 
UM = UEC IM) 


75 CONTINUE 
пати LT. 4) IM = 4 


NTK = IM 
XCTR = XC(NTR) 
GOTO 200 
C 
C | TRANSTION LOCATION WILL BE CALCULATED BASED ON MICHEL CRITERION 
C 
150 NIR = NXT + 1 


200 DO 210 I-1,NXT 

#10 GATRS(1)= 0.0 

C 

C TRANSITION LOCATTON PROVISIONALLY FROM PREVIOUS CYCLE 
C 


IF (TRFIND(ISF) ) THEN 

DO 211 1 =1 , NXT 

ОЗОТ) = XC7T) 

IF (1.11. MIN) XCS(I) - 
211 CONTINUE 

DO 215 I=1,NXT 

ШЕ ХОТЕ5( ТСЕ. ТЕ ESA) GOTO 217 
215 CONTINUE 
217 NTR = I-1 

XCTR = XCTRS(ISF) 

IF (NTR .LT. 3) THEN 


SACOG.) 


NTR = 3 
XCTR = XC(NTR) 
END IF 
END ТЕ 
C 
C CALCULATE GAMTR DISTRIBUTION 
C 


ТВТ. МР GOTO) 256 
FAC = (XCTR-XC(NTR))/(XC(NTR+1)-XC(NTR)) 


XIR = X(NTR) + FAC*(X(NTR+1)-X(NTR)) 
UETR = UE(NTR) + FAC*(UE(NTR+1)-UE(NTR) ) 
RXNTR = XTR * UETR * RL 

GGET = 1. 0/PGAMTRERL *2/ВХМТВ 1. 34 


DO 220 I=NTR+1,NXT 
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INT16400 
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INT16560 
L 116570 
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NOS 
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ALFAS(I) = 


220 GNTRS(I)= 1.0 
ALFAS(NTR) = 0.0168 
UEINTG = 0.0 
Ul = 0.5/UETR 
X1 = XTR 
DO 230 I=NTR+1,NXT 
U2 = 0.5/UE(I) 
X2 = х(1) 
UEINTG = UEINTG+(U1+U2)*(X2-X1) 
01 = 02 
X1 = X2 
GG 2 GGFT*UEINTG*(X(I)-XTR)*UE(I)s**3 
IF(GG .GT. 100m0 REBOEO 
GMTRS(I) = 1. 0-EXP(-GG) 
230 CONTINUE 
250 CONTINUE 
C 
C GENERATE B. L. ETA GRIDS + INTIAL VELOCITY PROFILES 
C 
CALL INTL(ETAE ,DETA1,VGP) 
DO 260 I=1,NXT 
UEO(I) = UE(I) 
260 CONTINUE 
С PRINT OUT INPUT DATA 
C 
IF (ICYCLE .GE. ICYTL-1 .OR. IP .GE. 0) THEN 
CALL HEADER( TITLE,SURFID,ISTRP ) 
WRITE(6,1002) NXT,ITR,IP,NS,NTR,ISWPT 
WRITE(6,1003) VGP,DETA1,RL,XCTR, OMEGA , PGAMTR 
ELSE 
IF (ISF.EQ.1) WRITE(6,1008) ICYCLE 
IF (ISF.EQ.2) WRITE(6,1009) ICYCLE 
END IF 
IF (NTR. LT. NXT) THEN 
IF (ITR.EQ.1) WRITE (6,1005) XCTR,XTR,NTR 
IF (ITR.EQ.2) WRITE (6,1006) XCTR,XTR,NTR 
IF (TRFIND(ISF)) WRITE(6,1007) XCTR,XTR,NTR 
END IF 
RETURN 
C 
@ -- - - o - - - m m m m mmm m въ еее es uU 0 
Е 
2 РОВМАТ( 2044) 
3 FORMAT( 1015) 
4 FORMAT(6F10. 0) 
1001 FORMAT(1H1,20X,20A4) 
1002 FORMAT(1HO,10H МХТ = ,15,7Х,10Н ІТК = ,15,7Х/ 
2 1H ,10H ПР = ,15,7X,10H NS  - ,I5,7X/ 
+ 1H ,10H NTR = ,I5,7X,10H  ISWPT= ,15) 
1003 FORMAT(1HO,10H VGP = ,E12.4,10H DETAI= ,E12.4/ 
P 1H ,10H RL = ,E12.4,10H  ХСТВ = ,E12.4/ 
+ 1H ,10H OMEGA = ,E12.4,10H PGAMTR= ,E12.4) 
1004 FORMAT(1HO,3X,2H 1,6X,2HXC,11X,2HYC,11X,2H X,11X,2HUE,11X,2HP1, 


0.0168 


+ 11X,2HP2,/(1H ЗОО Ва а 


1005 FORMAT(/3X,' BEGIN OF TRANSITION IS BEING INPUT AT X/C 


94 


. INT1713% 


=",F8.4,4X, INT17230 
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INT17000. 
INT17010 
INT17020 
INT17030 
INT17040 
INT17050 
INT17060 
INT17070: 
INT17080 
INT17090 
INT17100 
INT17110 
INT17120 


INT17140 
INT17150 
INT17 160 
INT17170 
INT17180 
INT1719 ) 
INT17200 
INT1721( 
INT17220 


1007 


1008 
1009 


C 
C 
C 


СС С СУС? 


10 
20 


30 


Сз С) + 
=> 


80 


90 


120 


+ ED US AX Z' 135 Оде 
1006 FORMAT(/3X,'BEGIN OF TRANSITION IS SET AT PRESSURE PEAK, X/C =', INT17250 
-- В S/C = p8.4.5X  NTR =',13/) INT17260 
FORMAT(/3X, ‘BEGIN OF TRANSITION IS PROVISIONALLY TAKEN FROM ', I7 270 
PREVIOUS CYCLE: X/C=',F8.4,4%, S/C = ,F8.4,4X%, NIR =' ,13/) INT17280 
FORMAT(//3X,'UPPER SURFACE CALCULATIONS OF CYCLE',I3) INT17290 
FORMAT(//3X,'LOWER SURFACE CALCULATIONS OF CYCLE',I3) INT17300 
END INT17310 
DATA SET KCBCINTL AT LEVEL 001 AS OF 08/24/84 INT17320 

ОБ СЕЕ ИИ AT T EVEL 001 AS"OF 08/24/84 INT17330 

DATA SET KCBCINTL AT LEVEL 009 AS OF 02/22/84 INT17340 
SUBROUTINE INTL(ETAE,DETA1,VGP) INT17350 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP INT17360 
Е ВСІ FC101,2),UC101,2) ,VC101, 2) W101, 2) ,BC 101, 2) ENS 0 
COMMON /BLC2/ DELF(101),DELU( 101) ,DELV( 101) , DELW( 101) INT17380 
ШЕШКІСІ /БІС6/ 51(101) ,52(101),63(101),58(101) 65(101),56(101), INT17390 
+ Sy 101 AI) AO Gone БЗ( 101) ,К4( 101) ІКТ17400 
COMMON /BONV/ ITMAX,EPSL,EPST,CONV INT17410 
COMMON /GRD / ETA(101),DETA(101),A(101) INT17420 
COMMON /GTY / X(101),UE(100),P1(100),P2( 100),CEL,CELH INT17430 
INT17440 

--------------------------------------.--.-.-..-..................- ІМТ17450 
INT17460 

GENERATE THE GRID INT17470 
INT17480 

DETA(1) = DETAIL INT17490 
IF(VGP. LT. 1.0) VGP = 1.0 пи QD 
IF((VGP-1.0) . LE. 0.001) GO TO 10 INT17510 
NP = ALOG((ETAE/DETA(1))*(VGP-1.0)+1.0)/ALOG(VGP) + 1.001 INT17520 
GO TO 20 INT17530 
NP = ETAE/DETA(1) + 1.001 INT17540 
IF(NP .LE. NPT) GO TO 30 7 550 
WRITE(6, 150 ) INT17560 
STOP INT17570 
A OO INT17580 
ВАА АО J=24NPT INT17590 
ЕТА(Ј) = ЕТА(Ј-1) + РЕТА(Ј-1) INT17600 
DETA(J)2 VGP*DETA(J-1) INT17610 
А(Ј) 2 0. 5*DETA(J-1) INT17620 
CONTINUE INT17630 
INT17640 

GENERATE INITIAL VELOCITY PROFILE INT17650 
DO 90 J=1,NP INT17660 
ЕТАВ = ЕТА(Ј) /ЕТА(М№Р) INT17670 
ETAB2 = ETAB**2 INT17680 
F(J,2) = 0. 25*ETACNP)*ETAB2*(3.0 - 0. 5*ETAB2) INT17690 
UCJ,2) = 0.5*ETAB*(3.0 - ETAB2) INT17700 
VCJ.2) - 1.5*(1.0 - ETAB2)/ETA(NP) INT17710 
Bowe) = 1.0 AO 
W(J,2) = 1.0 INT17730 
CONTINUE INT17740 
NX = 1 ING io? 50 
Т = 0 INT17760 
ІТ = ІТ +1 INT17770 
DANA LT НАХ) GO TO 130 INT17780 


s 
сл 


130 


150 


CA CU C) 


ОШОЛ СУ ООСО УОС СУС С С са СО КО 


STOP 


CONTINUE 

DO 140 J= 2,NP 

FB ПІЗ То шр 

UB - 0,55 ШІ, +H ssp 

VB = 045 GG. 2) ада леене, 

USB = ОБОО Е T oe 

DERBV = (B(J,2)*V(J,2) - B(J-1,2)*V(J-1,2))/DETA(J-1) 
FVB = 0.5%(F( МСт,2) + Е(1-1,2) О ЩЩ, 
S1(J) = 0. 5*PI(NX)*F(J,2) + ВКЛ. 2) реша В 
S2(J) = 0.5*P1(NX)*F(J-1,2) - B(J-1,2)/DETA(J-1) 
S3(J) = 0. 5*PIC(NX)*V( J, 2) 

S4(J) = 0. 5*PI(NX)*V(J-1,2) 

S5(J) = Вит, 

S6(J) = -P2(NX)*U(J-1,2) 

CRB = -P2(NX) 

R2(J) = CRB - (DERBV + P1(NX)*FVB - P2(NX)*USB) 


R1(J) = HUARAZ) то DE G тр 
R3(J-1)=W G.J 2) - ОСОЛ) ка ет p Е 
CONTINUE 

Ке) 

RAC) 

ROCNE) 
CALL SOLV3 

IF(ABS(CDEDVOID)) GT. БЕБІ ) 6009125 
CALL FILLUP( 2) 

CALL OUTPUT(1) 


WoW Il 
сос > 
ого 


RETURN 


FORMAT(1HO,37HNP EXCEEDED NPT -- PROGRAM TERMINATED) 
END 
DATA SET KCBCINTRP3 AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCINTRP3 AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCINTRP3 AT LEVEL 003 AS OF 04/05/84 
SUBROUTINE INTRP3 (N1,X1,F1,FP1,FPP1,FPPP1,N2,X2,F2) 


CUBIC INTERPOLATION 


GIVEN THE VALUES OF A FUNCTION (F1) AND ITS DERIVATIVES 
AT N1 VALUES OF THE INDEPENDENT VARIABLE (X1) 


FIND THE VALUES OF THE FUNCTION (F2) 
AT N2 VALUES OF THE INDEPENDENT VARIABLE (X2) 


X2 CAN BE IN ARBITRARY ORDER 


DIMENSION X1( 101) ,F1(101),FPiC 101) ,FPP1( 101) 5 FPPPiC Om) | 22a 
+ ECT) 


DATA EPS /1E-04/ 
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1ЧТ1817! 
INT1818f 
INT1819f 
ІМТ1820( 
ІМТ182 18 
INT1822( 
INT1823( 
INT1824( 
ІМТ1825( 
ІМТ1826 
ІХМТ1827 
ІМТ1828 
ІМТ1829( 
INT1830( 
ІМТ1831 
INT1832( 
ІМТ1833 
1NT1834( 





20 


30 


Creo) 


б р ш фр Ше» де ср 


Crt CO) 


10 


5 


ІМПЕХ 
ІМПЕХ 


JT = 2 
DO 10 I = 1,N2 

ХМ = X2(I) 

EEO J= JT,N1 

Ji =J -1 

ЖЕ (107) СЕ УМ ) со TO 30 
CONTINUE 
J = N1 
JT = J 
DXX 

DXX2 
DXX3 
F2(1) 


NOCT) SEIT) 
DXX * DXX / 2. 
в, о 
F1(J1) + DXX*FP1(J1) + DXX2*FPP1(J1) + DXX3"FPPP1(J1) 


RETURN 
END 
DATA SET KCBCJOIN 
DATA SET KCBCJOIN 
DATA SET KCBCJOIN 
SUBROUTINE JOIN( INDEX) 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
COMMON /BLC1/ F(101,2),U(101,2),V(101,2) ,W( 101,2) ,B( 101,2) 
COMMON /BLC7/ C(100,100),D(100),DB(100) ,DBP( 100) ,UEO( 100) ,GI 
COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR,GMTRS( 100) 


AT LEVEL 001 AS OF 08/24/84 
AT LEVEL 001 AS OF 08/24/84 
AT LEVEL 012 AS OF 02/20/84 


2 .ALFAS(100),FFS(100),RTS(100), IEDY,NXSPT 


COMMON /GTY / X(101),UE(100),P1(100),P2(100),CEL,CELH 
COMMON /GRD / ETA(101),DETA(101),A(101) 

COMMON /SAVE/ FS(101),US(101),VS(101),BS(101) ,WS( 101) 
COMMON /SMRY/ VW(100),ITP(100),ISL(100),DLS(100),CF(100), 


1 THT( 100) ,NPSTR( 100) 


meron (HE FIRST SWEEP 
E OR SUBSEQUENT SWEER 


CALL COMPGI 

eel = C(NX,NX) 

Wes) = GI / (1.0 = BDLS(NX) * SORTCRL) * CII ) 
IF( INDEX. EQ. 1) GOTO 15 


RETRIEVE PROFILES AT SPATION Ne BOR INVERSE В. L. 
CALCULATION 


DO 10 J=1,NPT 
F(J,2) = FS(J) 
U(J,2) = US(J) 
МСТ = vS(J) 
W(J,2) 2 WS(J) 
B(J,2) = BS(J) 


= UES/W(1,2) 
SOS # 1.0 


505 = 1.0 / SQRTCUES) 
DO 20 J=2,NPT 


pi) 


INT18350 
INT18360 
ИСУО 
J 165580 
INT18390 
INT18400 
INT18410 
INT18420 
INT18430 
INT18440 
INT18450 
INT18460 
INT18470 
INT18480 
INT18490 
INT18500 
Diao 
INT18520 
INT18530° 
INT18540 
INT18550 
INT18560 
INT18570 
INT18580 
J Tl 70 
INT18600 
INT18610 
INT18620 
INT18630 
INT18640 
INT18650 
INT18660 
INT18670 
INT18680 
INT18690 
INT18700 
INT18710 
INT18720 
INTIS730 
INT18740 
INT18750 
INT18760 
080015770 
TM OO 
INT18790 
INT18800 
INI15S10 
INT18820 
INT18830 
INT18840 
INT18850 
1ЧТ18860 
INT18870 
INT18880 
А“ шосето 
INT18900 


20 


30 


60 


65 
70 


са са со Сао) 


ЕТА(Т) - ЕТА( "9508 


DETAJIMIN = ЕЕ 


ACJ) 
CONTINUE 


DO ¿09 J=1 
@ 12; 
WCJ,2) 
E 
УОЛ?) 
CONTINUE 


Nea 


= ETACGJS 1) 
= 0, 5*DETA(J-1) 


U(J,2)*UES 

UES * W(J,2) 
F(J,2)*SQS*UES 
V(J, 2) /SQS*UES 


UE(NX) = W(1,2) 
UE(NX)*X(NX)*RL 


RX 
SQRX = 


SQRT( RX) 


ТЕСМХ. СТ. МЕБ) САВЕЗЕ 
CALL FILLUP( 2) 


ТЕ(ТКОЕХ. ЕВО СОТО И 


STORE PROFILES AT STAPION NS SPOR” [iE NE oer 


DO 65 J=1,NPT 
FS(J). = HON 
US(J) = U(J,2 
то, 
и) Ето, 
BS(J) = В(Ј,2) 
CONTINUE 

DO 80 J=1,NPT 
F(J,1) = F(J;2) 
ТӨРІ 00202 
WI D 99 WI, 2) 
W(J,1) = W(J,2) 
B(J,1) = B(J,2) 
CONTINUE 
RETURN 

END 


DATA SET KCBCMAIN 


PROGRAM MAIN 


AT LEVEL 005 AS OF 09/18/84 


SUBROUTINE CASBLP(K2,XP,YP,XMP,YMP,XS,YS,DLSP,VC ,DBPP 


,RN,NBL,ITRI,XCTRI,CASEID) 


COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 


COMMON/BLIN/ TITLE(20),XC(100),YC(100),ISG(100),DELS(100), 
XCTR,XTR,ISTRP, ICYCLE, ICYTL, XCTRS( 20m ITRFIND 


COMMON/EDDY1/RL,RX,SQRX,RXNTR,GMTR ,GMTRS(100), 


+ 


COMMON/BLOW/ VN( 100) 
COMMON/TRN/ PGANTR, OMEGA ,RTHETB ,RTRANB 
COMMON/PLOT/NVP(2) ,NXVP( 20,2) , ICC 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


CASEID( 
XP ( 
YME ( 
XS ( 
VT ( 
DLS ( 


20), 
100), 
1009 
100), 
100), 
100), 


ХСТКІ 
ҮР 

VC 

YS 

S 

XO 


98 


( 
( 
( 
( 
( 
( 


100), 
100), 
100), 
100), 
100), 


FERI 
XMP 


NBL 
DLSP 
YO 


ALFAS( 100) , FFS( 100) ,RTSC100) TEDNA E 


(22) 
( 100) 
( 100) 
О) 

( 100) 
( 100) 


INT189 1 
ТКТ 189 2( 
ТКТ 189 30 
INT18940 
INT18950 
INT18960 
INT1897( 
INT1898( 
INT1899( 
INT19000. 
INT19010. 
INT1902( 
INT1903¢ 
INT1904( 
INT19050 
INT19060 
INT1907( 
INT1908( 
INT1909( 
INT1910( 
ІХТ1911( 
INT1912( 
INT1913( 
INT1914( 
INT1915 
INT19161 
INT1917 
INT1918( 
INT1919( 
INT1920( 
INT1921( 
INT1922( 
INT1923¢ 
INT1924( 
INT1925( 
INT1926( 
ІХТ1927( 
INT1928( 
INT1929( 
INT1930( 
INT1931! 
INT1932( 
INT1933( 
INT1934( 
INT19351 
INT1936( 
ІМТ1937( 
INT19 38( 
INT1939( 
INT1940( 
INT1941( 
INT1942( 
INT1943( 
INT1944( 
INT19451 
INT1946( 










Un C2 C2 


Gee 


t5 


50 


60 
C 


3: 


a 


4 


+ 


+ + 


2HNX , 9X, 2HXP, 12X, 2HYP, 12X,1HS,14X, 2HVT) 
15,4Е14. 6) 


DIMENSION 01 (MODA 00—03 (оо) 

DIMENSION XPS (100), ҮР5(100) , DBPP( 100) 

FORMAT ( 20A4 ) 

FORMAT ( 415 ) 

FORMAT ( 315,3F10.0 ) 

FORMAT ( 6F10.5 ) 

PORMAT ( 1H1,5X,20A4 ) 

FORMAT ( 1HO,6X,'CYCLE NO. = ',I5 ) 

FORMAT ( 1H0,6X,'FINISHED TOTAL NUMBER OF CYCLES - ',I5, 
4X,' JOB COMPLETED. ' ) 

FORMAT ( 1H0,6X,'THE UPDATED DISPLACEMENT SURFACE',/1X,2HNX, 
9X,2HXP,12X,2HYP, 11X, 3HDLS, 10X, 4HDBPP, 12X, 2HVN) 

FORMAT ( 15,5E14. 6) 

FORMAT ( 1H0.6X,'READ IN CONTROL POINTS DATA’ ,/1X,2HNX, 9X, 

3HXMP ,11X, 3HYMP, 12X, 2HSM, 13X, 2HVC) 
FORMAT ( 15,4514. 6) 
FORMAT ( 1HO,6X,' INTERPOLATED ORIGINAL GEOMETRY DATA’ ,/1X, 
( 


FORMAT 


mec о, У, та X0)= (X0-X2)%(X0=-X3)/(X1- X2)/(X1- ВЕ 
+(Х0- xij :(Х0- X3)/CX2-X1) /CX2-K3)*Y2+(X0-K1)*CX0-X2) 
/(X3-X1)/CX3-X2)*Y3 

ISWPT 

TEDY 

NBL(1) = 

NBL(2) = 


E 


МОКЕ (К 6,100 ) CASEID 
CONTINUE 

ISTRP = ICYCLE 

Nome К2 - 1 

NT = K2 

КЕШТЕ ( 6,110 ) ІСҮСІЕ 


INTERPOLATE OUTPUT CONTROL POINTS TO ORIGINAL GEOMETRY 


= S(I- 1) + ооо KP SE 


+ Gre Gi) у А о); 


CONTINUE 

SM( 1) = SORT( (XMP(1)-XP(1))**2 Te CY IP AA 5522) 
DOSOO I 2929 ¿NN 

ЭМеТ) = SMÇI-l)+ SQRE( (NMP( IJSNMRGI - 1) 2+ 


+ (ҮМР(І)-ҮМР(1-1))%%2) 


CONTINUE 

CALL AMEAN(NN-10,NN,SM,VC, 1) 
CALL AMEAN(1,NN,SM,VC, 1) 

SNT = S(NT) 


99 


INT19470 
INT19480 
INT19490 
INT19500 
INT19510 
фҚ 19520 
РЫ 19530 
INT19540 
Ji 19550 
INT19560 
719570 
00719580 
19590 
ШТИ 2600 
ENS O O 
ШЕШЕ) 
и 51516; 
INT19640 
INT19650 
Ш 15610 
INT19670 
INT19680 
ый 19690 
INT19700 
ПӘЛЕ ТУ 
ГМТ 192780 
INT19730 
INT19740 
J ji SO 
INT19760 
INT19770 
NT 19 74810 
INT19790 
INT19800 
PIS LO 
PNT 19820 
90198320 
INT19840 
INT19850 
INT19860 
INT19870 
INT19880 
INT19890 
INT19900 
О КО 
ДАТУ ОРАО 
ТЫЛ 19980 
INT19940 
INT19950 
INT19960 
UNT 19940 
INT 19960 
INT19990 
INT20000 
INT20010 
INT20020 


SM(NT) = SM(NN)+SQRT((XMP(NN)-XP(NT) )**24+( YMP(NN) -YP(NT) )9*2) 
SMNT = S(NT) / SM(NT) 
DO 65 1= 1, NN 
65 SM(I) = SM(I) * SMNT 
CALL DIFF3(NN,SM,VC,D1,D2,D3,0) 
CALL INTRP3(NN,SM,VC,D1,D2,D3,NT, S,VT) 


C PRINT OUT INPUT DATA 

C 

C WRITE(6,150) 

C WRITE(6,160) (I,XMPCI),YMP(I),SMCI),VC(I),Iz1,NN) 
C WRITE( 6,170) 

C МКІТЕ(6,180) (І,ХР(І),ҮР(І),5(1),ҮТ(І),1=1,МТ) 

С 


XPMIN = XP(1) 
DO 44 I = 2 , NT 
IF (XP(I) .GT. XPMIN) GO TO 44 
XPMIN = ХР(Т) 
44 CONTINUE 
CHORD = XP(NT) - XPMIN 
DO 45 1 = 1, NT 
XP(I) = (XP(I)-XPMIN) / CHORD 
YP(I) = YP(I) / CHORD 
CONTINUE 
CALL COMPBL ( CASEIDYXP,YPYVTYS,DESP,DES /DBPP, NBL, DERG , XOmram 
+ RN,NT,ISWPT) 


£ 
(л 


SMOOTH THE CALCULATED DISPLACEMENT THINKNESS 
CALL SMFIT(1,NT,S,DLS,D1,2) 


ADD DISPLACEMENT THINKNESS ON THE ORIGINAL BODY 


Са СО СЗ СО СО СО СК? 


DO 70 I 
DLSP(T) 
CONTINUE 
CALL SMFIT (1,NT,XS,YS,D1,2) 
DO 80 т т 
XP(I) = XPS(]) 
YP(I) = YPS(I) 
80 CONTINUE 

IF (ICYCLE .GE. ICYTL-1 .OR. IP .GE. 0) THEN 


1, NT 
DLS(I) 


“J 
O 


б WRITE (6,130) 
C WRITE (6,140) (I,XP(I),YP(I),DLS(I),DBPP(I),VN(I),I-1,NT) 
WRITE ( 6,120 ) ICYCLE 
END IF 
С 
RETURN 
END 
С DATA SET KCBCMAIN2 AT LEVEL 001 AS OF 08/24/84 
C DATA SET KCBCMAIN2 AT LEVEL 001 AS OF 08/24/84 
è DATA SET KCBCMAIN2 AT LEVEL 010 AS OF 04/06/84 


SUBROUTINE MAIN2(ITR,ISWPT,SURFID) 

COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 

COMMON /BLC1/ F(101,2),U(101,2),V(101,2),W(101,2),B(101,2) 
COMMON /BLC2/ DELF(101),DELU(101),DELV(101),DELW( 101) 


100 


INT20030Q 
INT2004( 
ІМТ20054 
ІМТ20060 
INT20070% 
ІХТ20080 


ІМТ20090 
ІКТ20100! 


ІМТ20110 
INT20128 


INT20130 | 
INT 201405) 
INT20150% 


INT2016( 
ІМТ20170 


INT20180 | 
ІМГ2ОТЭ 


INT20200 
INT20210 | 


INT20220% 


ІМТ20230 
ІХТ20240 
INT20250B 
ІМТ20260 
ІКТ2027( 
ІМТ2028С 
INT20291 
INT20304 
ІМТ20310 
INT20321 
INT2033J 
INT2034C 
INT203 55 
INT2036( 
INT2037@ 
ІМТ2038( 
ІМТ20390 
ІХТ2040С 
INT2041( 
INT20420 
INT2043C 
INT20440 
INT2045( 
INT2046f 
INT2047C 
INT20480 
INT20490 
INT20500 
INT20510 
INT2052C 
INT2053J 
INT2054C 
INT20558 
INT20560 
INT20571 
INT20580 





10 


i 


ED 


25 


30 


70 


ШЕЕ T7 BLC//S6(1007 700) ,D( 100) 5 DBCITIOO) , DBP( 100) ,UEO( 100) ,GI 
COMMON /BONV/ ITMAX,EPSL,EPST,CONY 
COMMON/EDDW1/ RL,RX,SQRX,RNNTR , GHTR, GMTRS(I1OO0) 


+ 


r L52300) EPS БОО 5 11S(100) IEDY NXSPT 


COMMON /GTY / X(101) ,UE( 100) ,P1(100) ,P2(100) ,CEL,CELH 
COMMON /SMRY/ VW(100),ITP(100),ISL( 100) ,DLS( 100) ,CF(100), 


+ 


THT( 100) ,NPSTR( 100) 


COMMON /BLIN/ TITLE(20) ,XC( 100)3¥C( 190), 12SG(100) ,DELS( 100), 


4. 
COMMON/BLC9/ UEB(100) 


XCTR,XTR,ISTRP,ICYCLE,ICYTL,XCTRS(2) , TRFIND( 2) 
, CFS(100) 


COMMON/TRN/ PGAMTR,OMEGA,RTHETB , RTRANB 
COMMON /ISURF/ ISF 
COMMON/PLOT/NVP( 2) ,NXVP( 20,2) , ICC 
DIMENSION SURFID(4),RTSS(11) 


LOGICAL SMOOTH , 


SEPART , HOMOPY 


1350 199/70..0.0- 1,0. 2,0. 3,0.4,0. 5,0.6,0. 7,0. 8,0. 9,1. 0/ 


НЕНА X2.X3. Y12Y2, Y35 X0)7 (XO-X2)*0X00X9»/CXI* X25) /CX1 -X3)7*Y1 


+ + 


ISWP 
INDEX 
IGROWT 


*(X0-X1)*(X0-X3)/(X2-X1)/(X2-X3)*Y24(XO-X1)*(X0-X2) 
He SEO cS YS 
0 
| 


2 


ВЕТ = "NAT +*1 
CALL JOINC INDEX) 


NXSTOP = NXT-1 

ШЕШЕ СЕ. МІР ) GOTO 15 
ISWP = ISWP + 1 

NX = NX + 1 

HOMOPY = .FALSE. 

CEL = 0. 5*(X(NX)+X(NX-1))/(X(NX) -X(NX-1)) 
P1(NX) = 0.5 

EBENSO = 0.0 

ШЕШ) - 0 55СЕП 

TER = 0 

CALL COMPGi 

IGROW=1 

LT =IT+1 

RX = UE(NX)*X(NX)*RL 
SQRX = SQRT(RX) 

IF(IT . LE. ITMAX) GO TO 80 


ШИНОМОРҮ) СО TO 72 


TRG Ы 

RT = RTSS(IRC) 
HOMOPY = . TRUE. 
UEREF = UEO(NX-1) 
UESAVE = UEO(NX) 
UEO(NX) = RT*UESAVE+(1. O-RT)*UEREF 
DO 61 J=1,NP 
го, = F(J,1) 
Пака | = U(J,1) 
Сла) 
Wed ,2) = W(J,1) 


101 


INT20590 
INT20600 
INT20610 
Ра орао 
[5 120650 
INT206%0 
INT20650 
INT20660 
INT20670 
INT20680 
INT20690 
INT20700 
INT20710 
INT20720 
TNT207 50 
INT20740 
[11 20750 
INT20760 
TNT 207-70 
INT20780 
INIZ OU 
INT20800 
[51 0510 
INT20820 
INT20830 
INT20840 
INT20850 
INT20860 
INT20870 
INT20880 
[T0297 
INT20900 
О то 
INT20920 
INT20930 
INT20940 
INT20950 
INT20960 
INT20970 
INT20980 
INT20990 
Де 2 1000 
INT21010 
(0021020 
INT2 10590 
INT21040 
NT 21050 
INT21060 
INT21070 
INT21080 
INT2 1090 
INT21100 
INIZIO 
121120 
#121130 
INT21140 


БОЛО) = ВЯ 
61 CONTINUE 
COP TGs 


72 NASTOP NL 


CALL AMEAN(NS ,NXSTOP,X,CF,1) 
(с CALL AMEAN(NS,NXSTOP,X,VW,1) 
С CALL HEADER( TITLE,SURFID,ISTRP ) 


WRITE (6, 250 ) ISWP 
WRITE (6, 260 ) (M,XC(M),X(M),CF(M) ,DLS(M) , THT(M) ,UE(M), 
+  UEO(M),D(M) ,DB(M) ,GMTRS(M) , ITP(M) ,NPSTR(M) ,M=1,NXSTOP) 
WRITE(6, 270 ) NX 
STOP 
80 CONTINUE 
IF(NX .GT. NTR) GOTO 100 
С 
С LAMINAR FLOW CALCULATION 
С 
CALL COEF(GAMMA1,GAMMA2) 
CALL SOLV4(GAMMA1,GAMMA2) 
UE(NX) = U(NP,2) 
IF(ABS(DELV(1)) .GT. EPSL) GO TO 70 


CHECK ON LAMINAR FLOW SEPARATION. IF SEPARATION OCCURS, ASSIGN BEGIN 
OF TRANSITION TO THAT POINT AND RECOMPUTE THE CURRENT STATION NX 


2762 С?З 


IF(V(1,2)-GT.0.0 . OR: TIR- NE >) COOMER 
CALL TRNS(ICODE) 
COTO=25 


TURBULENT FLOW CALCULATION 


Б> С) С С 


00 CONTINUE 
CALL EDDY 
CALL COEF(GAMMA1,GAMMA2) 
CALL SOLV4(GAMMA1 , GAMMA2 ) 
ОЕСМО = Оро) 
VM = AMAX1(V(1,2),1.0) 
ТРАВ (ВЕНУ) ГУМ). СТ. ПЕРО o 


110 CONTINUE 
E 
С CHECK FOR B. L. GROWTH 
С 
IF(NP .GE. NP") GO TO 120 
IF(ABS(V(NP,2)) .LT. 0.0005 .AND. ABS(1.0-U(NP-2,2)/U(NP,2)) 
" . LT. 0. 0035. OR. IGROW. GT. IGROWT) GOTO 120 
CALL FILLUP(1) 
IGROW=IGROW+1 
l MEN 
GO TO 70 
E 


120 CONTINUE 
CALL FILLUP(2) 
CALL OUTPUT(2) 
IFCUNX.GE.NTR .OR. ТТВ. ЕО МОТО пи 


INT2115g 
INT2116§ 
INT2117@ 


19721180 


INT2119€ 
INT2120€ 
INT2 129 
INT21228 
INT2123¢ 


1М721240 


INT21238 
ІМТ2126( 
INT2127€ 
INT2128( 


INT21290 


INT2130( 


INT2131( 
ІКТ21320 
INT2133¢ 


INT2134( 


INT21350) 
ІКТ21360 
INT21370 
18721380 
13721390, 


INT2140€ 


INT21410% 


INT2142( 


INT21430 


INT2144€ 


1№Т21450 | 


INT2146( 
INT2147( 


ТКТ2 1480 


INT2149( 


INT21500% 
INT2151¢ 
INT2152( 


INT2153€ 
ІМТ21544 
INT2155( 
INT2156( 
INT2157Ü 
INT21588 
INT2159¢ 


ІКТ21600 
ІМТ2161{ 


ІМТ2162( 
INT21638 


INT21640) 


INT21658 
INT2166( 


INT2167( 
INT21680. 


INT2 1693 


A 9 p p 





ecc) 


RE 


Cane cane) 


150 


154 


ТН 
CH 
SW 


Б Со С) С С С С 


60 
-- 


+ 


ПРИ Шо . ОК: ITR.NE: 3) GOTOr150 


CALCULATE TRANSITION LOCATION USING MICHEL METHOD 


CALL TRNS(ICODE) 
IF(ICODE.EQ.0) GOTO 150 


TRANSITION OCCURS BASED ON MICHEL CRITERIOR AT STATION NX 


CALCULATE B. L. AT NX STATION ASSUMING THE FLOW IS TRANSITIONAL 


IT 

IGROW 
COTO 70 
CONTINUE 
МИА ЛОТ. HOMOPY ) GO TO 154 

ieee Gl 0.9999) GOTO 154 

IRC = IRC + 1 

RT = RTSS(IRC) 

UEC(NX) =RT*UESAVE + (1. 0-RT)*UEREF 
GO TO 30 

CONTINUE 

ae. LT. NASTOP) GO TO 20 


0 
1 


ЕНІ LG. CALCULATION FOR THE CURRENT SWEEP IS COMPLETED. 
Pere rOR THE CONVERGENCE AND , IT NOT, MOVE TO THE NEXT 
БЕР. 


CONTINUE 
MN). - GRANGCX(NXT-73), XCNXT22) , XCNXT* 1), DCNXT 73) , D(NXT-2), 


D(NXT-1) ,X(NXT)) 


LS( NXT) 
DLS(NXT-1),X(NXT)) 


UE(NXT) 2 GRANG(X(NXT-3),X(NXT-2) , XCNXT- 1) ,UE(NXT-3), 


B UE(NXT-2) ,UE(NXT- 1) , X(NXT)) 


DO 165 I = 1, 
UC CIN = Сг) 
CALL AMEAN(NS,NXSTOP,X, CF , 1) 

CALL AMEAN(NS,NXSTOP, X, VW, 1) 

CALL HEADER( TITLE,SURFID,ISTRP ) 

ПІ ЕСҮСІЕНШШЕ ECYTL21 AND. IP ІІ. 0980 120. 170 

WRITE (6, 250 ) ISWP 

Е Соол) ТСС) СТ) 0ЕСІ) Во) А1) ,0Е(1), 


NXSTOP 


+ UOC), 0. 0,GMTRS( 1) , epee NESTR( 1) 


WRITE (6, 264 ) (M,XC(M),X(M),CF(M) ,DLS(M) ,THT(M) ,UE(M), 


n UEO(M) , DLS(M)/THT(M) , GMTRS(M) , ITP(M) ,NPSTR(M) , M22, NXSTOP) 


800 
810 


170 


IF ((ICYCLE. EQ. ICYTL). AND. CIP. EQ. -2). AND. (NVPCISF). NE. 0)) THEN 
WRITE(12,800) NS+1,NTR,XCTR 
WRITE( 12,810) (XC(M),X(M) ,UE(M) ,CF(M) ,GMTRS(M),ISG(M), 
M=1,NXSTOP) 
FORMAT( 215 ,F10. 6) 
FORMAT(5E15.5,15) 
END IF 
CONTINUE 


DMAX = D(1) 
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CRANE ATAN TAO IES NAT -3) ,DLSCNAT-2), 


ІМТ21700 
ЉУТО 
INT21720 
ТАО 
INT21740 
INT2 SO 
TINT2 1760 
INT21770 
INT21780 
DISP 20 
INT21800 
INT21810 
INT21820 
ENTZ 19030 
INT21840 
INT21550 
INT21860 
Тао ло 
INC? 19909 
ENT2 1590 
ТАП ЭВО 
TODO 
INT21920 
EN D) 
INT21940 
ТОВЭ 
INT21960 
TNT ТОЛО 
ОТОС 
ШШІ2 1920) 
INT22000 
12 2010 
1512202090 
INT22030 
INT22040 
ЈЕ ОВО 
INT22060 
ооо 
INT22080 
INI 22090 
INT22100 
1122110 
INT22120 
INT22130 
INT22140 
INT P0 
INT22160 
TNI221 0 
INT22180 
INT22190 
INT22200 
IMP PPAR 
Is 6 
2280 
157222480 
ГМИІ?2250 


180 


de Өр Шен O) 


190 
195 


200 
22015 


210 


250 


260 


262 


264 
270 


rec) 


DDMAX 
DO 180 
DMAX 
DD 
DDMAX 
CONTINUE 
IP (AR 


II l Hoc | 


ABS(D(1) 
= 2 NXT 

MAX1( DMAX,D(I) ) 
ABS(D(I) - DB(I1)) 
AMAX1( DDMAX,DD ) 


Шер” 


DDMAX / DMAX ) 


UPDATE D FOR THE N БҮ ШЕ 


. LE. 0.0050 ) RETORN 


IF (ISWP TGT тусо TOMIS 
DO 190 I = NS , NXT 
D(I) s D(I)*(1. OXOMEGA*(UE( I) /UEO( I) - 1. 0)) 
со ШЕ 
IF CISKE EO ЕВО 
Do 200 1 = NE 
D(I) = D(I) * (1. 0+OMEGA*(UE(I)/UEB(I)-1. 0)) 
IF(ISWP .GE. ISWPT ) RETURN 
NX = NS 
NP = NPSTR(NX) 
INDEX = 2 
DO 210 I= 1,NXT 
DB(I) = 001) 
UEB(I) = UE(I) 
CONTINUE 
GOTO 10 
FORMAT( 1H0, ' 
1H0,4X,' ISXP Ta 


a 
T 
JB 


+ 


4 
J 


+ 


_FORMATCIHO, &X.2HNX,SX,3HX/C,9X, 1HX, 9X, 2HCF, 8X, 
3HDLS,8X,3HTHT,9X,2HUE, 8X,3HUEO,10X,1HD ,9X,2HDB,3X, 

4HGMTR, 4X, 2HIT, 1X,2HNP/(1H ,3X,13,F10.5,8E11.4,F8. 4,213)) 
FORMAT( 1HO, 4X, 2HNX, 6X, 3HX/C, 11X, 1HX, 10X, 2HCF, 9X, 


FORMAT( 1H 


FORMAT(1HO,' ** ITERATIONS EXCEEDED ITMAX AT NX 
z CALCULATIONS STOP. 


END 


3HDLS ,9X, 3HTHT, 10X, 2HUE , 9X, 3HUEO, 11X, 1HH, 8x, 
4HGMTR,4X,2HIT, 1X,2HNP/(1H ,3X,13,9E12. 4,213)) 


;24412,95087- 4.213) 


1Н0,! 


vege! ) 


=! » 15, 


DATA SET KCBCOUTPUT AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCOUTPUT AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCOUTPUT AT LEVEL 002 AS OF 02/22/84 
SUBROUTINE OUTPUT( INDEX) 


COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 


tt uud 
4% 4. 
5 


COMMON/BLIN/ TITLE(20),XC(100),YC(100), ISGC100),DEELSCIOO), 


XCTR,XTR,ISTRP,ICYCLE, ICYTL,XCTRS(2) , TRFIND(2) 


COMMON /BLC1/ F(101,2),U(101,2),V(101,2),W( 101,2) ,B( 101,2) 
COMMON /BLC7/ C(100,100),D(100),DB( 100) ,DBP( 100) ,UEO( 100) ,GI 


COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR , GHTRS( 100) 
.ALFAS(100),FFS(100),RTS(100), IEDY,NXSPT 
COMMON /GTY / X(101),UE(100),P1(100),P2(100),CEL,CELH 
COMMON /GRD / ETA(101),DETA( 101) ,A( 101) 

COMMON /SMRY/ VW(100) ,ITP(100),ISL(100) ,DLS( 100) ,CF( 100) , THT( 100) , INT22790_ 


NPSTR( 100) 


COMMON /TSURF/ TSE 
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»/ 


,/INT2254Q 

































INT222 
INT222 78 
INT22281 
INT22296 
INT2230g 
INT223 8 
15722324 
1INT2233 
1NT2 234% 
INT22 3% 
ІМТ22360% 
ІМТ22370 
INT223808 
1NT2 23908 
ІМТ22400 
INT22410 | 
INT22428 

INT22438 
INT22440 
1КТ22458 
INT2246( 
INT224 79 
INT22480 
INT2249C 
INT22500 
10722518 
1812252 
1812253 


INT2255¢ 
INT22560. 
INT2257 8 
INT22580 
INT22590. 
INT22600. 
INT22610 
INT22620 
INT2263¢ 
INT2264€ 
INT2265¢ 
INT22660 
INT2267¢ 
INT2268C 
INT22690. 
INT22700. 
INT22710. 
INT227207 
INT22730 
INT22740 
INT22750 
INT22760 
INT22770) 
INT22780. 


-— "ч AA P 


INT22800 
INT22810 


н со с л 


20 


С 
С 
100 


120 


15570] 


С 


COMMON/PLOT/NVP(2),NXVP( 20,2) ,ICC 


ITP(NX) = IT 
NPSTR(NX)=NP 
IF(NX. GT. 1) GOTO 5 
DLS(NX)- 0. 

VW(NX) = 
D(NX) 
THT( NX) 
CF(NX) 
VW(NX) = 
GOTO 150 

GOTO (10,100,200), INDEX 


Il l H ut 
оооооо 
ooo oa © 


CALCULATE B. L. PARAMETERS FOR TRANSFORMED COORDINATES 


CONTINUE 


ОС у - 2.0 * V(1,2) * B(1,2)/SQRX 


VW(NX) = UE(NX)#SQRT(UE(NX)/X(NX))#V(1,2) 
DLS(NX)= X(NX)/SQRX * (ETA(NP)-F(NP,2)) 
D(NX) = UE(NX) * DLS(NX) * 


SQRT(RL) 


01 A (10-01-20) 
SUM = 0.0 

DO 20 J=2 NP 

1% = II C1. 0979092299 
SUM = SUM + ACJ) * (U1 + 02) 
01 = U2 

CONTINUE 


THT(NX)= X(NX)/SQRX * SUM 
GOTO 150 


CALCULATE B. Г. PARAMETERS FOR SEMI-TRANSF COORDINATES 


CONTINUE 
Е SORT (XCNX)) 
SQRL = SQRT(RL) 


CEGNX) = 

VW(NX) = V(1,2) / SQXC 
UE(NX) = UCNP, 2) 

RX = RL * UE(NX) * X(NX) 
DLS(NX) = 

SUM = 0.0 

Ul = 


DO 120 J=2,NP 


U2 - U(J,2)/U(NP,2)*(1.0 -U(J,2)/U(NP,2)) 
SUM = SUM + ACJ) * (U1 + U2) 

Ul - U2 

CONTINUE 


THT(NX) = SUM /SQRL * SQXC 


D(NX) = (UCNP,2)*ETACNP)-FCNP,2)) * SQXC 


ВВ. СЕ. NXT) GO ТО 160 


По ЧЕХ EQ. O .OR. NX . LE. МВ) 6018 260 


2.0 * V(1,2) * B(1,2)/( SQXC*SQRL*W(NP,2)70*2) 


(ETACNP) -F(NP,2) /U(NP,2) ) /SQRL*SQXC 


=" (152) /MGENP, 2) C1. Oe ( 182 )/ 0B 2) ) 


C MODIFY ALFA USING SIMPSON'S ARGUMENTS 


C 


INT22820 
INT22830 
INT22840 
INT22850 
INT22860 
INT22870 
INT22880 
ШЫЛ 43-8 
ENT Z 2900 
Ооо 
225210 
INT22230 
INT22940 
ПО 02727950 
ІШІКІ Зор 
Т2 90 
INT22230 
INT222 720 
INT23000 
TNEZ5010 
J j 0209 
PN ЕО 
INT23040 
ENTZ3056 
INT23060 
INT23070 
INT23080 
ПУЛ 20510 
INT23100 
ШЫП > ІШ) 
P2320 
INTZ23 130 
INT23140 
ГА 
ІМТ23160 
2 0 
BIS 50 
TXEZ T90 
ИЕ 200 
222100 
2) 
о 
INT23240 
плаза 
INT23260 
I 1 22719 
INT 23260 
INT23290 
INT23300 
INT23310 
INT23320 
а 
INT23340 
INT23350 
ТУТ 23050 


CALL SMPSON 


160 
F(J,1) 
U(J,1) 
ШІ) 
MO) 
B(J,1) 
CONTINUE 


175 


ПО 175 l еи 


pe) 
ЛЕТ 22) 
ма) 
ү (ОШОО) 
gu 


IF (CIP. LE. 0). AND. (CIP. NE. -2). OR. CICYCLE. LT. ICI ВО 


C 


C PRINT OUT VELOCITY PROFILES 


200 TF (bee: 
IF (NX. LE. 
FAC] 
FAC2 
ELSE 
FAC1 
FAC2 
ENDIF 
NEST 


— NP 
WRITE(6,4001) 


1) -COTO TO 


NS) THEN 


Је 


SQRT(X(NX) /RL) 


1. 070909] 


Al 


WRITE(6,4000) 


WRITE(6,4100) 


+ 


WRITE(6,4100) 


a 


210 Io (Il ONE: 


NX,X(NX) 


SQRT(X(NX) /RL/UE(NX) ) 
0 


( J,ETA(J) ,F(J,2) ,U(J,2) ,V(J,2) ,W( 7,2) ,BCJ,2), 


ETA(J)*FAC1,U(J,2)*FAC2, J21,NPM1,3) 


NP,ETA(NP) ,F(NP,2) ,U(NP,2) , V(NP,2) ,W(NP,2) ,B(NP,2), 


ETA(NP)*FAC1,U(NP, 2)*FAC2 


-2) 


RETURN 


IF (CCNXVPCICC,ISF).NE.NX). OR. C1CC. GIO NVPCISE ) )) RE 
WRITE( 12,4200) NP 
WRITE( 12,4300)" (ETACH ВЫ 
WRITE( 12,4300) СП J 2 = D) 


ICC = 
RETURN 
4001 
4000 


ТОТ 


FORMAT(/1HO, МХ = ‚15, ' 
FORMAT(1H0,2H J,9X,3HETA, 15X, 1HF ,13X, 1HU,13X, 1HV, 13x, 1HW, Doe 


S/C =',F10.5) 


+ 13x; 3HY/G 10 X АН ДЕ 


4100 
4200 
4300 


FORMAT(IS) 


END 


С СУС? 


FORMAT( 8F10. 6) 


FORMAT(1H ,13,E14.5,2X,5E14.5,2X,2E14.5) 


DATA SET KCBCSMFIT AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCSMFIT AT LEVEL 001 AS OF 08/24/84 
DATA SET KCBCSMFIT AT LEVEL 001 AS OF 08/15/83 


SUBROUTINE SMFIT(NS,ND,X,Q,D,KS) 


NS 


Cec cy C2 C) C2 C2 


KS NO OF 


: BEGINNING POINT? 
X : INDEDENPENT COORDINATE? D 
Q : VARIABLE TO BE SMOOTHED 


SMOOTHING 


NI: z: 


THIS SUBROUTINE SMOOTHES DATA, Q, USING FIVE-POINT FORMULA. 


END POINT 
WORKING STORAGE 


DIMENSION X(101),Q(101),D(101) 


С) 


5М15(01,02,03, 04 ,05 


ОСОТ ООС) 


) 


0.0625*( 10. 0*03+4. 0%(092+04)-01-05) 
0.5%(02%(01%АВ5(Х3-Х2)%03%АВ5(Х2-Х1)) 
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INT23370 | 
INT2338( 
INT2339( 
INT2340( 
INT23410 | 
INT2342( 
INT2343( 
19723220 
INT2345( 

INT23460 
INT2347( 
ІМТ23480 
INT234908 
ІМТ23500 
INT23510% 
INT2352¢ 
ІМТ23530 
INT2354( 
INT23550% 
ІХТ23560 
INT2357] 
INT2358( 
INT2359( 

INT2360( 
INT2361( 
INT23620 | 
INT23630m 
INT23640. 
INT2365( 
INT23660 | 
INT23670 | 
INT23680 | 
INT23690 | 
1172 37007 
INT237108 
INT23720 | 
INT23730 | 
INT23740 | 
INT23750 
INT23760 | 
INT237708 
INT23780.. 
INT23790. 
INT23800 © 
INT23810 
INT238208 
INT23830 
INT23840 | 
INT23850 
INT2 3860 
INT23870 
INT23880 
INT23890. 
INT23900 
INT23910 
INT23920 




























20 


40 


100 


200 


220 


250 
590 


moo 


ЖАБЫСАЗ- ХІУ) 
Во. БЕ ӨШ КЕТПЕК 


МР1 = NS+1 
NSP2 = NS+2 
№01 = ND-1 
NDM2 = ND-2 
КОТЕ = ND-NS+1 


Ш СМОТР Г. З) KETURN 
Пи . 27. 5 ) 60 ТО 200 


DO 100 K=1,KS 

D(NS+1)= SMT3(Q(NS) ,Q(NS*1) ,Q(NS*2) ,X(NS) , X(NS* 1) ,X(NS2)) 
D(ND-1)= SMT3: Y(ND-2),Q(ND-1) ,QCND) ,X(ND-2) , X(ND- 1) ,XCND)) 
DO 20 IsNSP2,NDM2 

ПОНИ = 31T5090C1-2),QCI-1), QC I) ,QC I1) , QC I*2)) 

CONTINUE 

PO 40 I=NETi,NDM1 


ОСТ) DM 
CONTINUE 

RETURN 

DO 300 K = 1,KS 

DO 220 I = NSP1,NDM1 
ПО = ОСОО ООК СТ) ХО), Х(І+1)) 
CONTINUE 

DO 250 I = NSP1,NDM1 
ОСІ) = DCI) 
CONTINUE 

CONTINUE 

RETURN 

END 


DATA SET KCBCSOLV3 AT LEVEL 001 AS OF 08/24/84 

DATA SET KCBCSOLV3 AT LEVEL 001 AS OF 08/24/84 

DATA SET KCBCSOLV3 AT LEVEL 005 AS OF 02/21/84 
SUBROUTINE SOLV3 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
CUNNEN. /BUcT/ FC 101,2),UC101,2),V(101,2) ,W(101,2),B(101,2) 
COMMON /BLC2/ DELF(101),DELU(101),DELV(101) ,DELW( 101) 


COMMON /BLC6/ 51(101),82(101),83(101),84(101),85(101),56(101), 


222920 
INT23940 
тра 950 
Pai 23960 
ШЫ 39270 
15325990 
Т2 59010 
1724000 
INT24010 
INT24020 
INT24030 
INT24040 
INT24050 
INT24060 
INT24070 
INT24080 
INT24090 
INT24100 
INT24110 
INT24120 
INT24130 
INT24140 
INT24150 
INT24160 
INT24170 
INT24180 
INT24190 
INT24200 
INT24210 
INT24220 
INT24230 
INT24240 
INT24250 
INT24260 
INT24270 
INT24280 
INT24290 
INT24300 
INT24310 
INT24320 
INT24330 
INT24340 


S7(101),S8(101),R1( 101) , R2( 101) , R3( 101) , RA( 101) INT24350 


COMMON /GRD / ETA(101),DETA(101),A(101) 
COMMON /BLCB/ A11(101),A12(101),A413(101),A14(101), 
A21(101),A22(101),A23(101) ,A24(101) 


А11(1)= 1.0 
А12(1)= 0.0 
A13(1)= 0.0 
КОТ о.о 
Ето 
о 
611=-1.0 
НЕА 2) 
ШЕ- 010 
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INT24360 
INT24370 
INT24380 
INT24390 
INT24400 
INT24410 
INT24420 
INT24430 
INT24440 
INT24450 
INT24460 
INT24470 
INT24480 


G21= S4(2) 
G@G=-S2(2)/AC 2) 
622= 623+56(2) 


FORWARD SWEEP 


(59727202 


DO 500 J=2,NP 
IF ВО 2) СОЛО, 


108 











ІХТ24490 
ІМТ24504 
ІМТ24511 
ТКТ24 5 
INT2453f 
INT2454]f 
INT24558 
1ЧТ2456 


ПЕМ = (A13(J-1)*A21(J-1)-423(J-1)*A11(J-1)-A(J)* INT2457( 
т (A12(J-1)*A21(J-1) -A22( J- 1)*A11(J-1))) INT24581 
DBNi = А22(1-1)5 1) А (0-1; INT2459( 
611= (A23(J-1)*A(J)*(A(J)*A21( J- 1) -A22( J- 1)) ) /DEN INT24601 
C122 -(A(J)*ACJ)4611*(A12(J-1)*A(J) -A13(J-1)) )/DEN1 INT2461( 
G13= (G11%A13(J- D MOLD GENE NOR INT2462( 
G21= (S2(J)wA21(J-1) S&I FRJ D еу: INT24631 
i A22(J-1)-S6(J)*A21(J-1))) /DEN INT2464( 
G222 (-S2(J)486(J)*A(J) -G21*( AC J)*A12( J-1)-A13(J-1))) /DENI INT2465( 
G23= G21*A1 2 J угш ЗАД пера ee) INT2466( 
100  A11(J)= 1.0 INT2467¢ 
A) MEG INT2468( 
A13(J)= A(J)*G13 INT2469( 
АО ssp INT2470¢ 
A22(J)= S5(J)-G23 INT2471( 
A23(J)= S1(J)+A(J)*G23 INT2472¢ 
R1(J) = R1(J)-(G11=Ra@g- Wee" R2( Jancis elem INT24730% 
R2(J) = В2(Т)-ГЕ2 ВЕТ ое оО) 1722740 
500 CONTINUE INT247507 
C INT24760 
C BACKWARD SWEEP INT24770 
C | INT24780F 
DELU(NP) = R3(NP) INT24790 
El = R1(NP)-A12(NP)*DELU(NP) INT24800 ` 
E2 = R2(NP)-A22(NP)**DELU(NP) INT24810 | 
DELV(NP) 2 (E2*AI1(NP)-E1*A21(NP))/(A23(NP)*A11(NP) -A13(NP)** INT24820 | 
а A21(NP)) INT24830 | 
DELF(NP) = (E1-A13(NP)*DELV(NP) )/A11(NP) INT24840_ 
J = NP INT24850 
600 J = J-1 INT24860 | 
E3 = R3(J)-DELU( J+1)+A( J+1)*DELV( J+1) INT24870) | 
DEN2 = A21(J)*A12(J)*A(J-1)-A21(J)*A13(J)-A(J-1)*A22(J)*A11(J)* INT24880) 
Е A23(J)*A11(J) INT24890 | 
DELV(J) = (A1103)*(R2(J)+E3*A22(J))-A21(J)*R1( J) -E3*A21(J)*A12(J)INT24900 
~ ) /DEN2 INT24910 
DELU(J) =-A(J+1)*DELV(J)-E3 INT24920 
DELF(J) = (R1(J)-A12(J)*DELU( J)-A13(J)*DELV(J))/A11(J) INT24930 
IF(J .GT. 1) GO пе" 600 INT24940 
C INT24950 
DO 700 J=1,NP INT24960 
F(J,2)= FC(J,2)+DELF(J) 137224970. 
U(3,2)= U(J,2)4DELU( J) INT24980 
V(J,2)= V(J,2)+DELV(J) INT24990 
700 CONTINUE INT25000 
U Ee INT25010 
CALL EDGCHK(NP,ETA,F(1,2),U(1,2),V(1,2)) INT25020 
RETURN INT25030 
С----------------------------------- 15725040 




















Coen) 


10 
C 
C 


END INT25050 
DATA SET KCBCSOLV4 AT LEVEL 001 AS OF 08/24/84 INT25060 

DATA SET KCBCSOLV4 AT LEVEL 001 AS OF 08/24/84 INT25070 

DATA SET KCBCSOLV4 AT LEVEL 001 AS OF 02/21/84 INT25080 
SUBROUTINE SOLV4(GAMMA1,GAMMA2) INT25090 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP INT25100 
COMMON /BLC1/ F(101,2),U(101,2),V( 101,2) ,V( 101,2) ,B( 101,2) INT25110 
COMMON /BLC2/ DELF(101),DELU(101),DELV(101),DELW( 101) INT25120 
COMMON /BLC6/ S1(101),$2(101),$3(101) ,$4(101) ,$5(101) ,$6(101), INT25130 
+ S7(101).S8( 101), R1C 101) ,R2( 101) ,R3( 101) ,R4(.101) INT25140 
COMMON /GRD / ETA(101),DETA(101),A(101) INT25150 
COMMON /BLCB/ А11(101),412(101),413(101),А14(101), INT25160 
n A21(101),A22(101)99225(101) 4224( 101) INT25170 
INT25180 

сз О -—-—--:-2------- INT25190 
INT25200 

ПШ ШШ = 1.0 INT25210 
ШЕГІ) -о00 INT25220 
NEG) = 0.0 INT25230- 
NEU = 0.0 INT25240 
Mei) =0.0 15725250 
fee) -10 INT25260 
feed) -00 INT25270 
A24(1) = 0.0 INT25280 
DO 10 J = 2,NP INT25290 
AAl = A13(J-1)-A(J)#A12(J-1) INT25300 
AA2 = A23(J-1)-A(3)*A22(J-1) INT25310 
AA3 = §2(J)-A(J)*S6(J) INT25320 
DET = AA2*A11(J-1)-AA1*A21(J-1) INT25330 
AJS = A(J)#*2 INT25340 
611 = -(AA2+A21(J-1)*AJS) /DET INT25350 
EDS — (Al1(J-1)*AIS+AA1) /DED INT25360 
613 - Al2(J-1)*G114+A22(J-1)*G12+A(J) INT25370 
614 = A14(J-1)*G11+A24(J-1)*G12 INT25380 
621 = (S4(J)*AA2-A21(J-1)*AA3)/DET INT25390 
622 - (А11(1-1)%442-54(7У%ДАТУ/ПЕТ INT25400 
623 = Al2(J-*)*G214+A22(J-1)*G22-S6(J) INT25410 
ЕЙ) - 114(1-1):621:424(1-1)%622-58(7) INT25420 
ЕНІ) = го INT25430 
А12(Ј) = -*1”)-613 INT25440 
Ai3(J) = A(J)*G13 INT25450 
414(3) = -G14 INT25460 
421(Ј) = 83(Ј) INT25470 
А2207) = 55(Ј)-623 INT25480 
АТ) -- SICJOTÀCJ)*G23 INT25490 
424(1) = S7(J)-G24 INT25500 
R1(J) = В1(Ј) -G11*R1(J-1)-G12*R2(J-1)-R3(J-1)*G13 13725510 
+ -G14"R4(J-1) INT25520 
R2(J) = R2(J) -G21*R1(J-1)-G22*R2(J-1)-R3(J-1)*G23 INT25530 
E -G24"R4(J-1) INT25540 
CONTINUE INT25550 
INT25560 

BACKWARD SWEEP INT25570 
J = NP INT25580 
Gi = GAMMA1/GAMMA2 INT25590 
R3(J) = R3(J)/GAMMA2 INT25600 
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КАСЈ)-А12(Ј)“(КЕСЈЈЊЖЕЗСЈ)) - Ала Јука 


R2(J) = R2(J)-A22(J)*(R4(J)+R3( J) )-A24(J)*R3( 0) 
Gl = A11CJ)-61*0A12CJ)9 ATA ИП) 
C2 = A21(J)-G1*(A22(J)4A24( J) ) 
DET = Cira c NE 
DELF(J) = (R1(J)*A23(J) -R2( J)*A13(J) ) /DET 
DELV(J) = (C1*R2(J)-C2*R1(J))/DET 
DELW(J) = R3(J)-G1*DELF(J) 
DELU(J) = R4(J)+DELW(J) 
20 1 = 7-1 
дей = DELU( J+1)-R3(J)-A( J+1)*DELV( J+1) 
пей, = DELW( J+1)-R4(J) 
CC3 = A13(J) AG B rT) 
CC4 = R1(J)-A12Z2(J)"6e1-Al4(J)*ce2 
ces = A23(J)-A(J+1)%*A22(J) 
ССё = R2(J)-A220 ШЕТ 
DENO АТС) СЕО 
DELF(J) = (CC4*CC5-CC3*CC6)/DENO 
PELV(J) = (A11(J)*C0C6-A21(3)*CC4)/DENO 
DELW(J) = 662 
DELU(J) = CC1-A(J+1)*DELV(J) 
те Се № 20 10 20 
DO 30 J = 1,NP 
F(J,2) = F(J,2)+DELF(J) 
(7,2) - ШТ, 2У ПЕШІ 4%) 
УС]. 2) = VO ВЕТ 
W(J,2) = W(J,2)+DELW(J) 
30 CONTINUE 
(О 1527 ТОО 
CALL EDGCHK(NP,ETA,F(1,2),U(1,2),V(1,2)) 
RETURN 
E = = = = we = - = = = = = = Сы и a a Uc 
END 
Ç DATA SET KCBCTRNS AT LEVEL 001 AS OF 08/24/84 
C DATA SET KCBCTRNS AT LEVEL 001 AS OF 08/24/84 
C DATA SET KCBCTRNS AT LEVEL 005 AS OF 03/13/84 
SUBROUTINE TRNS(ICODE) 
E 
C CALCULATE TRANSITION LOCATION USING MICHEL CRITERION 
6 
COMMON /BLCO/ NX,NXT,NP,NPT,NTR,IT,INVRS,NS,IP 
COMMON /BLC1/ F(101,2),U(101,2),V(101,2),W(101,2),B(101,2) 
COMMON/BLIN/ TITLE(20),XC(100),YC( 100), I$G( 100) ,DELS(100), 
+ XCTR,XTR, ISTRP,ICYCLE,ICYTL,XCTRS(2) , TRFIND(2) 
COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR , GMTRS( 100) 
+ ,ALFAS( 100) ,FFS( 100) ,RTS(100) , IEDY ,NXSPT 
COMMON /GRD / ETA(101),DETA(101) ,A(101) 
COMMON /GTY / %(101),UE(100) ,P1(100),P2(100) ,CEL,CELH 
COMMON/TRN/ PGAMTR, OMEGA, RTHETB , RTRANB 
Ф - m < = s ss ав ss - = oe аз аз = = а а с а ме = oo 6 =] os = dm —--_" ___-_ __ eee 


100 FORMAT(/3X.'DOSIN OF TRANSITION HAS BEEN DETECTED BY MICHEL''S ', 
+'CRiTERION: X/C =',F8.4,4X,'S/C =',F8.46,4X. NOR =e) 

110  FORMAT(/3X,'BEGIN OF TRANSITION IS ASSUMED AT THE POINT OF ', 
+'LAMINAR SEPARATION: X/C =',F8.4,4X,'S/C =',F8.4,4X,' NTR =' ,13/) 
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С eee lS 


aa aa 


77272 


10 


929272 


ehe 


^^ 


20 


30 


TSEP 


il 


IE 2) РТ. 0.0) THBN 
TRANSITION PROCESS HAS BEGUN DUE TO LAMINAR SEPARATION *** 


FAC = 


GOTOSZ 


END IF 


в 


CHECK MICHEL'S TRANSITION CRITERION *** 


ISEP 
SUM 
F1 = 


0.0 
Vel UNE leo UC 1,2) /UCNP,Z)) 


Јо Ј=>2 МР 


p 
SUM 
li 
CONTINUE 
CONV = 


I NA. LE. 


THETA = 
RTHETA = 
RTRAN = 


UCJ 229 4UCNP 925: C109020(J, 2) /UCNP, 2) ) 
SUM + (F1 + F2 ) *A(J) 
F2 


SQRT(RL/X(NX)) 

NS) CONV = SQRT(RX)/X(NX) 

SUM / CONV 

RL * UE(NX) * THETA 

1.174 * (1.0+22400. 0/RX) * RX**0. 46 


IF(RTHETA. LT. RTRAN) THEN 
RTHETB = RTHETA 
RTRANB = RTRAN 


RETURN 


END IF 


TRANSITION PROCESS HAS BEGUN BECAUSE OF MICHEL'S CRITERION *** 


FAC = (RTHETB-RTRANB)/(RTRAN-RTRANB -RTHETA+RTHETB ) 


COMPUTE EXACT LOCATION OF TRANSITION BEGIN хе“ 


NTR = NX-1 

NTR1 = NIR + 1 

ХСТЕ = XC(NX-1) + FAC*(XC(NX)-XC(NX- 1)) 

XTR = X(NX-1) + FAC*(X(NX)-X(NX-1)) 

UETR 2 UE(NX-1) 9 FAC*(UE(NX) -UE(NX-1)) 

IF ( ISEP .EQ. 0 ) WRITE (6,100) XCTR,XTR,NTR 
IF ( ISEP . EQ. 1 ) WRITE (6,110) XCTR,XTR,NTR 
ICODE = 1 


RANTR = 
GGFT = 


сен CALCULATE INTERMITTENCY DISTRIBUTION w 


XTR ** UETR * RL 
RL**:2 / PGAMTR/RXNTR**1. 34**UETR**3 


DO 30 I=NTR1,NXT 


ALFAS(I) 


— 0.0168 


GMTRS(1)= 1.0 


CONTINUE 


ALFAS(NTR) = 0.0168 


UEINTG = 
01 = 
Х1 


XIR 


0.0 


ОРАО 


DO 40 I=NTR1,NXT 


U2 = 


VEINTG = 


Па ЈЕСТ) 
X2 = X(T) 


UEINTG4(U14U2)*(X2-X1) 
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40 


NE EDGEGHKCNP, ETA, в, ПО 
М ЕТА( 101), F( 101), UCG101) 0 0 0) 
NP 3 
ND EST 
JS, NPM1 
J 


NOTES: 
E 
2 


8 6S 6s 4N 4S 0S Фъ 


U2 
МА 
GGF 
ПИ Б 
GMTRS(I) 
CONTINUE 


ж 
| 
Hott IHH 


verte RESET FI 
50 


DO 60 J= 
ШО?) 
U(J,2) 
мт.) 
B( J,2) 


SUBROUTI 
DIMENSIO 


JS 
NPM1 
DO 10 J= 
T - 
ЖЕ AU 
CONTINUE 
RETURN 
JS = JJ 
lease cal 
CALL AME 
CALL AME 
DETAP = 
VJP = 
DO 30 J= 
DETAM 
VIM 
V(J) 
VJP 
DETAP 
CONTINUE 
V(NP) = 
RETURN 


PT LP m m Pm 
0% 7% 4% 4% 4% 


' EDDY' 


T*UEINTG*( XCI)-XTR) 
T. ТОО) СОВИ 
= 1. ОКР ЕО) 


NITE DIFFERENCE CALCULATIONS «ки 
ПИВА 

Е) 

UJ ы) 

Ve 

ВОЛИН) 

we s 


СЕ. КР) -OR УС. ЕЕ о она 


- 1 

„(КР-2)) Ј5 = МР-2 
AN(JS, NP, ETA, U, 1) 
AN(JS, NP, ETA, F, 1) 
ETA(JS) -ETA(JS-1) 
(UCJS) -U(JS-1))/DETAP 
JS,NPM1 

ETA( J+1)-ETA(J) 
(UCJ41)-U(J)) /DETAM 
(VJM*DETAP - VJP*DETAM) /( DETAP+DETAM) 
VJM 

DETAM 


m 
-V(NP-1) TZ Б 
eu oJ a nia n'a nl anta nta nte ufa nte uta efe efe afa ufa sta uta 
Фу ER EV ED ED ER EX ER ERX FR FR 2% въ къ бъ 4\ 


(FOR CHANGING FROM THE ORIGINAL PROGRAM) 


HAS BEEN MODIFIED BY ADDING 'FINT'. 


SUBROUTINE 'EDGCHK' HAS BEEN ADDED. 
3. GROWTH LIMIT HAS BEEN ADDED FOR 2 IN 'MAIN2'. 


ela nana ufa n'a nta ufa ufa ufa nia els ufa efe etant a ufa nta Je ъв ва а ъв ива а а uta eta ufa uta ufa uta afa afa uta ufa cfe ufa 
ву «ъ «Фу бъ гу 27 27“ Fe FR EL ER ER FR ER EX 06S FR FN 2% ov de er ed 2% Фу Фу чъ EN ER ER GR ER GER EB FX GY 


? 
че 
~“ 
2% 
eto 
“` 
ale 
гу 
of, 
77 
ales 

4% 


ale 
oe 


.... ~". 
т. 
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INT26950 
INT26960 
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INT27100 
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INT27120 


INT271309 


INT27140 
INT27158 
INT27160 
INT2717@ 


INT271808 
INT2719@ 
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INT2 72208 
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1М727260 
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OOOO aaa MAA CARACAS 


30 


E 


38 


C 


SUBE QUTINE SMPSON 


COMMON /BLCO/ NX,NXT,NP,NPT,NTR, IT, INVRS,NS, IP 

НІНІ 101) Е(101,2).1(101,2),У(101,2),М4(101,2),8(101,2) 
COMMON /BLC7/ C(100,100),D(100),DB(100),DBP( 100) ,UEO( 100) ,GI 
COMMON/EDDY1/ RL,RX,SQRX,RXNTR,GMTR , GMTRS( 100) 


n ,ALFAS( 100) ,FFS( 100) ,RTS(100) , IEDY ,NXSPT 


COMMON /GTY / X(101),UE(100),P1(100) ,P2(100) ,CEL, CELH 
COMMON /GRD / ETA(101),DETA(101) ,A(101) 

DIMENSION CRD(12),RTD(12) 

We RTD/O. 00,0. 05,0. 12,0. 20,0. 30,0. 40,0. 50,0. 60,0. 70, 


+ 0. 80,0. 90,1. 00/ 


JO ORD/5. 00,4975, 4. 3545980, 3/2552. 85,2. 58,2. 37,2. 25, 


+ 2215 2 o PIS 


ШЕР 1 CALCULATE (DU/DX)/(DU/DY) 


ШАМА ГТ. КХӘРТ) 6070 10 


ІМ THE SEPARATED REGION, ALFA SET TO BE CONSTANT 


ALFAS(NX)- ALFASP 
RETURN 


10 CONTINUE 


КОТ а). ст. 0.0) GOTO 20 


SEPARATION OCCURS. ALFA SET TO BE THE PREVIOUS ITERATED VALUE 


ALFASP = ALFAS(NX) 
NXSPT = NX 
RETURN 
MODIFY OUTER EDDY BASED ON SIMPSON SUGGESTION 
TM = 0.0 
IM = 1 
DO 30 J=2,NP 
TS = (B(J,2)-1.0)* V(J,2) 
ИН ӨС [T TH) COTO 30 
TM = TS 
JM = J 
CONTINUE 


ИЕ 99550 УС7М,22ҒУ(7М,1)) 
ПАРА LE. NS) GOTO 35 
PURSUIT, 2) =UCIMS 17). / ARENA) -X(NA-=1)) 


GO TO 38 
DUDX = CEL*(U(JM,2)-U(JM,1))+P2(NX)*U(JM,2)+0. 5*ETA(JM)* 
VNXM*( P2(NX) -1. 0) 
RU = RL 
IF(NX. LE.NS)RU = RL * UEO(NX) * X(NX) 
RL2 = SQRT(RU) 
RR = DUDX/VNXM/RL2 
STEP 2 : CALCULATE (UU - VV)/UV 


VNXM 7 0.5*(V(1,2)*V(1,1)) 


127260 
КІП 290 
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ТАЙ 2550) 
INTZ7 560 
1227520 
PNT 27550 
0127590 
INT27600 
INT27610 
INT27620 
PN 27650 
INT27640 
1NT27650 
INTZ7 060 
INT27670 
ENT CeO 
IND 27690 
INT27700 
ТО 
Ты 20 
PNT2 7780 
INT27740 
M27750 
INT27760 
INT2777@ 
NNT 2 7950 
INT2 7790 
1.127500 
JN OVO 
INT27820 


IF(FR .GT.0.35) FR = 0.35 
ТТ 0050-2-02 
FFS(NX)= (FFS(NX) + (I 0 FR) 


ALFAS(NX)- 0.0168/FFS(NX)**2. 5 


RT = VNXM/TM 
C 
IFCRT. GT. 1,0) С@ тен 
CR SOT 
C CR = 2.0 
C DO 40 I=2,1 
C IF( RT. LTƏRT 
C 40 CONTINUE 
C GOTO 70 
C 50 CE 
(10402-70 
60 CR 
C 
C STEER 3. >” CALCULATE Fr 
70 FR = CR * RR 
IF (FR 
RTS(NX)= RT 
RETURN 
(C (=-= 2 
END 
C 


SUBROUTINE XSPACE(NINRITE q, ALET, RAD, NEL INRI) 
DIMENSION XII( 200) ,T(200) 
DATA PIS. 141592653596 737 


RAD = 


NLEFT=NI-1-NRITE 


NRG=NR 


IFC CNRITE/2*2) „МЕ. ВВЕЛА l > 


РК1МТ' (3Х,215,3210.3) NX JM ЕВ 
ТЕ (ВТ „БТ. 00605 е Па 


= CED. i-1)+(CRDCI)-CRD(I-1))*(RT-RTD( I-1))/CRTDC1) -RTD( 1-1) ШЫ 
= (1.0 + @ /RT 


PI 


14702727) 


NL1=NR4+1 
NL2=NR4+NLEFT 


NR1I=NL 
NR2=NI 


ZE 


Р12-0. 5“РІ 


ҚАр2-(РІ-КАр)/2.0-РІ2 


RAD3=RAD2+RAD 


SRT =RAD2/FLOAT(NR4) 


SRT2=SRT 


IF((NRITE/2*2) .NE. NRITE) SRT2=RAD2/FLOAT(NR4-1) 
SLT = RAD/FLOAT(NLEFT) 


DO 10 


10 XII(I)=0. 5*(1. 0+COSCFLOATCI I ERE) 


DO 20 


20 XII(1)=0. 5+XLLT*COS(FLOAT(I-NL1)*SLT+RAD2) 


DO 30 


30 X1II(1)=0. 5*(1. 0+COS(FLOAT( I-NR1)*SRT2+RAD3)) 


Iz1,NR4 
I=NL1 ,NL2 


I=NR1,NR2 


NA=(NI+1)/2 


IF((NI/2*2) .EQ. NI) NA=NI/2+1 


FN1=FLOAT(NA-1) 


FN2=FN 


IF((NI/2*2) .EQ. NI) FN2-FLOAT(NA-2) 


DO 40 


40 T(I)=FLOAT(NA-1)/FN1 


1 


I=1,NA 


¡(150 2.0 ~IR F p 


2 
RTD(i)) GOTO 50. 
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C 


(72 


PEP P PLP 


60 


62 
65 


68 


ТЕ 
B 


CALL AMEAN(1,NA,T,XII,1) 
ES AIICI) - AII(C2) 
DEOR LI. 0.004) THEN 
DO 55 1-2,5 

iene) = XLICI-1)=xDIF*3.0 
CONTINUE 

CALL AMEAN(2,NA,T,XII,10) 
КОБ ТЕ 

DO 50 I=NA,NI 

=> XITI(NI-ITI) 
RETURN 

END 


ns I INE TRGRID (N1 , KO , YO,NI,NRITE,XLLI,N10,RAD, ID,NXSS) 
C THIS SUB. IS TO REGRID SPACING NEAR TRAILING-EDGE 


DIMENSION X0( 200) ,YO( 200) ,X1(200) ,YI(200) ,D1(200) ,D2( 200) ,D3( 200), 
X00( 200), Y00(200) ,XII(200) , YII( 200) ,WX( 200) ,WY( 200) , 
WXI: 200),WYI(200),T(200) 


N20 = N10 
ШЕГГГІІ7252) ЕС NI) N20 = N10+1 


IF((NI-(Ni,2)*2) .NE. 0) Nli= (NI-1)/2+1 


IF(CNI-(NI/2)*2) .EQ. 0) NIIZNI/2 
№21 = Nil 
IF((NI/2*2) .EQ. NI) N2I = N11+1 


CALL XSPACE(NI,NRITE,XI,XLLT,RAD,NL1,NR1) 
LEN NRITE- ,NRITE, XLLTsS ,XLLT 


WRITE (6, 290) 

WRITE (6, 300) (XO(I) ,I=1,N1) 
WRITE (6, 298) 

WRITE (6, 300) (YO(I) ,I=1,N1) 
IF(ID .EQ. 2) THEN 

DO 60 I=NL1,NXSS 

YI(I)=YO(I) 

XI(I)=X0( I) 

NXST=NXSS+8 

XM1=( XI(NXST) -XI(NXSS) )/8. 0 
XM2=(XI(NR1)-XI(NXSS))/(NR1-NXSS) 
ро 62 I=NXSS,NR1-1 

pone y=( XO( 1 )tX1(1))/2. 0 

DO 65 I=NXSS,NXST 

T( 1)=FLOAT( I-NXSS)*XM1+XI(NXSS) 
CALL AMEAN(NXSS,NXST, T, XI,8) 

DO 68 I=NXSS,NR1 

T( 1)=FLOAT( I -NXSS )**XM2+XI(NXSS) 
CALL AMEAN(NXSS,NR1,T,XI,12) 
N10=NL1 

N1I=N10 

N20=N1+1-NXSS 

N2I=NI-NXSS+1 

ELSE 

NXSS=N21 

END IF 


C FOR LOWER SURFACE 


DO 5 I1=1, N10 


INT28390 
INT28400 
INT28410 
INT28420 
INT28430 
INT28440 
INT28450 
INT28460 
INT28470 
INT28480 
INT28490 
INT28500 
INT28510 
| 123520 
J 1265 О 
INT28540 
INT28550 
INT28560 
DNT2 359 
INT28580 
182 8590 
INT28600 
INT28610 
525620 
INT28630 
INT28640 
INT28650 
INT28660 
INT28670 
INT28680 
INT28690 
INT28700 
[50/25 7 LO 
INE2 8720 
INT28730 
INT28740 
ENG 2 37 50 
INT28760 
ГОВ О 
INT28780 
Th i287 30 
INT28800 
INT28810 
INT28820 
INI 55 O 
INT28840 
INT28850 
INT28860 
INT28870 
INT28880 
INT28890 
INT28900 
INT28910 
INTDZ5920 
INI295990 
INT28940 
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с» 
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ТГ = МО ет 

UXCII 729 

СЕ = Your 

CONTINUE 

DO 7 I = 1 , N1I 

ТТ = №11 - 1 + 1 

и г; 

CONTINUE 

CALL DIFF3(N10,WX,WY,D1,D2,D3,0) 
CALL INTRP3(N10,WX,WY,D1,D2,D3,N11,WXI,WYI) 
DO9I=1, NiI 

ТТ = №11 - I+ 1 

CIO 

CONTINUE 


FOR UPPER ScURPAGKH 


DO 10 IT =1 , N20 
ІЛ = № - ОИ 


КОО ЕК 20... 
YOO(I) = YO“%IT) 
CONTINUE 


DO 20 I = 1 , N2I 

J Ni Io] 

WXI(I) = XI(II) 

CONTINUE 

CALL DIFF3(N20,X00, MONADA TD2 DS.0)) 
СА, ІКТКРЗ(М20,Х00,Ү00,Ы1, 84990521, МХ1 71) 
DO 25-T = ТЕК 

EL — N2I.- IS 

XIICII) = WICI) 

YII(II) = WYI(I) 

CONTINUE 


COMBINE TWO SURFACES INTO ONE CIRCLE 


NN = №! - №10 - №20 
DO 30 I =1 , NN 

E и 

111= N10 + I 

XI(II) = XOCIII) 
“ш = ОО, 
CONTINUE 
ро 40 1 = 1 
Il = NXSS + I -1 
12 = Коен A 
XII) = юш; 
УСТ) = ра) 
CONTINUE 

au) = E 
о 
YI(1) = YO(1) 

О ОВ 


N1 = 11 
DO 50 I 


= 1 , N1 
XOCI) = XI(1) 
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eo GaGa) O° Ca Un 


CONTINUE 
WRITE (6 , 295) 
WRITE (6 , 300) (XO(I) , I=1,N1) 
WRITE (6 , 298) 
WRITE (6 , 300) (YO{I) , I=1,N1) 
RETURN 
C БИСЕР КЕР 22 2а 2,2: 72 e e e e = - m ~ 
FORMAT(715) 


Фар Фо оо ФФ 


N 
© 


r T) = YI(I) 


FORMAT(6F10. 0) 

FORMAT(/,' ORIGINAL COORDINATES',/,' X/C') 
FORMAT(/,' INTERPLATED COORDINATES',/,' X/C') 
FORMAT( ! О, 

FORMAT(6F10. 6) 


SUBROUTINE STAGR(N,STAG,XO,YO,XSTGR,YSTGR) 
MENS TON XOCI00), YOCI00),XSTGR(CIOO0),YSTGR( 100) ,DS( 100) 


ХОТЕ = 0.5 * (XO(1)*XO(N)) 

YOTE = 0.5 * (YO(1)+YO(N)) 

DS€1) = SQRT((X0(1)-XOTE)**2 + (YO(1)-YOTE)**2) 
DSM = DS(1) 

pom T = 2 .N 

DS(I) = SQRT((XO(1)-XOTE)**2 + (YO(1)-YOTE)**2) 
E SC. LT. DSM) GOTO 10 

ВТ 

DSM = DS(I) 

CONTINUE 


ҮҮҮ = ҮОТЕ-ҮО(ІМ) 

ХХХ = XOTE-XO(IM) 

IF (YYY .EQ. 0.0 .AND. XXX .EQ. 0.0) THEN 
ANG = 0.0 

ELSE 

ANG = ATAN2(Y"7,XXX) 

END IF 

ANG = ANG + STAG 


COS( ANG) 

SINAN = SINC ANG) 

DOSMS I — 1, N 

YY = YO(1)-YOCIM) 

ХХ = ХО(І)-ХО(ІМ) 

IF (YY .EQ. 0.0 .AND. XX .EQ. 0.0) THEN 


COSAN 


я 


ANGCO = 0.0 

ELSE 

ANGCO = ATAN2(YY,XX) 
BND ТЕ 


XSTGRC(I)= XOC1I)*COSAN + YOCI)*SINAN 
YSTGR(I)2 YOCI)*COSAN - XOCI)*SINAN 
CONTINUE 


ме 


INT29510 
51029520 
| 1255 0 
INT29540 
[59029550 
22560 
9579 
INT29580 
ШАЙ 9540 
ENG 25600 
INT29610 
INT29620 
INTZ29630 
INT29640 
#129650 
INT29660 
INT29670 
INT29680 
1070229690 
DN 200 
[81295710 
0297720 
INT29730 
INT29740 
N29 750 
NT 297 60 
L [29770 
INT29780 
ENT 297390 
INT29800 
П 29010 
029620 
INT29830 
INT29840 
INT29850 
INT29860 
NAS а 
INT29880 
ENEZ 2890 
тр 2129900 
INT29910 
[NS 129920 
29950 
INT29940 
о 
|. 2006) 
INT29970 
IN IS 
22990 
INT30000 
INT30010 
INT30020 
TINTS 0020 
INT30040 
INT30050 
INT30060 


RETURN 
END 


118 





APRENDES DS? CASE ADE 


A. EXPERIMENTAL RESULTS 

The experimental results of the C4 cascade were obtained directly from professor 
G.J. Walker, University of Tasmania, Tasmania, Australia, who performed these exper- 
iments. 

The results of the boundary laver measurements of the C4 cascade are given below 
at four inlet angles: 34.1%, 36.3%, 45.6%, and 47.7”. The Reynold numbers, based on the 
chord and the upstream velocity, are 200000, 191000, 173000 and 171000 respectively. 
The results given in the following tables include the displacement thickness (0°), the 


shape factor (H) and the local free stream velocity (UE). 


Table 1. EXPERIMENTAL RESULTS AT INLET ANGLE OF 34,1" 


ó [07 FT] 











UE [FT/SEC] 








Table 2. EXPERIMENTAL RESULTS AT INLET ANGLE OF 36.3° 
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Table 3. EXPERIMENTAL RESULTS AT INLET ANGLE OF 45.67 


ARAS 
27.69 1.92 109.26 













Table 4. EXPERIMENTAL RESULTS AT INLET ANGLE OF 47.7? 












X c ó [107 FT] H LE [FT/SEC] 


сте st [|| ac [me 7 


Е 


2 





The results of the measurements of the velocity profiles in the boundary laver at two 


inlet angles. 34.1% and 36.5* at 509» chord are given below. 


Table 5. VELOCIT* PROFILES AT 50% CHORD. 


0.0 0.0 ЕБИ C Å 
2: 
57 





13.4 na ет 0.933 


1.000 


1.000 


B. C4 CASCADE COORDINATES 


DIMENSION X(0: 100), XU(0: 100) ,XL(0: 100) , YU(0: 100), YL(0: 100) C4 00010 
рАТА А1,А2,А3,А4/0. 15492,0. 06563,0. 2528,0. 2811/ C4 00020 
DATA B1,B2,B3,B4/0. 03866,0. 07871,0. 1467,0. 03448/ C4 00030 

PI = ACOS(-1. 0) C4 00040 

С READ (5,800) NMAX C4 00050 
800 FORMAT (15) C4 00060 
NM£X=33 C4 00070 

Ç READ (5,810) (X(1),I=0,NMAX) C4 00080 
810 FORMAT (6F10. 6) C4 00090 
DO 50 I=0,NMAX C4 00100 
X(I) = (1. 0-COS(PI*I/NMAX))/2. C4 00110 

50 CONTINUE C4 00120 
ро 100 I=0,NMAX C4 00130 

ЗЕТ = SQRT((O. 5/SIN(PI/12))***2-(0. 5-X(I) )***2) C4 00140 

УС (0 5/ТАМ(РІ/12) + SRT C4 00150 

DY = ATAN((0. 5-X(1I))/SRT) C4 00160 

IF (X(I). LT. 0.3) THEN C4 00170 

YT = Al*SORT(X(1)) - A2*X(I) - AS3*X(I)**2 4 A4*X(I)**3 C4 00180 

ELSE C4 00190 
AP Вот) - B372K(1)%%2 + BAZX(I)2#23 C4 00200 

END IF C4 00210 


ет) МС ОБС П) >и 


YL(I) = YC - COS(DY)*YT 
XU(I) s X(I) - SIN(DY)*YT 
XL(I) 2 X(I) + SIN(DY)*YT 
CONTINUE 


WRITE (6,900) (I,XCD ,XUC D) , YU( D), XLCD) , YLC I) , I0, NMAX) 
FORMAT (15,4X,F10. 6,4X,2F10. 6,4X,2F10. 6) 

WRITE (1,910) (XL(1), I=NMAX,0,-1),(XU(I) , I=1,NMAX) 
WRITE (1,910) (YL(1),I=NMAX,0,-1),(YUCI), I=1,NMAX) 
FORMAT (6F10. 6) 

STOP 

END 


ІР? 


0022 
002% 
002: 
0023 
0026 
0027 
002: 
002€ 
003€ 
003] 
0032 
0038 
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